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Abstract 



(N . The stability or lack thereof of nonrelativistic fermionic systems to interactions is studied 

^ ! within the Renormalization Group (RG) framework, in close analogy with the study of 
O 

o 

cn 



critical phenomena using 0^ scalar field theory. A brief introduction to 0^ theory in four 
dimensions and the path integral formulation for fermions is given before turning to the 
problem at hand. As for the latter, the following procedure is used. First, the modes on 
either side of the Fermi surface within a cut-off A are chosen for study in analogy with 



■ the modes near the origin in 0^ theory and a path integral is written to describe them. 

O ■ Next, an RG transformation which eliminates a part of these modes, but preserves the 

action of the noninteracting system is identified. Finally the possible perturbations of this 

•i-H . 

^ ' free- field fixed point are classified as relevant, irrelevant or marginal. A d = 1 warmup 
calculation involving a system of fermions shows how , in contrast to mean-field theory, which 
predicts a charge density wave for arbitrarily weak repulsion, and superconductivity for arb 
itrarily weak attraction, the renormalization group approach correct ly yields a scale invariant 
system (Luttinger liquid) by taking int o account both instabilities. Application of the 
renormalization gr oup ind = 2 and 3, for rotationally invariant Fermi surfaces, automatically 
leads to Landau's Fermi liquid theory, which appears as a fixed point characterized by an 
effective mass an d a Landau function F, with the only relevant perturbations being of 
the superconducting (BCS) type. The functional flow equations for the BCS couphngs are 
derived an d separated into an infinite number of flows, one for each angular momentum. It 
is shown that similar results hold for rotationally non-invari ant (but time-reversal invariant) 
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Fermi surfaces also, with obviou s loss of rotational invariance in the parametrization of the 
fixed point interactions. 

A study of a nested Fermi surface shows an additional relevant flow leading to charge 
density wave formation. It is pointed out that for small A/Kp, a 1/N expansion emerges, 
with N — Kp/A, which explains why one is able to solve the narrow cut-off theory. The 
search for non-Fermi liquids in d — 2 using the RG is discussed. Bringing a variety of 
phenomena (Landau Theory, charge density waves, BCS instability, nesting etc.) under the 
one unifying principle of the RG not only allows us to better understand and unify them 
, but also paves the way for generalizations and extensions. The article is pedagogical in 
nature and is expected to be accessible to any serious graduate student. On the other hand 
its survey of the vast literature is mostly limited to the RG approach. 
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I INTRODUCTION 

This article is an expanded version of a short paper (Shankar 1991) in which the apphcation 
of the renormahzation group ( RG ) methods to interacting nonrelativistic fermions in more 
than one spatial dimension was considered. It contains more technical details than its pre- 
decessor and is much more pedagogical in tone. Several related topics are reviewed here so 
that readers with a variety of backgrounds may find the article accessible and self-contained. 
Consequently each reader is likely to run into some familiar topics . When this happens he 
should go through the section quickly to ensure that this is indeed the case and get used 
to the notation and conventions. Upon reading this article, the readers with a condensed 
matter background will see how the RG allows us to synthesize a variety of seemingly un- 
related phenomena in condensed matter theory such as Landau's Fermi Liquid Theory, the 
BCS instability. Charge Density Wave and Spin Density Wave instabilities, ne sting and so 
on. Readers familiar with the RG but not these topics , will see that by following a route 
parallel to the one that led to a very successful treatment of critical phenomena, we are 
automatically le d to many known results in the above mentioned topics and newer ways of 
understanding them. However there are also many fascinating differences between critical 
phenomena and the phenomena considered here which make it very interesting from the 
point view of the RG . At the time of wr iting, there are relatively few new results and 
the emphasis is on the deeper understanding and unification of the older results the RG 
aff ords us. However the machinery developed here, especially for anisotrop ic systems, has 
the potential for changing this state of affairs in the not too distant future. The author is 
working on a few new applicatio ns and hopes the readers will find many more. 

The concept of the RG was first introduced by Stuckelberg and Petermann (1953). Its 
implications for quantum electrodynamics were explored in a sem inal paper by Gell-Mann 
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and Low (1954). These concepts were extended and ge neralized by Callan (1970) and 
Symanzik (1970). 

What is the RG ? When we speak of a group in quantum mechanics we are thin king 
of symmetry operations, i.e., transformations that leave the physic s invariant. What is 
the transformation here? Let us consider quantum electrodynamics. When we compute a 
physical quantity like the scattering ra te between electrons, in a power series in the coupling 
constant a, we find that the coefficients of the series are given by integrals over par tide 
momenta k, and that these in turn diverge because the allowed value s of /c go up to infinity. 
These ultraviolet divergences are at vari ance with experiment which gives finite answers for 
all physical quantities . Renormalization is the way to reconcile these two facts. In this 
scheme, one first cuts off all integrals at the cut-off A. This gives answers that are finite, but 
dependent on A, which is an artifact in continuum theory. To get around this, one asks if it 
is possible to choo se for each cut-off A a corresponding couphng q;(A) so that the physical 
quantities like scattering amplitudes come out A- independent. (In quantum electrodynamics 
one must also renormahze the mass of the electron with the cut-off. The word couphng shall 
mean all s uch parameters that define the theory.) It is by no means obvious that this can 
be done in every field theory. However in the case of quantum electrodynam ics or any 
renormalizable field theory, one can prove that to any giv en order in perturbation theory, it 
is possible to choose a handful of pa rameters of the model in a cut-off dependent way so as to 
make physics at momenta mu ch smaller than the cut-off independent of it. Since the cut-off 
is eventua Uy sent to infinity, this means physics at any finite momentum. This change in 
the cut-off by a factor s, accompanied by a suitable change in couplings is an invariance of 
the theory. Th ese transformations form a group with the composition rule that a ch ange 
by a factor si followed by a change by a factor S2 should e qual a change by a factor siS2- If 
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we write 

(1) 

so that 

m = Aoe-* (2) 

where Aq is some fixed number, the group composition law is that when two transformations 
are implemented in sequence the param eters t add. 

A central quantity in this approach is the (5-function defi ned as follows: 

m - f (3) 

where g is the generic name for the coupling constant(s). Our convention is the one used in 
condensed matter physics wherein i ncreasing t decreases the cut-ofi^. The field theorists use 
the opposite convention and differ by a sign. To avoid all confusion let us consider the case 
of Yang-Mills theory in which 

^ = cg^ + higher orders c > 0. (4) 

If we integrate this equation from t — Q to t — t, (so that the cut-off changes from Aq to 
AqC"*) we find 

' ^ l-2gmcf 

What this equation tells us is that as we send the cut-off to infinity {t to — oo), we must 
reduce the coupling to zero logarithmically: 

g\t):^l/\t\. (6) 

The point to notice in all of the above is that one is interested in how to vary the cut- 
off only with the intention of eventually sending it to infinity, which is where it belongs in 



a continuum theory. Had the theory been free of ultraviolet divergences, the question of 
changing the couplings and the cut-off, keeping the physics invariant, may never have come 
up. It is clear that from this vantage point the RG has no place in condensed matter physics 
where the degrees of freedom live on lattice and there is a natural cut-off on all momenta: 
A ~ 1/a, where a is the lattice constant. 

This point of view was dramatically altered following the work of Kadanoff (Kadanoff 
1965) and Wilson(1971) who gave a different and more physical interpretation of renormal- 
ization. In this modern view one contemplates changing the cut-off (and the couplings) even 
in a problem where nature provides a natural cut-off such as the inverse lattice spacing and 
there are no ultraviolet infiniti es. We will now discuss an example from statistical mechanics 
wh ere the value of such a procedure is apparent. The discussion will be schematic since a 
more detailed one follows in the next section. 

Let us consider a cubic lattice (in d dimensions) with a real scalar field 0(n) at each site 
labeled by the vector n with integer coefficients. The classical statistical mechanics of this 
system is described by the partition function 



This is just the usual sum over configurations with the Boltzmann factor e~^^ written in terms 
of the action S, also called the hamiltonian. (Both terms will be used inter changeably to 
prepare the reader for what happens all the time in the literature.) As long as the number of 
sites is finite, S is just a r egular function and Z is just a multiple integral. In the limit of i 
nfinite sites S becomes a functional and Z becomes a function al integral or a Feynman path 
integral. Feynman introduced his path integral to describe quantum mechanical problems in 
d spatial d imensions as a sum over classical configurations in d + 1 dimensions. Thus our 
Z could very well stand for Feynman's representation of a quantum problem in one lower 




(7) 



n 
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dimension and the following consideration apply to it. [] For the problems with bosonic 
operator s the derivation of the path integral can be found in Section IIIB of the review by 
Kogut ( Kogut 1979). The derivation for the fermionic problem will be given in section III. 

A typical quantity one is interested in is the average of the correlation between the 
variables at two different sites, also called the two-point function: 

G{ni,n2) = Glui — 122) [assuming translation invariance) (8) 
= <0(ni)0(n2)> (9) 

/nnt^0(n)e^W")) ■ ^ ^ 

For long separations this correlation function typically falls off exponentially as 

G{ni - 712) - e-l"i-"2l/« (11) 

where ^ is the correlation length. The exception is when the parameters are such that the 
system is at a critical point , as in the case of a magnet undergoing a Curie transition from 
the ferromagnetic to paramagnetic state. In this case it falls like a power 

G{n, - n2) ^ (12) 

where x is a critical exponent. Other critical exponents characterize other power laws at the 
critical point. A remarkable feature that we will address shortly, is that several systems with 
microscopically distinct hamiltonians (or actions) have the same critical exponent. 

In the case of quantum problems written as path integrals the correlation length is related 
to m, the mass gap, or the lowest excitation energy above the ground state as per 

e = 1/m (13) 

^The d+l-th dimension is of course time. One works with imaginary time, with the option to analytically 
continue to real time at the end if needed. In this case one often finds that all d+ 1 dimensions are equivalent. 
In this discussion of the scalar fiel d we will assume that this so. 

^In the fermion problem we are going to study, the (imaginary) time and space directions are not equiv- 
alent. The correlation length in these discussions refers to the time direction. 



and the critical case corresponds to m = 0. 

An equally complete description of the above system is possible in terms of the Fourier 
transforms: 

m-^T.e''Mn), (14) 

where V is the volume of the system. The allowed momenta k lie within a Brillouin cube of 
sides 27r/a in all directions. The partition function becomes 

Z = J n d(t>{k)e^^^^^'^\ (15) 

|fe|<7r/a 

and G{k), the Fourier transform of G{ni — 712) is given by 

< ct>{k^)(t>{k2) > = {2T:Y5^''\k^ + k2)G{ki) (16) 

/ng#(fc)0(fci)0(fe2)e^W^^» 
/nfc#(^)e^(^(^)) 

Let us now imagine that we are only interested in the physics at long distances, (compared 
to the lattice spacing a), for example in G{r) for large separations r. In momentum space 

— * 

this translates to small k. To be specific let us say we are interested only in correlations of 
modes that lie within a tiny ball of size A/s (with s very large) centered at the origin, and 
not interested in the modes that lie in the rest of the Brillouin zone, taken to be a sphere 
(rather than a cube) of radius 27r/a. (This modification makes no difference to the small k 
asymptotics.) We will refer to the small k modes as "slow modes" and the others as "fast 
modes" . Let us define two sets of variables: 

0< = (t){k) for < k < A/s {slow modes) (18) 
0> = 0(A;) for A/s <k<A {fast modes). (19) 

By assumption we are going to be interested only in correlations of 0<. However at present 
we are computing these objects by doing an integral over fast and slow modes. The first 
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step in the RG program is to ask if there is an effective action or Boltzmann weight e^'^'^<^ 
such that when integrated over just the slow modes, it will reproduce all the slow correlation 
functions. We shall now see that the answer is affirmative. 
Let the action be expressed as follows: 

5(0<, 0>) = 5o(0<) + 5o(0>) + Si{cj>^, 0>) (20) 

where 5*0 is a quadratic function of its arguments that separates into slow and fast pieces 
and Si, called the interaction, is the part which mixes the two. Then 



Z ^ j W d(t){k) n ci0(A;)e^°('^<)e^°('^>)e^^('^<'*>) (21) 

Q<k<K/s A/s<k<A 

= j[d4)<] j[d4)y\e^^^^<^e^^^^>^e^'^^<''^>^ (22) 

= J[d(P<]e^°^^<^ y[d0>]e^°('^>)e^^('^<'^>) (23) 

^ /[d0<]e^'(^<) (24) 

which defines the effective action 5''(0<). Let us manipulate its definition a little: 



^ eSo(0<) ^ g5, (</.<,</->) (^25) 

where <>o> denotes averages with respect to the fast modes with action 5*0 and where the 
corresponding partition function Zo> has been dropped in going to the last line, since it will 
merely add a constant to the effective action independent of 0<, which in turn will make no 
difference to any correlation function of slow modes. 

Although 5"(0<) provides a good description of the slow mode physics, the RG trans- 
formation has two more steps besides the above mode elimination. These steps will now be 
motivated. 
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One aim of the RG is to see the various parameters in the interaction evolve or flow as 
the cut-off is reduced, i.e., to compute the /9-function. Suppose before mode ehmination we 
had 

5(0) ^ + + . . . (26) 

and after, 

5'(0<) = ry< + u'ct>'^ + . . . (27) 
(The action above is schematic. For example u(j)^ could be the shorthand for 

j dkidk2dk^dki5{ki + k^ + k2 + ki)u{ki, ks, k2, /ci)0(/c4)0(^3)0(^2)0(^i) (28) 

where u{k4, ...ki) is a coupling function and not just a coupling constant.) In any event, we 
are trying to compare r to r', u to u' and so on. The problem with doing that is that we are 
comparing apples to oranges. The old and new theory are defined on two different kinematical 
regions. For example the couphng m(A, A,A,A) has no counterpart in the effective theory 
which has all its momenta below A/s. (In field theory where the old and new cut-off are 
both sent to infinity, this point does not come up.) To remedy this defect, we will define 
new momenta after mode elimination: 

k' = sk (29) 

which run over the same range as k did before elimination. 
There is just one more problem. Consider two actions: 

5(0) = r^^ + T/c/.^ (30) 
5'(0) = 4r0^ + 16ii0^ (31) 

which seem different. They are however physically equivalent because we can simply define 
2(f) — (f) in the second action (and ignore the Jacobian in the functional integral since it is a 
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0-independent constant) and reduce it to the first action. In other words, certain changes in 
parameters are not of physical importance since they can be absorbed by field rescaling. To 
weed these out, we will follow mode elimination and momentum rescaling by a field rescaling, 
defining new fields: 

= rV<(^'A) (32) 

and choose C such that a certain coupling in the quadratic part of the action has a fixed 
coefficient. The final action S' will then be expressed in terms of this new field. Thus the 
three stages in the RG transformation are as follows: 

• Eliminate fast modes, i.e., reduce the cut-off from A to A/s. 

• Introduce rescaled momenta: k' — sk which now go all the way to A. 

• Introduce rescaled fields (f)' {k') = (~^4'^{k' /s) and express the effective action in terms 
of them. This action should have the same coefficient for a certain quadratic term. 

With this definition of the RG transformation , we have a mapping from hamiltonians or 
actions defined in a certain phase space to actions in the same space. Thus if we represent 
the initial action as a point in a coupling constant space, this point will fiow under the RG 
transformation to another point in the same space. This definition of the RG opens up a 
possibility that did not exist without all three steps: a fixed point S* of the group action, that 
is to say the action function which reproduces itself after the three step RG transformation. 
Geometrically this means that the point S* does not move or fiow under the RG action. If 
the system had a correlation length before the RG (in the old units), in the new units (in 
which momenta get boosted by a factor s) it would decrease to ^/s. On the other hand, at 
the fixed point, it must remain the same under the RG . This means that the correlation 
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length at a fixed point must have been either zero or infinite. We shall be dealing with the 
latter case here. Fixed points will dominate our analysis. 

In summary, we see that in the modern viewpoint, the cut-off is not to be viewed as an 
artifact to be sent to infinity but as the dividing line between the modes we are interested 
in and the modes we are not interested in. The preceding discussion explains how we may 
change the cut-off and the couplings without affecting the slow mode physics even in a 
problem where there were no ultraviolet infinities. Let us now understand why we would 
want to do such a thing. 

Consider the remarkable phenomenon of universality. How can systems with different mi- 
croscopic hamiltonians have the same decay exponent x in their critical two-point functions? 
The RG explains this as follows ( Kadanoff 1965,1977; Wilson 1971, Wilson and Kogut 1974, 
Fisher 1974,1983; Wilson 1975 ). Let Sa and Sb be two critical hamiltonians defined in the 
full k space. Each is described by a set of coupling constants. Let us represent each as a 
point in a space S (in the notation of Wilson and Kogut 1974, Section 12) of hamiltonians 
where along each axis we measure one coupling constant. The fact that Sa Sb implies 
that they are given by distinct points in coupling constant space and that there are many 
observables that differ in the two cases. Consider however extreme long distance physics, 
in particular the long distance decay of two-point functions. To calculate these we need 
just the slow modes. To this end let us trade each hamiltonian for its equivalent one after 
renormalization down to a very small cut-off. What we will find is that they both asymp- 
totically approach the same fixed point hamiltonian S* where the flow stops. This explains 
why they share the same long distance physics and in particular the exponent x. Although 
the coupling constant space is infinite dimensional let us consider a toy model in which it is 
three dimensional. Let all critical hamiltonians (in particular Sa and Sb) he in the x — y 
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plane. Under the RG they all flow to S*, which lies at say the point (1, 1, 0). Let us shift the 
origin of coordinates to the fixed point. Any deviation S — S* that lies in the critical plane 
is termed irrelevant in the RG terminology since it renormalizes to zero and hence makes no 
difference to long distance physics. The fixed point, being a special case of a critical point, 
will of course have power law decay of correlations. 

What we see is that if to this fixed point an irrelevant perturbation is added, the perturbed 
system will also have the same power law decay. If the functional integral stands for some 
quantum system wrtten as a path integral, this means that a gapless system will remain 
gapless if an irrelevant perturbation is added. This idea will be invoked later in this article. 
By contrast, any deviation off the critical (x—y) plane is called relevant and will get amplified 
by the RG transformation. The long distance behavior of correlations in this problem is 
unclear; it is controlled by the ultimate destination of this flow, and typically (but not 
always) corresponds to exponential decay. In the general problem there can also be marginal 
perturbations, which neither grow nor decay under the RG transformation . They play a 
major role in the nonrelativistic electron problem to which we now turn our attention. A 
truly marginal perturbation does not cause a gap. 

Table I summarizes some of the above concepts. 

The preceding discussions have set the stage for introducing our main topic. Consider 
a system of noninteracting fermions ^ at zero temperature (T = 0) either in the continuum 
with a dispersion relation 

E = K^/2m (33) 

We use the term fermion instead of simply electron to accomodate spinless fermions which do not exist 
in nature but simplify the analysis by obviating the need for spin indices. While the spin of the electron is 
certainly not ignorable when comparing theory to experiment, it will be seen that it really is an inessential 
complication in the RG program to be described here and maybe incorporated readily. 
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or on a lattice with some energy function E{K) defined within the Brillouin zone. Qln all 
cases, one particle states with E < fi, where fi is the chemical potential, are filled in the 
ground state of the many-body system. The filled states are bounded by the Fermi surface 
. In the continuum in ci = 2 or d = 3 the Fermi surface is a circle or sphere respectively of 
radius 

Kp = ^J2mfi. (34) 

This ground state has gapless excitations corresponding to the promotion of fermions from 
just below the Fermi surface to just above it. The central questions we ask in this paper are 
the following: 

• If some perturbation is added to the free theory, will the system develop a gap at once 
or will it remain gapless? 

• If it remains gapless what is the natural way to describe the low energy physics of 
the system, in particular its response to "soft probes", probes of low frequency uj and 
momentum Q7 

The answer to these questions are clearly dictated by the modes near the Fermi surface, at 
least for the case of weak perturbations. For example in any kind of perturbation theory 
of the ground state, these modes will come with the lowest energy denominators. We will 
therefore focus on modes within a bandwidth A of the Fermi surface (the slow modes of this 
problem) and get rid of all the modes out side this cut-off (the fast modes). In the case of 
fermions in free space we define a lower case momentum 

k = \K\- Kf (35) 

"^We use upper case letters to denote momenta measured from the origin in contrast to the preceding 
discussion where lower case symbols were used. This is a deliberate departure from convention and reflects 
the different physics that emerges here. 
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and work with modes obeying 

\k\ < A. (36) 

In the case of electrons on a lattice, the wave vector \K\ is no longer a measure of energy and 
we must keep those modes that whose energy lies within some cut-off. This complication 
will be discussed in the sections devoted to lattice problems. For the present let us focus 
on electrons in free space and imagine an annulus or shell (in two or three dimensions 
respectively) of thickness 2A with mean radius Kp within which reside the slow modes of 
this problem. 

Let us now turn to the elimination of the fast modes outside the cut-off. This may be 
done within the operator formalism by the use of projection operators to define an effective 
hamiltonian restricted to the subspace of slow modes. This effective quantum hamiltonian 
depends on the cut-off in such a way as to produce cut-off independent results for the 
surviving slow modes and the fixed point, if any, is unaffected by this transformation. This 
is the approach used by Wilson (1975), Anderson and Yuval (1970) , Nozieres (1974) and 
Krishna- Murthy et al (1980) in their treatment of the Kondo problem. However this problem, 
which is a paradigm for how the RG is to be used in quantum problems in many-body 
physics, is essentially one dimensional. By contrast the problems we deal with here are 
truly two and three dimensional and the application of the RG to these has a short history. 
Although Anderson had suggested this possibility even in his book (Anderson 1984) no 
detailed analysis was carried out for some time. Benfatto and Gallavotti (1990) and Feldman 
et al (1990, 1991, 1992) then combined the RG with rigorous bounds to study (to all orders 

in perturbations) the stability of gapless Fermi systems to perturbations. Shankar (1991) 

^Here one deals with a conduction band of electrons interacting with a single fixed impurity. By using 
spherical waves (instead of plane waves) centered on the impurity and keeping just the s-wave, one reduces 
it to a quasi-one dimensional problem in the radial coordinate. 
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developed the method to be described here, which is less rigorous, more intuitive, covers other 
instabilities like charge or spin density waves, rotationally-noninvariant systems and maybe 
easier to use for people with a background in critical phenomena or modern field theory. 
More recently Polchinski (1992) employed a very similar approach to the nonrelativistic 
fermion problem to better understand the concept of effective field theories in particle physics. 
Weinberg (1993) recently derived the effective low energy action and RG flow equations for 
superconductors with Fermi surfaces that obeyed time- reversal symmetry and nothing else. 
All these approaches are fundamentally different in spirit from the method used by Hertz 
(1976) who completely integrated the fermions in favour of some bosonic variables. In 
particular he integrated the modes at the Fermi surface . This is analogous to integrating 
the k — modes in critical phenomena. The effective theory for the bosons then has 
singular parameters. Hertz handled this problem carefully and found a way to analyze phase 
transitions that can be described by the bosonic variables. 

The approach described in this paper is as follows. To heighten the analogy with critical 
phenomena one first shifts from the operator approach to a path integral approach and 
writes down the path integral. First the noninteracting problem is considered. Since it is 
gapless one expects it to be the fixed point of a RG transformation in which the cut-off is 
reduced. Such a transformation is found. With respect to this transformation perturbations 
are classified as relevant, irrelevant or marginal. In the last two cases the system will remain 
gapless while in the first case one can only make statements if one assumes that the behavior 
seen for small perturbations persists at strong coupling also. In all cases considered, this 
corresponds to the opening of a gap in the spectrum. While all this sounds like critical 
phenomena, (and is meant to), there are crucial difference which can be traced back to the 
nature of the phase space for slow modes. In critical phenomena all the long distance physics 
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comes from a tiny ball centered at the origin of k-space and the fixed point is characterized 
by a few couphngs. Q 

The same is true for continuum field theories hke quantum electrodynamics or quantum 
chromodynamics: both fermion and boson momenta are restricted to a sphere of radius A 
centered at the origin. In the problem at hand, we renormalize not towards a single point, 
the origin, but a surface, the Fermi surface (which may itself change under renormalization 
in the nonspherical case.) Unlike in critical phenomena where all momenta and momentum 
transfers are small (bounded by the cut-off) here only k = \K\ — Kp is small and large 
momentum transfers of the order of Kp are possible within the slow modes. Renormalization 
only reduces the dimension normal to the Fermi surface, the tangential part survives. As for 
the fixed point, it is characterized by a surface and coupling functions defined on it. Q Notice 
that (i = 1 is special: here the Fermi surface is a set of two disjoint points. Apart from 
this doubling (which converts nonrelativistic fermions into Dirac fermions) we have the same 
situation as in a continuum field theory in one space dimension and there are once again 
just a few coupling constants. This is why there has been a lot of activity and a lot success 
(Bourbonnais and Caron 1991, Solyom 1979) in applying the RG to one dimensional fermion 
problems in condensed matter and a lot of resistance to going to higher dimensions. 

We now turn to the details. In Section II the reader is given a very brief review of how the 

RG works for a scalar field theory in four dimensions. This will serve to remind the readers 

thank Pierre Hohenberg for pointing out to me an exception: a problem not involving a Fermi surface, 
which has nonetheless a similar phase space after any amount of renormalization: the condensation of a 
liquid into a nonuniform state, studied by Brazovskii (1975). See Swift and Hohenberg (1977) for the study 
of fluctuations on an equivalent model. It is an open question whether the methods developed here can be 
applied to Brezovskii's problem. 

^AU this can be stated in another way. In field theories or in critical phenomena one also runs into 
coupling functions. But these are functions of just k. When Taylor expanded in fc, only a few terms are 
marginal or relevant. In the present problem, the coupling functions depend on k as well as the coordinates 
of the limiting Fermi surface . The latter never get eliminated and all terms in the Taylor series for the latter 
will be important. This point will be discussed further as we go along. 
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familiar with the subject the highhghts that we will recall frequently in our progress by 
analogy. As for the newcomers, it will give them the minimum required to follow this article. 
References for more details will be given. Section III explains how a path integral can be 
written for fermions and how one is to extract correlation functions from it. This will require 
the introduction of Grassmann variables. Readers not used to these should kill two birds 
with one stone by using the pedagogical review provided here to learn this tool which is often 
used in condensed matter theory. In Section IV we study the problem of spinless fermions 
in one dimension at half-filling: with one particle per every other site on the average. This 
section serves as a warm up for the RG program, since, as explained above, it resembles the 
run of the mill field theory in one dimension. It also shows the power of the RG: whereas 
mean field theory ( a self-consistent approximation to be detailed later) predicts a gap for 
the smallest repulsion, and superconductivity for the smallest attraction, the exact solution 
tells us that the system remains gapless for a finite range of coupling of either sign. It will 
be seen that the RG gives results in agreement with the exact solution. Sections V and VI 
deal with circular and spherical Fermi surfaces. To lowest order in a perturbative expansion 
(Section V) one finds that there exists a fixed point described by two marginal coupling 
function F and V which depend on the angles on the circle or sphere as the case may be. 
To the next order (Section VI) one finds that F is still marginal while each coefficient in 
the angular momentum expansion of V grows to produce the superconducting instability 
if attractive and renormalizes to downwards if repulsive, a result originally discovered by 
Morel and Anderson (1962). No new surprises come at higher orders. This is explained 
in the next section. The fixed point theory which exists in the absence of V is what is 
known as Landau's Fermi Liquid Theory. The Kohn-Luttinger effect which destroys the 
Fermi liquid at low temperatures is derived in the RG language. Section VII provides a new 
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way of understanding why it is possible to solve the fixed point theory characterized by the 
interaction F even though F is not necessarily small. This is tied to the fact that certain 
theories with a large number of fields can be described by an expansion in 1/N, N being the 
number of components. (In other words the coupling need not be small as long as 1/N is.) It 
is shown that the Fermi system with cut-off A has a 1/N expansion with = Kp/A.. Thus a 
given problem in the full momentum space can initially be renormalized to a small A theory 
(without running into any singularities) and then when N is large enough, solved in the 
1/N approximation. Section VIII has a discussion of Landau's Fermi Liquid Theory. Only 
some aspects of this extensive field are brought up. In Section IX we consider non-circular 
Fermi surfaces with no special features other than time reversal invariance: if K lies on it so 

— * 

does —K. It is found that one must stop using \K\ as a measure of energy and use actual 
equal-energy contours to define the fast and slow modes. The net result is exactly as in the 
rotationally invariant case except for the fact that F and V now depend on more variables 
due to the lack of rotational invariance. Section X deals with the very interesting case of 

— * — * — * 

nested Fermi surfaces in d — 2: surfaces such that if K lies on them, so does K + Qn where 
Qjv is a fixed nesting momentum. We choose to illustrate the ideas with spinless fermions on 
a rectangular lattice in which case the nesting vector Q^v has components (tt, tt). (Readers 
unfamiliar with nesting may wish to peek at Fig. 17 in Section X for an example.) What we 
find is that to lowest order a third couphng function W insinuates itself at the fixed point. At 
next order it begins to fiow. One can show that there are definitely some relevant directions 
if this force is repulsive and these tend to produce charge-density waves: the ground state 
has a nonuniform charge density which oscillates with momentum Q^. In Section XI we use 
the methods developed here to look for non-Fermi liquids in two dimensions. Regrettably 
the results are negative for the case of weakly coupled problems with a circular Fermi surface 
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. Section XII contains the summary and outlook. Many of the remarks made in this preview 
will be repeated there, and the reader will have a clearer picture of their significance. The 
Appendix deals with two special topics: Coulomb screening and the Kohn-Luttinger effect 
as they appear within the RG framework. 



The problem chosen to illustrate the RG at work involves a complex scalar field in d — A. The 
functional integral can be viewed either as describing the quantum field theory of a charged 
scalar field in three space dimensions or as describing the classical statistical mechanics of 
a system with one complex field or a two real fields at each point on a lattice. The Ising 
model which is described by a real field is not chosen here since the fermion problem we will 
study later involves charged Fermi fields. Readers new to the problem should be aware that 
this section has the very limited objective of making the rest of the paper comprehensible 
to them. For a deeper introduction to critical phenomena, the reader is directed to any 
of the excellent reviews (Wilson and Kogut 1974, Fisher 1974, Kadanoff 1977) or books 
(Goldenfeld (1992), Itzykson and Drouffe(1989), Le Bellac(1991), Ma (1976) Plischke and 
Bergersen(1989), Zinn-Justin(1989)). Readers famihar with the subject are still urged to 
skim through this section to get acquainted with the notation and as well as to refresh their 
memory, since our approach to the interacting Fermi problem will rely heavily on analogy 
to this problem where the RG approach has been spectacularly successful. 
The partition function for this problem is 
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AN EXAMPLE OF THE RG FROM d = 4 
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[ci0ci0*]e^('^''^*\ where 
dRe(j){k)dIm(j){k) 



(37) 
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(38) 



|fe|<A 
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5(0,0*) = / 0*(fc)J(fc)0(^)^, and (39) 

J\k\<A [^'^) 

J{k) = 2 [(cos k^-l) + (cos A;j/ - 1) + (cos k^-l) + (cos A;^ - 1)] . (40) 



This action is obtained by Fourier transformation of the following nearest neighbor interac- 
tion in coordinate space: 

S{<t>A*) ^ -\Y.\^{n) - <t>{n + ^^ (41) 

n,i 

where n is the vector with integer coordinates used to label the sides of the hypercubic lattice 
and i is any of the eight unit vectors in the direction of increasing or decreasing coordinates. 
Notice that this action favors the ahgnment of neighboring fields, i.e., is ferromagnetic. 

Since we are interested in small k physics, let us hereafter approximate J{k) by its leading 
term in the Taylor series and write 

-50(0,0*) = - / , ^4>\k)k'm^^ (42) 

J\k\<K (,27rj* 

Z = / [c^0d0*]e-^i<A^*(^^)^^^^^^)(^ (43) 
J\k\<k 

[^A^A*^ TT dRe(t){k)dIm(t){k) ^ dcp* {k)d(t){k) 

[d<t>d<t> ] = 11 = 11 (44) 

|fe|<A ^ |fc|<A "^^^ 

This is called the gaussian model. The corresponding functional integrals is product of ordi- 

— * 

nary gaussian integrals, one for each k. This makes it possible to express all the correlation 
functions in terms of averages involving a single gaussian integral. The only averages that do 
not vanish are products of an even number of variables, wherein each (j){k) is accompanied 
by its complex conjugate. This is because the action and measure are invariant under 

0(fc) ^ 0(fc)e*^ 0*(fc) ^ 0*(fc)e-^^ (45) 

where 9 can be difi^erent for difi^erent fc's. Thus the only integral we will ever need follows 
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from the simple problem involving just a pair of complex conjugate variables z and z*: 



roo dzdz* ^- 
J-oo 2m 



< zz* >= °1 . = -. (46) 



roo dzdz* ^-azz* 
— oo 

The other two bilinears have zero average: 

< zz >=< z*z* >= (47) 
because the action and measure are invariant under 

z^ze'^ z*^z*e-'^ (48) 

while the bilinears are not. The reader wishing to verify the above results is asked to switch 
to X and y, the real and imaginary parts of the integration variables and to use 

dzdz* dxdy 



2m TT 
If there are two sets of variables we have 



(49) 



poo dzidzl dz2dz* -aizizi-a2Z2Z^ r 
■^ -oo 2-Ki 2-Ki ^i^j'^ _ "jj 

.oo dzid^dz2dzl aizizl-a2Z2Z* 
J -co 2m 2m 



^ ^*^. J-oo 2ni 2m — l!:^ 1 ^ > (r.r\\ 

^ ^i^J roo dzidzl dz2dz* ~ „ — ^ ''J " ^"^^^ 



As for the four point function, the reader may verify that 

< z*Zjzlzi >=< ij >< kl >< il >< kj > (51) 

This result makes sense: it demands that for the answer to be nonzero, the fields must come 
in complex conjugate pairs. Since this can happen in two ways, the result is a sum of two 
terms. The generlization to more variables and longer strings is obvious. 

In view of the above, the reader will not be surprised that the two-point function in our 
gaussian model is 

< ,.C,)m > = (52) 

= {2Tr)^S\k - k')G{k) (53) 
= < 21 > . (54) 
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and likewise 



< 0*(A;4)0*(4)0(^2)0(A?i) >=< 42 >< 31 > + < 41 >< 32 > . 



(55) 



This is a case of Wick'sTheorem for bosons. For the case of 2n fields, the answer is a sum 
over all possible pairings, each term in the sum being a product of n 2-point functions. The 
result follows from the preceding discussion, upon making the change from Kronecker deltas 



gaussian model is an integral (over k) rather than sum over variable labels. 

Note that G has power law behavior in momentum space ( l/k"^) and hence will do so in 
coordinate space ( 1/r^). Thus the action of the gaussian model is critical and must flow to 
a fixed point under the action of the RG. We will now see that it is itself a fixed point. 

In the first stage of the RG transformation, we integrate out 0>. Since Sj = here, we 
see from Eqn. (P^D that 



to Dirac delta functions in Eqns.(^T]) to take into account the fact that the action in the 




(56) 



We now carry out the last two steps by rewriting the action in terms of 



0'(^') = rv<(^'A) 



(57) 



and obtain 




(58) 



(59) 




(60) 



If we now make the choice 




(61) 
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we find that the gaussian action is the fixed point: 

S'i<P') = Si<P) = S*. (62) 

Having found the fixed point, we next classify its perturbations as relevant, irrelevant or 
marginal. We will only consider perturbations involving an even number of fields. Let us 
start with the quadratic case: 

SS = -[ <P*{k)r{k)m^^ (63) 

J\k\<A l^'^j 

where the coupling function r is assumed to have a Taylor expansion 

r{k) = ro + r2k'^ + ... (64) 

which refiects the short range nature of the perturbation in coordinate space. One often 
writes 

ro = ml (65) 

and refers to ml as the mass term since in the quantum field theory interpretation of the 
functional integral, adding this term to the gaussian model yields the quantum field theory 
of a particle of mass mo- 

Since this perturbation does not mix slow and fast modes, all we have to do is replace (f) 
by 0< everywhere and reexpress the result in terms of new momenta and fields. This gives 

6S'{<p'{k')) = -s' 1^ J*{k')T{k'/s)4>'{k')^^ (66) 

J\k\'<K \^^) 

^ - [ /*{ky{k')<i>'{k')-^^ (67) 

J\k\'<A (^"^J 

where in the last equation we have invoked the definition of the renormalized coupling r'{k'). 
By comparison, we find 

r\k') = s\{k'/s) (68) 
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which imphes that the Taylor coefficients obey 



r. 








(69) 



(70) 




(71) 



and so on. Thus we find that tq is relevant, r2 is marginal and the rest are irrelevant. This 
is a concrete example of how in the low energy physics the coupling function r{k) reduces to 
a few coupling constants. (In fact r2 makes no difference since it can be absorbed by field 
rescaling.) In quantum field theory, where we send the cut-off to infinity, all momenta are 
small compared to the cut-off and the theory is defined by a few coupling constants. We shall 
see that the same thing will happen for the quartic interaction: a coupling function of four 
different momenta will reduce to a single coupling constant. We may understand all this as 
follows. In the original Brillouin zone, of size 1/a, all these functions are nontrivial and we 
need them in their entirety. As we eliminate modes, we need their behavior in a smaller and 
smaller ball near the origin, see Eqn.(^S|). Not surprisingly, the function is well described by 
a few terms in the Taylor series. This is the picture in fixed or "laboratory units" . In the RG 
one uses sliding units that constantly change to keep the cut-off (ball size) fixed at A and the 
same phenomenon appears as the rapid shrinkage of higher coefficients in the Taylor series. 
(Of course as we renormalize, we are not just rewriting the original coupling function in new 
units, the function itself changes due to eliminated modes. But it is expected nonetheless to 
be smooth in k. This is one of the points emphasized in the modern RG theory: elimination 
of modes does not introduce new singularites into the couplings. As we shall see, this is 
because the effect of mode elimination may be expressed in terms of integrals which are 
convergent at both ends.) 
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Let us now consider the quartic perturbation 

= - f 0*(4)0*(3)0(2)0(1M4321) (72) 
where the couphng function obeys the symmetry condition 

m(4321) = u(3421) = u(4312). (73) 

In other words the couphng function is invariant under the exchange of the first two or 
last two arguments. Even if we started with a function that did not have this symmetry, 
the invariance of the measure and the rest of the integrand under this symmetry would 
automatically project out the symmetric part and annihilate the antisymmetric part. The 
factorials up front are conventional and are put there to prevent similar factors from arising 
in subsequent calculations. 

The renormalization of the quartic interaction is complicated by the fact that unlike the 
quartic perturbations, it mixes up the slow and fast modes. Thus we have to use the formula 

= e^°+^^' (75) 



Next we invoke the cumulant expansion which relates the mean of the exponential to the 
exponential of means: 

^gH^ ^ g[<n>+i(<n2>-<Q>2)+...]_ (7g) 

The reader may wish to verify the correctness of this expansion to the order shown. Using 
this expansion we find 

5S' =< 5S > +h< SS^>-<SS>^) + ... (77) 
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Since 5S is linear in u, this a weak coupling expansion. It is now clear what has to be 
done. Each term in the series contains some monomials in fast and slow modes. The former 
have to be averaged with respect to the Boltzmann weight 5'o(0>) by the use of Wick's 
theorem. The result of each integration will be to give a monomial in the slow modes. 
When reexpressed in terms of the rescaled fields and momenta, each will renormalize the 
corresponding coupling. In principle the reader has been given information to carry out this 
process. There is however no need to reinvent the wheel. There is a procedure involving 
Feynman diagrams which automates this process. These rules will not be discussed here since 
they may be found, for example in Sections 3-5 of Wilson and Kogut (1974) or in any good 
field theory book (Itzykson and Zuber 1980, Zinn- Justin 1989). Instead we will go over just 
the first term in the series in some detail and comment on some aspects of the second term. 
Readers familiar with Feynman diagrams should note that while these diagrams have the 
same multiplicity and topology as the field theory diagrams, the momenta being integrated 
out are limited to the shell being eliminated, i.e., A/s < k < A. 
The leading term has the form 



The sixteen possible monomials fall into four groups: 

• 8 terms with an odd number of fast fields. 

• 1 term with all fast modes. 

• 1 term with all slow modes. 

• 6 terms with two slow and two fast modes. 

We have no interest in the first two items: the first since it vanishes by symmetry and the 
second since it makes a constant contribution, independent of 0< to the efi'ective action. 




(78) 
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Consider next the third term with all slow modes, distinguished by the fact that it requires 
no integration (or averaging) over fast modes. This is called the tree level term in field theory. 
The tree level term is obtained from the original perturbation by simply setting = 0<. 
Rewriting it in terms of new momenta and fields we find it leads to the following quartic 
renormalized interaction 



SS',,ree = / , <t^'* {k:)<t^'%k',)<P\k',)<P' {k[)u{k'j S, . . . k[ / s) {27:)' S'ik'.+g-g-k',) U — ^ 

2\2\ J\k\<A [2n) 

(79) 

The reader should note that the field rescaling factor s^^ has been exactly cancelled by 
rewriting the delta function and integration measure in terms of new moments. (Note that 
the delta function scales oppositely to the momenta.) 

It is evident that the renormalized four point coupling is given by 

u {k'^, . . . k'l) = u{k'j s,. . .k'l/s). (80) 

Carrying out the Taylor expansion 

u^uo + 0{k) (81) 
we see that the constant term is marginal 

u'q = uo (82) 

and the rest are irrelevant. This is why the scalar field theory in four dimensions is described 
by a coupling constant and not a coupling function. Hereafter we will replace the coupling 
function by the couphng constant. The effect will be irrelevant in the RG sense. 

We now pass from the tree level term to the six terms which have two slow and two fast 
modes in them. Of these, two with 0>0> or their conjugates are zero. The others clearly 
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renormalize the quadratic interaction: 



1 



SS2{<P<) = -^^0 ^^^^^^(0;(4)0*<(3) + 0*<(4)0;(3))(0>(2)0<(1) + 0<(2)0>(1))^ . 

^(83) 

If we now evaluate the averages of the fast modes we will find that all four terms give the 
same contribution (which takes care of the factorials in front) and we end up with 

55;(0<) = -Mo/, , ^,<f><ik)Mk) r (84) 
J\k\<\/.s [2tt)^ Ja/s [2iTpk^ 

SS',i<p\k')) = -uos' 7l;u<^'*(^')0'(^')A^:^(l-3)7^ (85) 

where in the last step we have used the fact that the area of a unit sphere in four dimensions 
is 27r2. 

Eqn. (P^D gives us the change in tq: 

^''■> = fi^'^' - 

Let us agree to measure tq in units of the cut-off squared and drop the from now on. 

Notice that the quartic coupling has renormalized the quadratic coupling. This is more 
the rule than the exception. The quadratic perturbations were special in that they did not 
generate new couplings. In view of this, we must really study the problem in which both 
tq and Mo are present from the outset. This amounts to replacing the propagator l/k"^ by 
1 / (/c^+ro) in Eqn.(P^. However this only modifies the result to higher order in the expansion 
in To and mo. The flow to this order is 



^0 = s\r, + ^^{l-l/s^)) (87) 



Mq = Mq. 



If we take s = 1 + 1, with t infinitesimal, we find the differential equations 



drn Uf) 
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This completes our analysis of the first term in the cumulant expansion. Let us see briefly 
how the above results follow in the diagrammatic approach. First we associate with each 
quartic perturbation 6S a four pronged X as in Fig. la. The incoming arrows correspond to 
and the outgoing ones to 0*. Each prong can stand for a 0< or a 0>. Next we do the 
average over the fast modes. The prongs corresponding to the matching pairs that give a 
nonzero contribution are joined and correspond to the propagator. The diagrams in Fig. I 
tell us what happens. The first one corresponds to all slow modes and there is nothing to 
average, i.e., no lines to join. Figure lb corresponds to the eight terms with an odd number 
of fast lines. These average to zero. Figure Ic describes the case with two fast and two 
slow lines with both sets coming in complex conjugate pairs. The two fast lines are joined 
by the averaging and the line joining them is the propagator and this corresponds to the 
renormalization of the quadratic term as per Eqn.(|84D. This is called the tadpole diagram. 
Finally Fig. Id describes the case where all lines are fast and come in pairs. We now have 
two propagators. We did not consider this above since it is a constant independent of 0<. 

Notice that although all terms are of order uq, they have very different topologies. The 
tree level term has no loops or sum over fast modes. Figure Ic has one loop and Figure Id 
has two loops. Now the correct way to organize the cumulant expansion is by counting loops. 
The reason is best seen in the language of quantum field theory where the action has the 
prefactor 1/h and the number of loops measures the powers of h. In critical phenomena this 
fact becomes very clear when one works in 4 — e dimensions (Wilson and Fisher 1972). One 
finds then that the loop expansion is an expansion in e. The reader who wants to know more 
should consult references given at the beginning of this section. 
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At zero loops, or tree level, the equations are 



r'o = sVo ^ ^ = 2r (91) 
^ du , . 



Equations ( [89| - pOD are halfway between zero and one loops: they are good to one loop for 
ro and to tree level for uq. To be consistent, we must evaluate the flow of mq to one loop 
also, which means going to second order in uq via the next term in the cumulant expansion 
namely 

I [< {6Sy >-< {6S) . 

Here we draw two crosses and do the usual pairing. All diagrams in which no line runs 
from one cross to the other, i.e., the disconnected diagrams may be dropped since they get 
cancelled by — < [bS^ Of the rest the only graphs that affect u are shown in Figure 2 
and correspond to the following analytical expression: 



,, 2r f d^k 1 f d^k 

M(fc4,...,fci) = Uo-Uq[ . ,, h / 



dA (27r)4 k^\k - k'-i/s + fc'i/sP JdA (27r)4 _ g^js + 



2 



ZSgraph ZS'graph 

+ - f ^ 1 (93) 

2 idA (27r)4 p| _ + fc-'/^/s + 
" V ' 

BCSgraph 

Several remarks are in order. First note that even though we started with a constant u = uo, 
the renormalized coupling has acquired momentum dependence. If we expand the renormal- 
ized coupling in a Taylor series, keeping just the lowest term, we will get the renormalized 
Uq. This is what we will do, and ignore the irrelevant higher terms in the series. This in 
turn means that we can set all external momenta to zero. Before so doing, let us look at 
the three one loop diagrams. Since we need to refer them individually many times we need 
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a system of nomenclature. The one used here is by no means standard. Consider the first 
diagram labeled ZS which stands for "zero sound." In this diagram lines labeled 1 and 3 
meet at a vertex. In Fermi liquid theory a graph with the same topology occurs and is very 
important when Q = ki — is small. The physics of the present problem couldn't be more 
different: the lines here stand for bosons and unlike in Fermi liquid theory, the internal loop 
momenta are restricted to lie at the cut-off rather than take all values within the cut-off. In 
the second ZS' diagram lines 1 and 4 meet at a vertex. Usually when Q is small, Q' = ki — k^ 
is large and this diagram is not very important in Fermi liquid theory, and does not have a 
name. However in problems with nesting this diagram can be important if Q' is the nesting 
momentum. The BCS diagram with lines 1 and 2 meeting at a vertex has the topology as 
one that will appear later in our description of the superconducting instability. The reader 
is once again cautioned that the names of these diagrams are based solely on the topology 
and do not generally imply the corresponding physics. 

Readers familiar with Feynman diagrams could have easily written them down. They 
must however pay attention to the symbol /^^ which reminds us that all internal propagator 
momenta (corresponding to integrated fast modes) are summed only over the band being 
eliminated, which we take to be a shell of thickness dA at the cut-off. 

Readers new to the subject are strongly urged to work out the combinatorics and derive 
this result. They will then see why the factorials were included in the definition of the 
perturbation and why an extra factor of 1/2 appears in the BCS diagrams. 

All readers should note that the one loop correction has a minus sign in front of it, 
reflecting the decrease of the interaction strength as we go to the infrared modes. (Although 
the one loop graphs have a positive value, they reduce mq since the latter is defined to occur 
in the action with a negative sign, see Eqn.(|7DD 
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Let us now set all external momenta to zero, since we are interested in just uq. We are 
now assured that if the loop momentum k lies in dA, so does the other momentum which 

— * — * 

either equals k in the ZS and ZS' diagrams, or equals —k in the BCS case. All the integrals 
are now equal and we get 



duo 5ul 



dt 167r2 
where in the last step we have recalled 

\dA\ 



(95) 



A 

and the area of a unit sphere in four dimensions (Stt^). 
To one loop accuracy we have the following flow: 



dro 



= dt (96) 



= 2ro + auo (97) 

dt 

where a and b are positive constants whose precise values we are no longer interested in. 

We shall now analyze these equations. First observe that besides the gaussian fixed point 
at the origin, there are no other points where both derivatives vanish. Next, the equation 
for Uq is readily integrated to give 

uo{t) = ''f] , . (99) 

This means that if we start with a positive coupling Uq{Qi) and renormalize, the effective 
coupling renormalizes to zero as 1/t = l/(ln(Ao/A)). One says uq is marginally irrelevant. 
In the case of bosons a negative Uq is unphysical since the functional integral over fields 
will then diverge for large fields. In some fermion problems one gets the same equation and 
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negative uq is allowed. In that case the coupling is marginally relevant and grows. The 
above equation, derived for weak coupling, will soon have to be abandoned, in contrast to 
the postive uq case, where it gets more and more reliable at larger and larger t. Notice that 
the fate of marginal couplings (unlike relevant or irrelevant couplings) depends on the sign. 

The statement that mq is marginally irrelevant at the gaussian fixed point needs to be 
understood properly. In particular, it does not mean that if we add a small positive Uq to the 
gaussian fixed point, we will renormalize back to the gaussian fixed point. This is because the 
small Mo will generate an ro and that will quickly grow under renormalization. What is true 
is that ultimately Uq will decrease to zero, but tq can be large. In fact to flow to the gaussian 
fixed point, we must start with a particular combinations of tq and Uq which describes the 



critical surface. All this comes out of the equation (|9^ ) for tq which is integrated to give 



roit) 



„2t 



(100) 



Let us consider large t. Typically tq will flow to infinity exponentially fast due to the 
exponential prefactor, unless we choose tq such that the object in brackets vanishes: 



ro(0) 



which, for very small Mo(0) translates to 



ro(0) 



1 + buoiO)t 



auo{0) 



(101) 
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which defines the critical surface (a line in this case) in the tq — uq plane. Any point on this 
approaches the origin as follows: 



Uo{t) ^ a/t 
ro{t) ^ 



2t 



t 1 + buoiO)t 



a 

'2bt 



(103) 
(104) 
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Figure 3 depicts the state of affairs. Table II summarizes the results for the gaussian fixed 
point and its leading perturbations. 

The analysis of couphngs with more powers of the fields is similar. All of them are 
irrelevant even at the tree level and higher loops cannot change that. For example the 
constant part of the sextic (0*0)^ coupling falls like 

A The Field Theory Approach to the /3-function 

We just derived the flows in the modern approach, which is intuitively very appealing and 
consists of integrating out fast modes. We will now rederive the one loop flow of uq the 
old way, where the aim is to banish all cut- off dependence from physical quantities. The 
two approaches will then be compared and contrasted. The reason we even bring up the 
field theory method is that at higher loops it is more tractable than the modern approach. 
In the Appendix we perform two calculations involving interacting fermions for which the 
field theory method proves more convenient. The present discussion will rather succinct and 
readers new to diagrams will have to work that much harder. 

Consider a field theory with two coupling constants: a mass term tq, a quartic coupling 
Uq and a cut-off A. The physical quantity we wish to hold fixed is r(A;4 . . .ki), called the 
irreducible four point vertex or four point function. (Arrows on vectors will be suppressed.) 
It is defined as follows. Let us define the action of a massive free field Sq 



Consider now (suppressing the momentum integration measure for variables labeled 5 — 8 in 
the quartic coupling) 




(105) 
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((-(/.*(fc4)0*(fc3)0(fc2)0(fcl))e-^./A^*(*^8)^*(fc7)^(fe6)</.(fe5)^^ 
^ Xv0*(fc8)^*(fe7)^(fe6)<^{fcB)^ 

where in going to the last equation, we have multiphed and divided by the partition function 
with action 5*0 and <>o stands for averages with respect to this measure. Notice that aU 
momentum integrals go from to A. This is because we are not eliminating modes, we are 
carrying out a calculation of some correlation function in a given theory. We now calculate 
the answer in a power series in uq by expanding the exponential. We then throw out all 
disconnected diagrams (diagrams in which some lines are not connected to the others) and 
delete the four propagators that hnk the external fields (whose momenta are labeled 1 — 4) 
to the vertices that come from the exponential, and the 5-function for overall momentum 
conservation. This defines r(4321) which is the object we want to be cut-off independent. In 
field theory r(4321) is the scattering amplitude for the process in which 1 + 2 — 3 + 4, and 
is a measure of the interaction between particles. In the gaussian model it will vanish since 
all diagrams will be disconnected, the disconnected diagrams describing the independent 
propagation of noninter acting particles. 
Let us now calculate T to order uq. 

If we expand the exponential in the numerator to first order, we get a connected piece 
in which the external fields numbered 1 to 4 get paired with the quartic interaction fields 
numbered 5 to 8. The factorials get neutralized by the number of ways to pair and the 
propagators get dropped and the net result is that to this order 

r(4321) = Uq. (107) 

(The denominator is set equal to unity since expanding it to order uq will change the answer 
to order Uq.) The reader new to this subject is very strongly urged to carry out the steps 
using Wick's theorem and paying attention to the combinatorics. Since to lowest order in 
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perturbation theory r(4321) = m(4321), we will sometimes refer to u as the amplitude for 
scattering. To be exact it is just a coupling constant in the theory which equals F in the 
weak coupling limit. It does however have the all (anti)symmetries of F under the exchange 
of momentum labels of the external (fermions) bosons. 

The above answer for F is clearly cut-off independent and one may choose uq to match 
scattering experiments. Notice that it is also independent of external momenta. 

To next order we must expand the numerator to order Uq and the denominator to order 
Uo (since F starts out at order uq in the numerator.) The Feynman diagrams are exactly as 
before except for the fact that all loop momenta go up to the cut-off. The result is 



F(4321) = uo-% 
+ 



A d^k 1 /-A d'^k 1 

+ 



(27r)4(F + ro)(|A: + /c3- fciP + ro) Jo (27r)4 (P + ro)(|fc + /C4 - /ciP + tq) 

1 d^k 1 , , 

(108) 



2 Jo (27r)4(P + ro)(| - A; + fc2 + A;i|2 + ro) 
Let us now demand that F(OOOO) be independent of cut-off as the latter goes to infinity. 
In this limit we find 

F(OOOO) =Uo- ul-^ In — . (109) 

Let us now act on both sides with d/dt = —Ad/dA and demand that they vanish. This gives 
us 

from which we find 

Two points need clarification here. First: why did we not take the implicit t derivative of 
the Mq term in Eqn.( p!09| )? The answer is that (3 is of second order in the coupling and 
this will give a third order term. Next one may wonder about F with all /c's not equal to 
zero. Will they also be cut-off independent if we choose mo(A) as above to make F(OOOO) 

39 



cut-off independent? The answer is yes. If we expand the above integrals in the external 
momenta, the integrals will become convergent and cut-of independent if we send the cut-off 
to infinity. Thus the external momenta must be much smaller than the cut-off for the field 
theroy renormalization to work. If we want the physics to be cut-off independent for external 
momenta comparable to A, we will need to introduce new couplings besides uq. ) 

In the same way one can derive the flow for tq be demanding that the pole in the 
full propagator (the two-point function in the theory with uo ^ 0) have a certain cut-off 
independent location. 

Suppose we add just the quartic coupling but no mass term tq to the gaussian model. 
Then we will find that r(OOOO) has an infrared logarithmic divergence. (Send tq to zero in 
Eqn. (p!09D . )This is a physical divergence in a massless theory, analogous the infinite cross 
section for Rutherford scattering in electrodynamics. However the /3-function, which involves 
the derivative with respect to the upper limit of momentum integration is still well defined 
and has the same value quoted above. In the modern approach, even if ro = 0, we will never 
see any infrared divergence in the calculation of the /9-function since the loop momentum will 
now go from A/s to A. This was the meaning of the earlier statement that mode elimination 
does not produce singularities in the parameters that appear in the effective action because 
the flow is given by integrals that are well behaved at both ends. 

Although the two methods gave the same answer, this is a fact that needs some expla- 
nation since the methods are very different. In the modern approach a change in cut-off is 
compensated by a change in an infinite number of couplings while in the latter one tries to 
compensate by changing just ro and uq- How can this be possible? The answer is that in 
the field theory approach one always sends the cut-off to infinity (or equivalently looks at 
correlation functions with external momenta very small compared to the cut-off), while in 
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the modern approach we can ask for quantities involving momenta right up to the cut-off. 
If in the modern approach we hmit ourselves to momenta much smaller than the cut-off, we 
could trade the complicated hamiltonian for a simpler one at low momenta dominated by a 
few marginal and relevant couplings. 

At a graphical level there are differences in the range of integration in the loop graphs 
that contribute to the flow. In the modern approach we demand that each internal line lie 
in the shell of width dA near the cut-off. In the field theory approach, where we take the 
A-derivative of momentum integrals going up to A, the answer is a sum of terms in which 
one of the internal propagator momenta is at the cut-off and the others go up to the cut-off. 
In our flow equation for Uq this difference was suppressed. This was because we argued that 
only the lowest term in the Taylor series for the coupling was marginal and the rest were 
irrelevant, allowing us to set all external momenta equal to zero. This meant that if one line 
in the loop was at the cut-off, the other being either equal to it (ZS and ZS') or opposite to 
it (BCS) also had to be at the cut-off. Thus both lines were at the cut-off in both approaches 
to the flow. Had we been interested in the renormalization of irrelevant operators, we would 
have had to consider non-zero external momenta and the two schemes would have yielded 
different answers. 

Notice that the two schemes do not have to give the same flows, they just have to 
give the same physics (at momenta much smaller than the cut-off). The book-keeping can 
be very different. Consider a more general graph in the field theory approach with many 
internal lines and 4 external lines, so that it contributes to the renormalization of the 4-point 
coupling. Some of these internal momenta may be at the cut-off and the rest below it. Such 
a graph is forbidden in the modern approach. The effect of these graphs (with slow and fast 
momenta in the loops) will appear as follows in the modern approach. First all internal lines 
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(propagators) with slow momenta are snipped and the danghng hnes are made into external 
lines. This graph is then used to renormalize a higher point function with that many more 
external lines, say 6 lines in all if two new external lines were produced by snipping an 
internal line. Suppose we now stop renormalizing and compute an object, say the 4-point 
vertex of slow modes, using these couplings. The answer will be given as integrals over 
slow momenta. The 6-point coupling that was generated by snipping will contribute to the 
4-point function when the two lines that got snipped get joined again. In the field theory 
approach such a contribution would already be sitting inside the effective 4-point coupling 
which got renormalized by graphs with slow and fast lines in the loop. 

Ill PATH INTEGRALS FOR FERMIONS 

In this section the reader is introduced to the path integral representation of fermion prob- 
lems. Some elementary problems involving dynamics and thermodynamics will be first solved 
by operator methods and then the same results will be rederived using path integral methods 
reviewed here. For a more detailed treatment, the reader is asked to turn to the standard 
references (Berezin 1966, Itzykson and Drouffe 1989, Schwinger 1970, Negle and Orland 
(1988)). 

A The fermionic oscillator: dynamics and thermodynamics via 
operators 

Let ^ and be two fermionic operators obeying antzcommutation relations: 

= {*t*^} = 0. (113) 
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Note that the last equation tells us 

^^2^0. (114) 
This equation will be used all the time without explicit warning. The number operator 

N = (115) 

obeys 

= ^t^^t^ ^ ^t(i _ ^t^)^ ^ ^t^ ^ (116) 

Thus the eigenvalues of N can only be or 1. The corresponding normalized eigenstates 
obey 

A^|0> = 0|0> (117) 

N\l> = 1|1>. (118) 

We will now prove that 

*^|0> = |1> (119) 

^|1> = |0>. (120) 

As for the first, 

TV^tjo >= ^t^^t|o >= ^t(^i _ ^t^)|o >= ^t|o >^ (121) 

which shows that ^^|0 > has = 1. Its norm is unity: 

||*t|o > ||2 ^< o|^*t|o >=< o|(l - *t^)|o >=< o|0 >= 1. (122) 

It can be similarly shown that >— |0 > after first verifying that ^|1 > is not a null 
vector, that it has unit norm. 
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There are no other vectors in the Hilbert space: any attempts to produce more states 
are thwarted by \1>^ = = 0. 

Consider now a Fermi oscillator with hamiltonian 

Hq = fio^^^ (123) 

whose eigenvalues are clearly and Qo- 
We will work not with Hq but with 

H = Ho-ijN (124) 

where /i is the chemical potential. For the oscillator, since 

H = {Vlo - fi)^^"^ (125) 

this merely amounts to measuring all energies relative to the chemical potential. ^ 
Let us now turn to thermodynamics. The grand partition function is 

Z = Tre-^(^«-'^^) = e-^^(^'^) (126) 

where the trace is over any complete set of eigenstates, (3 is the inverse temperature 1/T and 
A is the free energy. The latter is clearly a function of fi and jS and its differential is 

dA = - < N >diJ- SdT (127) 

where S is the entropy and < N > stands for the average particle number. Let us verify 
that 

8 A 

<^>=-^ (128) 



®The eigenvalues of H are T = free energies rather than energies. We shall however often refer to H as 
the hamiltonian. 
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as follows: 

_dA _ ^InZ 

- J^^z (1^°) 

= ^ < /?7V >=< N > . (131) 
The free energy A is the double Legendre transform of the internal energy £^(5*, < N >): 

A{i^,p) = E{S,< N >) - ST - ii< N > . (132) 
Thus E must equal < Ho >. This is indeed so: 

_ 9(-lnZ) «91nZ 

= <i/o>- (137) 

where the steps leading to the last hne are left to the reader. It also follows from the 
definition of the Legendre transform that 

dE , , 

so that \i is the (minimum) energy needed to add an extra particle. 

The partition function of the Fermi oscillator is easily found (by doing the trace over 
eigenstates of A^) to be 

Z = 1 + e-'^(^°-'') (139) 

from which it follows that 

A= -iln(l + e-^(^°-'^^), (140) 
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which in turn imphes 



< ^ MO V- (141) 



We shall be interested in the limit /3 — > oo in which case 

< AT >= e{n - Qq) (142) 

which means the fermion is present if its energy is negative (relative to the chemical potential) 
and absent if it is positive. This is to be expected since atT — 0, A — E — ijl<N> and 
minimizing the free energy is the same as minimizing < H >. 

We now consider the dynamics. Prom the Schrodinger operators we can form Heisenberg 
operators: 

^(i) = e^^**(0)e-^^* = *(0)e-*("°-'')* (143) 
^t(^) = giHt^t(o)e-^-f^* = ^t(o)e^(^o-'^)*. (144) 

We will study imaginary time quantum mechanics for which the time evolution operator is 

U{t) = e-"^ (145) 

and in which 

*(r) = *(0)e-(^°-'^)^ (146) 
^t(^) = ^t(o)e(^o-^)^. (147) 

Note that despite the notation, ^{t) and ^'^(t) are are not adjoints except at r = owing 
to the fact that U (r) is not unitary. 

Readers not familiar with imaginary time quantum mechanics merely have to observe 
how the functional formalism reproduces the results of the operator formalism. They may 
also wish to learn about imaginary time quantum mechanics using this simple example. 
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Next consider the time-ordering symbol T whose action on a pair of fermionic Heisenberg 
operators is: 

r(*(T)*^(0)) = *(t)*^(0) t > (148) 
= -^^(O)^(t) t < 0. (149) 

Note that 

Ar= hm -r(*(r)*'^(0)). (150) 

r— »0~ 

In field theory and many-body physics one is interested in the Green's function: 

G{t) =< T(*(t)*^(0)) > (151) 

where <> denotes the average with respect to Z. For our problem we find 

^ < Oir(vP(r)vpt(0))|0 > + < l|r(vP(r)vF(0))|l > c-'"'"-^) 
9{r)e~^^°~>'^'^ - 6'(-r)e"(^0"'')('"+^) 

In the zero temperature limit this reduces to 



(153) 



G{t) = ^(r)e-(^°-'^)^ iKflo (154) 
= -e{-T)e-^^°-^^^ H>no- (155) 



Let us define the pair of transforms: 



/oo 
G{T)e'''^dT, (156) 
-oo 

G(r) = /_^G(u;)e— ^. (157) 



We find that 



Gioo) = ^— (158) 

\lo — IJ, — lU 
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independent of which of Qq or /i is greater. 

Let us calculate < N > using Eqn. ( P^ ) and the above Green's function: 

<A^> = -^(0") (159) 



(160) 

oo 2tt lU — fi — ilo 

= e{f,-Qo) (161) 

where in the last step we have argued that unless the 6'-function is satisfied, the contour must 
be closed in the upper half plane (as dictated by the exponential) will be free of singularities. 
Consider finally a toy "Hubbard model" with two fermions and a repulsive interaction 

U: 

Ho = l]o(^l^l + *2^2) + t^^l*1^2*2 (162) 

= noN+^NiN-1) (163) 

where 

N = Ni + N2 = + ^5^2- (164) 

(Readers wishing to fill in the missing steps should note that they have to use Nf = Ni and 
likewise for N2. In any event they should check the correctness of the final result for various 
choices of A^i and A^2-) Each fermion has the usual anticommutator of unity with its adjoint 
and anticommutes with everything else including all member of the other set. 
By summing over the eigenstates of A^i and N2, 

Z = l + 2e-^(^°-'^) + e-'3(2(no-M)+c/)) ^^g^) 

where the factor of 2 in the middle term reflects the degeneracy of the one fermion states. 
From the above, we obtain by differentiation 

2(1 + g/3(M-^^o-C/)) 

< iV > = lim — sr-rn ^7 — k — ttt (166) 
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= // < Qo 



(167) 



= 1 Qo < l^<^o + U 



(168) 



(169) 



Table III summarizes results from this subsection. In the next subsection they will be 
rederived using path integrals. 

B Fermion Coherent States 

In this section we will be using Grassmann numbers. Here are the rules for manipulating 
them: 

• All Grassmann numbers anticommute with each other and with all fermionic operators. 

• As a result of the above, the square of any Grassmann number is zero and the product 
of an even number of Grassmann numbers will commute with anything. Likewise any 
Grassmann number will commute with an even number of fermion operators such as 
N — When a Grassmann number is taken through a ket or bra containing an 
even (odd) number of fermions it will not (will) change sign. 

• Do not associate a numerical value to Grassmann numbers. There are no big or small 
Grassmann numbers. All you will need are the above definitions. 

Consider the state 



where is a Grassmann number. This state, called a fermion coherent state is an eigenstate 
of with eigenvalue ■0: 



1^ |0 > --0|1 > . 



(170) 



>= > 



(171) 
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as is readily verified: 

^|?/^> = *|0>-^^/'|l> (172) 

= + V'^|l> (173) 

= ?/^|0 > (174) 

= ^(|0>-^|1>) (175) 

where we have appealed to the fact that ip anticommutes with \l/ and that ip'^ = 0. If we 
act on both sides of Eqn. (|171|) with the left vanishes due to \E'^ = and the right due to 

ij^ = 0. 

It may be similarly verified that 
where 

<^| =< 0|- < 1|V^=< 0| 1|. (177) 

Please note two points. First, the coherent state vectors are not the usual from a complex 
vector space since they are linear combinations with Grassmann coefficients. Second, ip is 
not in any sense the complex conjugate of and < ip\ is not the adjoint of IV' >• You should 
therefore be prepared to see a change of Grassmann variables in which ip and ip undergo 
totally unrelated transformations. 

The inner product of two coherent states is 

<^\iP> = (< 0|- < 1|V^)(|0 > -^|1 >) (178) 
= < 0|0 > + < ll^^Z-ll > (179) 
= l + i^^p (180) 



181) 
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Any function of a Grassmann variable can be expanded as follows: 

F(V') = Fo + FiV' (182) 

there being no higher powers possible. 

We will now define integrals over Grassmann numbers. These have no geometric signifi- 
cance and are formally defined. We just have to know how to integrate 1 and since that 
takes care of all possible functions. Here is the list of integrals: 



ipdip = 1 (183) 
J Idil: = 0. (184) 

That's it! As you can see, selling tables of Grassmann integrals is no way to make a living. 
There are no hmits on these integrals. Integration is assumed to be a hnear operation. The 
differential ^-0 is also a Grassmann number. Thus / dil^il) = — 1. The integrals for -0 or any 
other Grassmann variable are identical. A result we will use often is this: 



j ipipdipdip = 1. (185) 



Note that if the differentials or variables come in any other order there can be a change of 
sign. For example we will also invoke the result 



^^dil^dil^ = -1. (186) 

Let us now consider some gaussian integrals. You are urged to show the following: 

j e-'^'^'^d^dilj = a (187) 
J e-^^^[##] = detM (188) 

where in the second formula M is a 2-by-2 matrix, ■0 is a column vector with entries ipi 
and ip2, ip ^ column vector with entries ipi and ip2 and [dipdip] = d'4'idipidip2dil'2-T^i^is result 

51 



is true for matrices of any size. To prove these simply expand the exponential and do the 
integrals. 

Consider next the "averages" over the gaussian measure: 

The proof is straightforward and left as an exercise. 

Consider now two sets of Grassmann variables (labeled 1 and 2). The reader should 
verify that 

= ^=<a>. (191) 

We now have a Wick's theorem for fermions: 

< w-'ih-tbi.'iJji > = — = = = = 192 



(193) 

< il X jk > - <l.k >< jl > . (194) 



The reader not familiar with such objects is urged strongly to prove this simple case of 
Wick's theorem for fermions. Note the strong similarities to the bosonic case. Once again 
we find that the answer is zero unless each Grassmann is accompanied by its partner. The 
answer is once again a sum over all possible parings. The only difference comes from the 
minus signs which are determined as follows. We first move each variable till it is next to its 
partner. In the example above, if j = k and i = I, the middle two Grassmanns are already 
next to each other and the ones at the ends can be brought together without any minus 
signs since they are separated by a pair of Grassmanns. This is why the first term is positive 
in Eqn.(|T9^). On the other hand if i = k and j = I, we must move j through k to meet 
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its mate and this produces a minus sign. When more than four variables are averaged, an 
obvious generahzation holds: pair the fields in all possible ways and put in a minus sign for 
every time a variable crosses another. Although we did not see it here, the following thing 
can happen and does happen in a calculation that comes later in this section: the variable 
and its partner are next to each other, but in the wrong order, with ip to the left oi ip. In 
this case an extra minus sign is needed to rearrange these. 

Finally note that Jacobians behave counterintuitively for Grassmann variables. Consider 

j ai/jdil) = a. (195) 

In terms of 

aV- = X, (196) 

JxJ{i^/x)dx = a. (197) 

Assuming J is a constant we pull it out of the integral (with no minus signs since it involves 
an even number of Grassmann variables) , use the fact that the integral of x is unity to 
obtain 

J{^/X) = a (198) 

while one might have expected the inverse. 

As an application of this resut, the reader may wish to rederive Eqn. (|187l) by making the 
change of variables from ifj to x = flV"- Note that there is no need to transform tp at the same 
time. 

We need two more results before we can write down the path integral. The first is the 
resolution of the identity: 

J = j \ip>< '^le-^^di^dijj (199) 
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In the following proof of this result we will use all the previously described properties and 
drop terms that are not going to survive integration. (Recall that only ■^V' = —■0-0 has a 
non-zero integral.) 

J >< ^|e-^^#d^ = J\i/J>< ^1(1 - i^i/j)cfipdi/j (200) 

= y(|0>-V'|l>)(<0|-<l|V^)(l-V^V)## (201) 
= y(|0 >< 0| +V'|1 >< 1|V^)(1 - V^V')#ciV^ (202) 
-i^i/j)(hpdi/j + |1 >< 1| j ifip(kpdip (203) 

(204) 



= 0X0 



The final result we need is that for any bosonic operator (an operator made of an even 
number of Fermi operators) 

TrVL = j < > e'^'^'d^dip. (205) 

The proof is very much like the one just given and is left to the reader. 
C The fermionic path integral 

Consider the partition function for a single oscillator: 

Z = 2^^g-/3(no-M)*t* (206) 
= y < -^le-^^^^o-^)*'*!^ > e-^'^di^di/j. (207) 

You cannot simply replace ^'^ and ^' by —ip and ip respectively in the exponential. This is 
because when we expand out the exponential not all the ^'s will be acting to the right on their 
eigenstates and neither will all ^'^'s be acting to the left on their eigenstates. (Remember that 
we are now dealing with operators not Grassmann numbers. The exponential will have an 
infinite number of terms in its expansion.) We need to convert the exponential to its normal 
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ordered form in which all the creation operators stand to the left and all the destruction 
operators to the right. Luckily we can write down the answer by inspection: 

g-/3(no-/.)*t$ ^ ^ ^ (g-/3(no-/x) _ ;^^^t^ (208) 

whose correctness we can verify by considering the two possible values of ^f^^. (Alternatively 
you can expand the exponential and use the fact that N'^ = N for any nonzero k.) Now we 
may write 

Z ^ J + (e-/^("o-^) - l)^t^|^ > e-^'^## (209) 

= y < > (1 + {e-f^(^°-i''> - l)(-V^V))e"^^## (210) 

= J{1- (e-^^^^o-^) - l)V^V^)e-2^'^#d^ (211) 

= l + e-^("°-'') (212) 

as expected. While this is the right answer, this is not the path integral approach. As for 
the latter the procedure is this. Consider 

Z = Tre-^" (213) 

where if is a normal ordered operator H{^\ ^). We write the exponential as follows: 

e-^^ = lim(e-^^)^ (214) 

AT— »oo 

= (1 - c H) ... (1 - cH) e = (3/N (215) 



N times 

and introduce the resolution of the identity N — 1 times: 

Z ^ y < -^^1(1 - sH)\il:^_, > e-^~-i^^-i < V^^,_i|(l - eH)\ilJr,-2 > e"^^-^^^"^ < i:N-2\ 

_ _ N-l 

. . . > e-^^V-a eH)\jj;i > e'^i^^ J] (216) 
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Now we may legitimately make the replacement 

< - eH{<il\ >=< - eH{^,^„^i)\^i >= e^»+i'/'>e-^^(^'+-'^») (217) 

where in the last step we are anticipating the limit of infinitesimal e. Let us now define an 
additional pair of variables (not to be integrated over) 

i^N = -^1 (218) 
V^AT = -^1. (219) 



The first of these equations allows us to replace the left-most bra in Eqn.( ^16'[ ), < —ipi\, by 
< ipi^^l- The reason for introducing ip^ will follow soon. 

Putting together all the factors (including the overlap of coherent states) we end up with 

/.N-l _ _ 

/ 1 I [ e'^'+i'^'e-"^('^'+i''^')e-'^^'^>#i#i (220) 



i=l 



I n 



'd^idijji (221) 



i=l 

/o''^(-)(-i^-f^o+M)V'Wrfrj^^^^| ^222) 



where the last step need some explanation. With all the factors of e in place we do seem to get 
the continuum expression in the last formula. However the notion of replacing differences by 
derivatives is purely symbolic for Grassmann variables. There is no sense in which — ip^ 
is small, in fact the objects have no numerical values. What this really means here is the 
following. In a while we will trade ipi'^) fo^' "^i^) related by Fourier transformation. At 
that stage we will replace — ^ hy iuj while the exact answer is e*"^ — 1. If we do not make 
this replacement, the Grassmann integral, when evaluated in terms of ordinary numbers, 
will give exact results for anything one wants to calculate, say the free energy. With this 
approximation, only quantities insensitive to high frequencies will be given correctly. The free 
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energy will come out wrong but the correlation functions will be correctly reproduced. (This 
is because the latter are given by derivatives of the free energy and these derivatives make 
the integrals sufficiently insensitive to high frequencies.) Notice also that we are replacing 
-f^(V'i+i5 V'i) = H{ip{T + e), Tpir)) by H{iP{t), iP{t)) in the same spirit. 
Now turn to the Fourier expansions alluded to above. Let us write 

= E-5-V^(^) (223) 

n P 

p-iuJnT 

= E— (224) 

n P 

where the allowed frequencies are chosen to satisfy the antisymmetric boundary conditions 
in Eqn.(|18|-|l|). Thus 



= (225) 

where n is an integer. Note that we have chosen the Fourier expansions as if ip and ip were 
complex conjugates, which they are not. This choice however makes the calculations easy. 
The inverse transformations are 

^iu) = [%(T)e''^"^dT (226) 
Jo 

^(w) = f^ir^e-'^^-^dT, (227) 

JO 

where we use the orthogonality property 

/ e'^-^e-'^-^dT=— - = l35mn. (228) 

Jo l{(jJn-UJm) 

Performing the Fourier transforms in the action and changing the functional integration 
variables to ip{uj) and ip{uj) (the Jacobian is unity) and going to the limit /5 — >■ oo, which 
converts sums over discrete frequencies to integrals over a continuous u;, we end up with 

J-°l I^V^HC— f^o+/.)^(-) (229) 
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Although (3 has disappeared from the picture it will appear as 27r(5(0), which we know stands 
for the total time. (Recall Fermi's Golden Rule calculations.) An example will follow shortly. 

Let us first note that just as in the case of the scalar gaussian model, the correlation 
function is related to the integral over just a single pair of variables (Eqn.([TH^)) and is given 
by: 

/V^(o.i)^(c..)en^^(-)^-"°+^)-^(-)[#(c.)^^(a;)] 
^ /#(.).^(.)eni^^H(^^-o^.)^H ^'''^ 



iuji — f2o + /i 

In particular 



(231) 



< ^(-)^(-) >= . = -A (232) 

lUl — \ Iq + 11 lUJ — \ Lq + li 



Let us now calculate the mean occupation number < N >: 



1 dZ 

= 7^/ —<H^m^)> (234) 



doj e 



ILOO ' 



(235) 

-00 277 ioJ — + fJL 

= 9ifi-no) (236) 

as in the operator approach. Notice that we had to introduce the factor e*^°^ into the u 
integral. We understand this as follows. If we had done the calculation using time r instead 
of frequency u, we would have calculated the average of \E'^\E'. This would automatically have 
turned into iIj{t + e)tp{T) when introduced into the path integral since the coherent state 
bra to the left of the operator would have come from the next time slice compared to the 
ket at the right. (Remember how if(\E'"'^, \E') turned into H[ip{i + l)'tp{i)).) Notice that the 
integral over uo was not convergent, varying as duj/uo. It was therefore sensitive to the high 
frequencies and we had to intervene with the factor e^^°"^ . Later we will deal with integrals 
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that have two or more powers of lo in the denominator and are hence convergent. We will 
not introduce this factor in those cases. 

Our final calculation will be the determination oi < N > for the toy Hubbard model 
using path integrals. The partition function is 

Z = y[#i(iV'i#2t^V'2]e^°e^/^i'^2'^i^2 ^;^g^g (237) 

/oo 2 
— ^^,(a;)(^a;-l]o + /^)^iM, (238) 
-oo ZTT ^ 

/"OO ^ duj- 

^i^2^i^2 = / n l7^^i('^4)^/'2('^3)V'i('^2)V^2(<:^i)27r5(a;4 + a;3 - a;2 - a;i).(239) 

J—oo ZTT 

Note that the hamiltonian already was in normal ordered form: each creation operator 
was to the left of its destruction operator. This allowed us to replace the operators by 
their eigenvalues when the coherent states were introduced. For notational uniformity we 
have further arranged to have all creation operators to the left of all destruction operators. 
This merely introduces an extra minus sign here since operators corresponding to different 
oscillators anticommute. 

Prior to calculating < > let us calculate < iVi >. This is given by 

— < M^)M^) > e'^"^ (240) 

-oo ZTT 

where < ■0^(a;)'0i(a;) > stands for the correlation function with the full action and not just 
Sq. We may however express it in terms of averages over the gaussian measure Sq using 
the same trick we used for bosons: we multiply and divide the exact Z hy Zq, the partition 
function with U — to obtain 

^/ N,/ N < V^i(o;)^i(o;)e^/^i'^2V'i^2 > 

In principle the reader has all the information needed to evaluate this expression to any 
order in perturbation theory. The exponential is to be expanded (in the numerator and 
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denominator) to the desired order and all averages done using Wick's theorem. We shall 
carry out the calculation to order U to show the details. The result to all orders will then 
simply be stated and the sceptics are encouraged to check it to higher orders. The integration 
measure and delta functions will be occasionally suppressed. 
To first order we find 



_ < ?/'i(a;)?/'i(0j) >0 +f//(V'l('^)^l(^)^l(t^4)^2('^3)^l(^2)V'2(^l),„ ^ ^ 

< ^iH^i(cu) >= = ^ ^(242) 





Both in the numerator and denominator we have just one pair of fields with label 2; these 
must obviously be paired, though it take a sign change to bring them next to each other. 
The four fields with label 1 can be paired in two ways. One of the ways, in which the 
external fields iPi{uj)i/ji{lj) get paired is precisely cancelled by the denominator expanded 
out to order U. This corresponds to the cancellation of disconnected diagrams. Reader new 
to this concept are urged to verify this. What remains may be written as follows: 

< ^1 u; ViM > = FT^ — + T- o / 7^ (243 

too — \lo + fi [lio — \Iq + jj,)"^ J-oo ivr tuji — \lo + jj, 

^ ^ (244) 



itu - ^lo - fx - U < N2 > 

In going to the last step we have taken two terms of the power series in U assumed that they 
represent a geometric series and summed the series. This result is undoubtedly correct to 
the order we are working in. We have also replaced the integral over ui with < N2 > which 
is also good to this order. It turns out that both these approximations are in fact exactly 
what we would get if we went to all orders, as will be explained shortly. Let us accept this 
for the present and see what follows. Using the above correlation function into the formula 
for < A^i >, Eqn.(plO|), we obtain 

< Ni>=e[iJ-no-U < N2>]. (245) 
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Obviously we can similarly derive another equation 

< N2>=e[i^-no-U < Ni>]. (246) 

Let us explore these equations for various cases. First if /x < Qq, neither 6 -function can 
be satisfied since < Ni >> 0. Thus we get < >= for this before. Likewise if 

fi > Qq + U, both ^-functions will be satisfied satisfied since < Ni >< 1. This gives us 
< >= 2 as before. Finally consider Qq < fi < Qq + U. Since each < Ni > equals a 
^-function, it can equal only or 1. It is readily seen from these equations that the only two 
consistent choices are < A^i >= 1,< N2 >=0 and vice versa, once again in agreement with 
the operator solution. 

Now for the higher terms in the expansion. These are best seen in diagrammatic terms. 
Consider Figure 4a. The dark line stands for the full propagator and the thin line for 
the one computed in the gaussian measure. To order U we kept the two diagrams shown 
and these correspond to the expressions in Eqn.(2^). The disconnected diagram that got 



cancelled by the denominator is shown in Figure 4b. If we go to higher orders we can run into 
either iterates of the one loop or embellishments of it. The embellishments convert the loop 
integral over the free propagator of species 2 to the integral over the full propagator which 
then reduces to < N2 >. The iterations produce the remaining terms in the geometric series 
that was presumed in going from Eqn. (|243|) to the next one. This leaves us with diagrams 
such as the "sunrise" diagram in Figure 4d. (Once again the nomenclature is from field 
theory.) These diagrams and all the rest vanish in this problem because the corresponding 
frequency integrals are convergent and have all the poles on the same half-plane, allowing us 
to close the contour the other way. 

The generalization of Grassmann integrals to many-body problems is straightforward. 
The labels 1 and 2 from the toy Hubbard model can run over, say the modes in the Brillouin 
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zone. The action is once again obtained by replacing the normal ordered hamiltonian by the 
corresponding function of Grassmannian coherent state labels. As for the coupling functions, 
just as the coupling in the bosonic it(4321) was symmetric under the exchange of the first or 
last two labels among themselves, the fermionic couplings will be antisymmetric under such 
an exchange due to the anticommuting nature of the Grassmann variables. 

Finally a matter of notation. We will be switching from upper case letters ^ and for 
fermion operators to lower case. It will be clear from the context whether we are referring 
to the operators or the Grassmann variables. 

Table IV summarizes the results from the discussion on fermionic path integrals. 

IV MOTIVATION AND WARMUP: RG IN = 1 

We are now ready to turn to the main topic: application of the RG to interacting nonrela- 
tivistic fermions. The best way to explain the method is to deal with specific problems to 
which it applies. We begin with the problem of charge density wave (CDW) formation in a 
system of spinless fermions at half-filling. It has relatively simple kinematics and illustrates 
the RG approach very nicely. In fact the methods explained here were originally developed 
(Shankar 1991) to deal with this problem in two dimensions. Let us begin with a discussion 
of the various terms used above in describing the model. 

We begin with the justification for the study of spinless fermions, which are admittedly 
a theorist's construction. As we progress with the paper, it will become apparent that the 
RG is primarily concerned with the the symmetry properties of the Fermi surface of the 
noninteracting fermions. For example the superconducting instability for arbitrarily small 

— * 

attraction is due to the invariance of the Fermi surface under time reversal: if K lies on the 
surface, so does —K. Likewise the charge density wave (CDW) instability on the square 
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lattice, towards a ground state in which the charge density is nonuniform and oscillates 
between the sublattices, is due to nesting: if K lies on the Fermi surface, so does K + Q]\f, 
where Qat is a fixed nesting vector. Now these very same properties of the Fermi surface 
will also destabilize a system of real electrons. The actual nature of the instability can 
of course be different in the two cases. For example the nested Fermi surface in d = 2 
will cause electrons to go to an antiferromagnetic state in which the magnetization (rather 
than charge density) oscillates in magnitude between the two sublattices. Likewise the time 
reversal invariant Fermi surface will lead to superconductivity but (Cooper) pairing can 
take place for any angular momentum, while for spinless fermions only odd orbital angular 
momentum states are allowed due to antisymmetry requirements. To summarize, spin really 
is an inessential complication if we are simply tying to understand how the RG works. When 
comparing theory to experiment, spin will of course have to be included, but this really will 
be straightforward. 

Consider now the requirement of half-filling, which means that the system has half the 
maximum number allowed by the exclusion principle. For spinless fermions this means one 
particle for every other site, while for electrons it means one particle per site. Despite this, 
the two problems will have the Fermi surface in the noninteracting case. In both cases, the 
Fermi surface will enclose half the Brillouin zone and have the same shape, decided by the 
lattice parameters. However each filled momentum state will carry two electrons (of opposite 
spin) but just one spinless fermion. In other words, the condition of half-filling implies a 

different particle density in the two cases, but the same Fermi surface . ^ 

^An aside for readers with a different background, say particle pliysics, wlio are troubled by the following: 
why bother with effects at such a special filling ? Is there any chance that a generic system will have filling 
factor of exactly half, as compared to say .51? Yes! Consider a square lattice with an atom at each site. Since 
each atom contributes an integer number of electrons to the conduction band, the filling factor (the ratio of 
the number of electrons to the maximum allowed number of two per site) is bound to be a half-integer or 
integer. For a more complicated unit cell, or intercalated compounds, one can have other simple fractions 
like 1/4. 
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We start here with the one dimensional version of the spinless fermion problem for two 
reasons. First, the one dimensional problem is very interesting in itself and shows the power 
of the RG approach. Secondly, as mentioned earher, due to the fact the Fermi surface in 
d — 1 consists of just two points, the problem resembles quantum field theory with a few 
coupling constants (rather than a theory with coupling functions) and affords a painless 
introduction to the use of the RG for fermions. We will then be better prepared for the 
d — 2 version in Section X. 

A The d = 1 model: definition and mean-field analysis 

Let us consider the following specific hamiltonian for a spinless fermion system on a d = 1 
lattice labeled by an integer n: 

H = Hq + Hi (247) 

= -IJI^Kj + I)^(J) + h.c. + UoY.ii^\mj) + 1)V'(J + 1) -^)(248) 
j j 

where the fields obey 

{7P^j),^{m)}^5^j (249) 

with all other anticommutators vanishing. 

The first term represents hopping. The hopping amplitude has been normalized to 1/2. 
The second term represents nearest neighbor repulsion of strength Uo- The role of the 1/2- 's 
subtracted from the charge densities rij (= ip^jipj) and n^+i is this. When we open up the 
brackets, it is readily seen that they represent a chemical potential 

= Uq. (250) 

This happens to be exactly the value need to maintain half-filling in the presence of the 
repulsion Uq. To see this, make the change ip ^ ip"^ at all sites. This exchanges the site 
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occupation number n — ip^ip with 1 — n or changes the sign of n — 1/2. Thus both brackets 
in the interaction term change sign under this, and their product is unaffected. As for the 
hopping term, it changes sign under ip ip^ . This can be compensated by changing the 
sign oi ip on just one sublattice (which preserves the anticommutation rules and does 
not affect the other term) . Thus H is invariant under exchanging particles with holes. This 
means the ground state (if it is unique) will satisfy < n >=< 1 — n > which in turn means 
< n >— 1/2. (If there is a degeneracy of the ground state, the result still holds, but takes a 
little more work to establish.) 

Let us understand this model in the extreme limits Uq — and Uq — oo. 

As for the first case let us introduce momentum states via 

^ ^^(i^)e^^^' (251) 

and the inverse transform 

m)=J2e-'^^m- (252) 

3 

Using J2j — 27r(5(0), we can verify that 

{i>{K),i)'i{K')} = 27r5{K - K'). (253) 

In terms of these operators 

—i;\K)ij{K)E{K) (254) 

-TT iTl 

E{K) = -cosK. (255) 

The Fermi sea is obtained by filling all negative energy states, i.e., those with \K\ < Kp — 
7r/2, which corresponds to half- filling. The Fermi surface consists of just two points = 
±7r/2. It is clear that the ground state is a perfect conductor since we can move a particle 
just below the Fermi surface to just above it at arbitrarily small energy cost. 
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Consider now the situation in the other extreme Uq = oo. We now ignore the hopping 
term and focus on just the interaction. It is evident that the lowest energy states are those 
in which no particle has a neighbor: thus either the A-sublattice consisting of even sites 
is occupied, or the B- sublattice, made up of the odd sites, is occupied. This makes the 
product {rij — l/2)(nj+i — 1/2) negative on every bond. These two states, which break 
the translational symmetry of the lattice are the CDW states. The order parameter, which 
measures the difference between the mean occupation of the odd and even sites is maximal 
(unity). In the CDW state, the system is an insulator. Any excitation of the ground state 
requires us to move the charge and this will cost an energy of order [/q. (This is clearest 
as Uq — >■ oo.) One expects that even for large but finite Uq, the symmetry would still be 
broken, but with a smaller order parameter. 

Here is the question we want to answer: will the system develop the CDW order and gap 
for arbitrarily small repulsion, or will it remain a conductor up to some finite Uq ? 

We will use the RG to answer it. But first let us see what a very standard tool, namely 
mean field theory, can tell us. In this approach one assumes a CDW order parameter in 
the ground state and asks if the assumption is self-consistent. The self-consistency check is 
approximate, as will be explained. Mean-field theory predicts that CDW will set in for any 
repulsion however small. Here is a short description of the calculation. 

Let us begin with Eqn. (|248|) and make the ansatz 

<n,>=\ + \{-iyA (256) 

where A is the CDW order parameter. We will now see if the ground state energy of the 
system is lowered by a nonzero value of A. To this end, we will find the ground state energy 
as a function of A and minimize it and see if the minimum occurs at a nonzero A. However 
this last step will be done approximately since this is an interacting many body system. The 
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approximation is the following. We start with the interaction 



^ = E ^'PKj + mj) + h.c + Uo - - ^) 

j 3 



(257) 



and make the substitution 



(258) 
(259) 



where the normal ordered operator : Uj : has no expectation value in the true ground state 
and represents the fluctuations in number density. Upon making these substitutions and 
some rearrangements we find 



j 

+Uo[\e^' - AE(-i)'%] + UoE-- % ■■■■ rij+i ■■ 

j j j 



(260) 
(261) 



In the mean field approximation we ignore the last term. The rest of the hamiltonian is 
quadratic and solved by Fourier transformation . Due to the factor (—I)-' which multiplies 
rij, states with momentum K and K' — K + tt will mix. The hamiltonian becomes 

A2 



Jo iTT 



E{K) -UqA 
-UoA E{K') 



+ Uo27r5(0)- 



(262) 



Notice that we have halved the range of K integration, but doubled the number of variables 
at each K. The two-by-two matrix, which is traceless due to the relation 



E{K') = - cos{K + tt) = -E{K), 



(263) 



is readily diagonalized. The one-particle energy levels come in equal and opposite pairs and 
we fill the negative energy states to obtain the following ground state energy per unit volume: 

^0 A^Uo r dK 



27r5(0) 



(IK I 



(264) 
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Uj '-^ , ' (266) 



where the integral comes from the filled sea. Minimizing with respect to A we obtain the 
relation 

A=r^, ""^ . (265) 

Assuming A 7^ 0, we cancel it on both sides. It is clear that C/q < is not acceptable since 
the two sides of the equation would then have opposite signs. For positive Uq, a nontrivial 
solution requires that 

dK 

which is called the gap equation. On the left hand is the number 1 and on the right hand 
side, something of order Uq. It appears that we will get a solution only above some minimum 
Uq. This is wrong. The integrand becomes very large at the Fermi points = Kp = 7r/2, 
where E{K) vanishes. Writing 

E(K) = k (267) 
k = \K\- Kf (268) 

we approximate the gap equation as follows: 

l = uj''-t^ ' (269) 

J-A TT ^}.2 + A2[/2 TT AUq 

where A, the upper cut-off on is not very important. What is important is that due to 
the logarithmic behavior of the integral near the origin in k, i.e., near the Fermi surface, 
there will always be a solution to the gap equation given by 

A = ^e-'^/'^". (270) 

The logarithmic divergence is also reflected in the divergent susceptibility of the noninteract- 
ing system to a probe (or perturbation) at momentum tt. (At second order in perturbation 
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theory, the perturbation will link the ground state to states of arbitrarily low energy in which 
a particle just below the right (left) Fermi point is pushed to just above the left (right) Fermi 
point. The small energy denominators, summed over such states will produce the logarithm.) 

Mean-field theory also predicts that the same thing will happen in o? = 2. In this case the 
nesting condition (readers unfamiliar with this concept should consult Figure 17, Section X) 
ensures that the perturbation at (tt, tt) will excite particles just below the Fermi surface to 
states just above it on the "other side" no matter where the starting point is on the Fermi 
surface. Now for any Fermi surface, if we take some perturbation of any momentum (which 
is not too large compared to the size of the surface ), there will always be some points just 
below sea level, that will get knocked to points just above sea level. But these points would 
have to come from some special angular region of the Fermi surface which are connected by 
this momentum (whose angular dimensions will decrease with the energy denominators) and 
the integral over the small energy denominators will converge. If we use a coordinate 9 on 
the surface and a coordinate e normal to it, the integrals will be off the form 



The cut-off A focuses on the small energy denominator region. In the absence of nesting, for 
a perturbation of some given momentum Q, the range of 6 will shrink with e and the integral 
will converge. On the other hand, at the nesting momentum, the angular integral will over the 
entire Fermi surface (since every point on or near the Fermi surface get scattered to another, 
as shown in Figure 17, Section X) and the integral will have a logarithmic divergence. 

Returning to d = 1, mean- field theory also predicts that the system will a have non-zero 
superconducting order parameter for the smallest attractive coupling. In the corresponding 
calculation the instability will stem from the time-reversal symmetry of the problem: E(K) — 




(271) 



E{-K). 
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Unfortunately both these predictions are wrong. The error comes from the neglected 
quartic operator : rij :: %+i :. We know all this because the present spinless hamiltonian 
Eqn. ( |248D can be solved exactly (Yang and Yang 1976). Q The exact solution tells us that 
the system remains gapless for Uq of either sign until it exceeds a minimum value of order 
unity. We will now develop the RG approach to this problem and obtain results in harmony 
with this exact result. 

B The RG approach for d =1 spinless fermions 

Our goal is to explore the stability of the noninteracting spinless fermions to weak interac- 
tions. We are not interested in the fate of just the model in Eq.( p48| ) but in a whole family 
of models of which this will be special case. Our strategy, stated earlier on, is the following. 
First we argue that at weak coupling, only modes near -^Kp will be activated. Thus we will 
linearize the dispersion relation E{K) = —cosK near these points and work with a cut-off 
A: 

/A f]U 
-4[kmk)k (272) 

where 

k = \K\- Kf (273) 
i = L,R {left or right). (274) 

Notice that Hq is an integral over fermionic oscillators which we studied in Section III. The 
frequency Qq of the oscillator at momentum k is simply k. 

Next we will write down a T = partition function for the noninteracting fermions. 

This will be a Grassmann integral only over the degrees of freedom within a cut-off A of 

^'^ The reader consulting this reference wiU find that these authors solve a problem of quantum spins on 
a line. This XXZ chain is related to the spinless fermions by a Jordan Wigner transformation . 
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the Fermi surface. We will then find an RG transformation that lowers the cut-off but 
leaves the free-field action, Sq, invariant. With the RG well defined, we will look at the 
generic perturbations of this fixed point and classify them as usual. If no relevant operators 
show up, we will still have a scale-invariant gapless system. If, on the other hand, there are 
generic relevant perturbations, we will have to see to which new fixed point the system fiows. 
(The new one could also be gapless.) The stability analysis can be done perturbatively. In 
particular, if a relevant perturbation takes us away from the original fixed point, nothing at 
higher orders can ever bring us back to this fixed point. The fate of the nearest neighbor 
model will then be decided by asking if it had a relevant component in its interaction. 
Let us then begin with the partition function for our system of fermions: 

^0 = / n n d'ilJi{ujk)d'^i{uk)e^° (275) 

i=L,R\k\<A 

/A fjh POO rji J 
IT ^U^muj - k)U^k) (276) 

This is just a product of functional integrals for the Fermi oscillators at each momentum 
with ^o{k) — k. 

The first step in the RG transformation is to integrate out all ip{ku!) and ip{ku!) with 

A/s < |A;| < A (277) 

and all lu. Thus our phase space has the shape of a rectangle, infinite in the uj direction, 
finite in the k direction. This shape will be preserved under the RG transformation. Since 
there is no real relativistic invariance here, we will make no attempt to treat uu and k on 
an equal footing. Allowing cu to take all values allows us to extract an effective hamiltonian 
operator at any stage in the RG since locality in time is assured. 

Since the integral is gaussian, the result of integrating out fast modes is just a numerical 
prefactor which we throw out. The surviving modes now have their momenta going from 
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—A/s to A/s. To make this action a fixed point we define rescaled variables: 

k' = sk 
ijj' — sou 

ip[{k'uj') = s-^'^iPi{kuj) (278) 

Ignoring a constant that comes from rewriting the measure in terms of the new fields, we 
see that 5*0 is invariant under the mode elimination and rescaling operations. 

We can now consider the effect of perturbations on this fixed point. Rather than turn on 
the perturbation corresponding to the nearest neighbor interaction we will perform a more 
general analysis. The result for the particular case will be subsumed by this analysis. 

C Quadratic perturbations 

First consider perturbations which are quadratic in the fields. These must necessarily be of 
the form 

/A rjl^ fOO fj/,^ 
-/ —^^{kuj)^,{ujk)^,{ujk) (279) 
z=j., u 27r J-oo 27r 

assuming symmetry between left and right fermi points. 

Since this action separates into slow and fast pieces, the effect of mode elimination is 
simply to reduce A to A/s in the integral above. Rescaling moments and fields, we find 

H'iou', k', i) = i). (280) 

We get this factor as a result of combining a factor from rewriting the old momenta 
and frequencies in terms of the new and a factor which comes from rewriting the old fields 
in terms of the new. 

Let us expand in a Taylor series 

li{k, uj) = iiQQ + iiiok + iiQiiuj H h i^nmk'^iiuj)"' H (281) 
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The constant piece is a relevant perturbation: 

/^oo — ' snoo- (282) 

This relevant flow reflects the readjustment of the Fermi sea to a change in chemical potential. 
The correct way to deal with this term is to include it in the free field action by filling 
the Fermi sea to a point that takes /ioo into account. As for the next two terms, they 
clearly modify terms that are already present in the action, can be absorbed into them and 
correspond to marginal interactions. When we consider quartic interactions, it will be seen 
that mode elimination will produce terms of the above form even if they were not there to 
begin with just as 0^ theory. The way to deal with them will be discussed in due course. As 
for higher order terms in Eqn.( p8TD , they are irrelevant under the RG mentioned above. This 
is however a statement that is correct at the free-field fixed point. We shall have occasion 
to discuss a term that is irrelevant at weak coupling but gets promoted to relevance as the 
interaction strength grows. 

D Quartic perturbations: the RG at Tree Level 

We now turn on the quartic interaction whose most general form is 

SS, = ^ |^^V^(4)V^(3)^(2)^(1)«(4,3,2,1) (283) 

where 

V^(0 = i^iKi,iu,) etc., (284) 

^ dK- duj- — 

' 27r5(is:i + fs:2-^3-^4)27r5(tUi + tU2-cj3-tU4) (285) 



Kuj 



n 



i=l 



'-IT 27r J-oo 27r 

and 5 enforces momentum conservation mod 27r, as is appropriate to any lattice problem. A 
process where lattice momentum is violated in multiples of 27r is called an umklapp process. 
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The delta function containing frequencies enforces time translation invariance.Tlie coupling 
function u is antisymmetric under the exchange of its first or last two arguments among 
themselves since that is true of the Grassmann fields that it multiplies. Thus the coupling u 
has all the symmetries of the full vertex function T with four external lines. 

To get a feeling for all these ideas let us consider the nearest-neighbor repulsion from 
Eqn. (P^5D and ask what u it generates in the action. Let us first begin with the operator 

Hi = UoYl ^^-^i^^'+i^i+i = -Uo E (286) 
j j 

and make a Fourier transformation to get 
dK 



Hi = -Uo 



n 



i=l 



2n5iKi + K2-Ks-Ki)i;^{Ki)i;^{Ks)i'{K2)ijiKi)e'^'''-'''l (287) 



27r 

We now antisymmetrize e^^^'^~^^'' with respect to 1 ^ 2 and 3^4 since the rest of the 
integrand is antisymmetric under either of these operations. This gives us the result 

We next use the fact that due to the (5-function, 

Ki + K2 = + Ki + Q Q = or27r (289) 

where Q is the umklapp momentum. It cannot be any higher multiple of 27r given that 
\Ki\ < Tx. This means 

sinQ/2 = 0,. (290) 

This result and some simple trigonometry applied to Eqn. ( |288| ) lead to the following coupling: 

^(4,3,2,1) . ,Ki-K2. . ,K^ + K2-K^-K, 

— = f/osm( )sm( )cos( ). (291) 

In arriving at this equation we have used Eqn. ( |290|) and gone from the interaction hamilto- 
nian to the corresponding action in the the path integral the usual way. The integral of —Hi 
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(with Fermi operators replaced by Grassmann numbers) from r = to r = oo becomes the 
integral over four frequencies with one overall delta function. Notice that u has no frequency 
dependence. 

H - 1285|), and restrict the momenta 



Let us now return to the general interaction, Eqn. (Ee 
to lie within A of either Fermi point L or R. Using a notation where L (left Fermi point) and 
R (right Fermi point) become discrete a label i = I or R and 1-4 label the frequencies and 



momenta (measured from the appropriate Fermi points). Eqns. ( |283| - p85| ) become 



6Sa 



E 



where 



A 

Klu 



2!2! . . , „jklo 



^ dki ■ ■ ■ dk^ f°° dui ■ ■ ■ du^ 



V^,^ (4) (3)^,, (2)^,, (4, 3, 2, 1) 



(292) 



A (27r)^ 



(27r) 



[2tt5{lUi + LU2- LU4)] 



2tx 5{ei,{KF + h)+ Si^iKp + ^2) - Si^iKp + ^3) - + ^^4)) (293) 



and 



Si = ±1 for R , L. 



(294) 



Let us now implement the RG transformation with this interaction. This proceeds exactly 
as in 0^ theory. Let us recall how it goes. If schematically 



g 4 +0>)'' 



(295) 



is the partition function and we are eliminating 0>, the effective u for 0< has two origins. 
First, we have a term —ucj)'^^ which is there to begin with, called the tree level term. Next, 
there are terms generated by the 0> integration. These are computed in a cumulant ex- 
pansion and are given by Feynman diagrams whose internal momenta lie in the range being 
eliminated. The loops that contribute to the flow of u begin at order v^. 
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Let us first do the order u tree level calculation for the renormalization of the quartic 
interaction. This gives us just Eqn. ( |293| ) with A — A/s. If we now rewrite this in terms of 
new momenta and fields, we get an interaction with the same kinematical limits as before 
and we can meaningfully read off the coefficient of the quartic-Fermi operators as the new 
coupling function. We find 

(^i' ^0 = '^uhi2n {Kl ^'il (296) 

The reader who carries out the intermediate manipulations will notice an important fact: 
Kp never enters any of the b functions: either all Kp^'s, cancel in the nonumklapp cases, or 
get swallowed up in multiples of 27r (in inverse lattice units) in the umklapp cases due to 
the periodicity of the 5-function. As a result the momentum 5 functions are free of Kp and 
scale very nicely under the RG transformation: 

l>(k) 6{k'/s) (297) 
= s6{k') (298) 

Turning now to Eqn.( |296D , if we expand m in a Taylor series in its arguments and compare 
coefficients, we find readily that the constant term uq is marginal and the higher coefficients 
are irrelevant. Thus u depends only on its discrete labels and we can limit the problem to 
just a few coupling constants instead of the coupling function we started with. Furthermore, 
all reduce to just one coupling: 

uo = ulrlr = urlbl = -urllr = -ulrrl- (299) 

Other couplings corresponding to LL RR are wiped out by the Pauli principle since they 
have no momentum dependence and can't have the desired antisymmetry. 
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As a concrete example consider the u that comes from the nearest- neighbor interaction 
Eqn.(2^) reproduced here for convenience: 



^(4,3,2,1) . ,K^-K2. . ,K3-IU, ,K^ + K^-K^-K, 

— = f/o sm( ) sm( ) cos( ) (300) 

and ask what sorts of couphngs are contained in it. 

If 1 and 2 are both from R, we find the following factor in the coupling 

. Ki - K2. . fh- 

sm = sm 



2 

.ki - k2) 



(301) 



which leads to the requisite antisymmetry but makes the coupling irrelevant (due to the 
k's). There will be one more power of k from 3 and 4 which must also come from near just 
one Fermi point so as to conserve momentum modulo 27r. For example the umklapp process, 
in which RR <-> LL, has a coupling 

UNN{umklapp) ~ {ki — k2){k3 — k^) (302) 

and is strongly irrelevant at the free-filed fixed point. 

On the other hand if 1 and 2 come from opposite sides, 

K — K 

sin(^ -) ~ sin(7r/2 + 0{k)) (303) 

2 

and likewise for 3 and 4, and we have a marginal interaction uq with no fc's in the coupling. 

The tree level analysis readily extends to couplings with six or more fields. All these are 
irrelevant, even if we limit ourselves to constant {ui and k independent) couplings. 

To determine the ultimate fate of the coupling mq, marginal at tree level, we must turn 
to the one loop RG effects. 
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E RG at one loop: The Luttinger Liquid 



Let us begin with the action with the quartic interaction and do a mode ehmination. To 
order u, this leads to an induced quadratic term represented by the tadpole graph in Figure 
5. We set a; = A; = for the external legs and have chosen them to he at L, the left Fermi 
point. The integral given by the diagram produces a momentum independent term of the 
form S/iipj^ipL. But we began with no such term. Thus we do not have a fixed point in this 
case. Instead we must begin with some term Sn^ipiipL such that upon renormalization it 
reproduces itself. We find it by demanding that 



5/1* — s 



SfjL* 



Ur 



dk 



1 



oo 27r A/s<|fc|<A 27r iu — k 



(304) 



where we have used the zeroth order propagator and the fact that to this order any uq = u^. 
The exponential convergence factor is the one always introduced to get the right answer for, 
say, the ground state particle density using < ipip >. Doing the uj integral, we get 



6^* 



dk 

lA/s<\k\<A 27r 



27r 



(305) 
(306) 



It is evident that the fixed point is given by 



(307) 



Alternatively, we could just as well begin with the following relation for the renormalized 
coupling 



Sji — s 
which implies the flow 



^ ^ I-oo 27r Ia 



oo 2n JA/s<\k\<A 271 ILO — k 



dfi Uq 
'dt ~ ^~ 27r 



(308) 



(309) 
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assuming we choose to measure /z in units of A. The fixed point of this equation reproduces 
Eqn.( 



We can find in yet another way with no reference to the RG . If we calculate the 
inverse propagator in the cut-off theory to order u, we will find 

G-^ = ruj-k-^ (310) 
2tx 

indicating that the Fermi point is no longer given by = 0. To reinstate the old Kp as 
interactions are turned on, we must move the chemical potential away from zero and to the 
value = Thus the correct action that gives us the desired Kp, for this coupling, to 
this order, is then schematically given by 

S = ^{iu -k)^ + ^^^^ + (311) 

An RG transformation on this action would not generate the tadpole graph contribution. 

A very important point which will appear again is this: we must fine tune the chemical 
potential as a function of u, not to maintain criticality (as one does in 0'^ where the bare 
mass is varied with the interaction to keep the system massless) but to retain the same 
particle density. (To be precise, we are keeping fixed Kp, the momentum at which the 
one-particle Greens function has its singularity. This amounts to keeping the density fixed 
( Luttinger(1960).) If we kept /i at the old value of zero, the system would fiow away from 
the fixed point with Kp = 7t/2, not to a state with a gap, but to another gapless one with 
a smaller value of Kp. This simply corresponds to the fact if the total energy of the lowest 
energy particle that can be added to the system, namely fi, is to equal 0, the kinetic energy 
at the Fermi surface must be slightly negative so that the repulsive potential energy with 
the others in the sea brings the total to zero. 

Now, we do not have to work with fixed density; we could take the given /i and accept 
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whatever Kp it leads to. However, we will frequently work at fixed density for two reasons. 
First, there is a simple experimental way to fix the density, namely send in the desired number 
of particles and seal the system, i.e., work with the canonical and not grand canonical system. 
(This simple and viable procedure looks rather complicated in the grand canonical language.) 
Secondly we will keep Kp fixed in many two and three dimensional cases to make contact 
with the pioneering work of Landau, done at fixed density. 

Let us now turn our attention to the order graphs that renormalize uq. These are 
shown in Fig. 6. The increment in uq, hereafter simply called m, is given by the sum of the 
ZS (zero-sound), ZS' and BCS graphs. The analytical formula for the increment in u is 

du{A?,2l) = m(6351)m(4526)G(5)G(6)(5(3 + 6 - 1 - 5)d5d6 
- j ii(6451)ii(3526)G'(5)G'(6)5(6 + 4 - 1 - 5)d5d6 
-^j m(6521)m(4365)G'(5)G'(6)(5(5 + 6 - 1 - 2)dbdQ (312) 

where 1 to 4 stand for all the attributes of the (slow) external lines, 5 and 6 stand for all the 
attributes of the two (fast) internal lines: momenta (restricted to be within the region being 
eliminated), and frequencies; G are the propagators and the 5 functions are for ensuring the 
conservation of momenta and frequencies and / dhdQ stands for sums and integrals over the 
attributes 5 and 6. (In the figure the momenta 1 to 6 have been assigned some special values 
(such as 5 = K in Fig. 6a) that are appropriate to the problem at hand. The formula is very 
general as it stands and describes other situations as well.) The couplings u are functions 
of all these attributes, with all the requisite antisymmetry properties. (The order in which 
the legs are labeled in u is important due to all the minus signs. The above equations have 
been written to hold with the indicated order of arguments. In their present form they are 
ready to be used by a reader who wants to include spin.) 
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This is the master formula we will invoke often. It holds even in higher dimensions, if we 
suitably modify the integration region for the momenta. 

To derive this formula , we do exactly what we did in Section II with bosons. We split 
the modes into slow and fast ones and do the fast integral, using the cumulant expansion to 
collect the terms that feed back to the quartic coupling. The calculation uses the fermionic 
Wick's theorem and exactly the same three diagrams (ZS,ZS' and BCS) that we saw in 
the scalar field example appear. The major difference compared to the scalar case is in 
the extra minus signs in the fermionic Wick's theorem. Of course the propagators now 
have different forms and the range on loop variables is different. Readers familiar with 
Feynman diagrams may obtain this formula by drawing all the diagrams to this order in 
the usual Feynman graph expansion, but allowing the loop momenta to range only over 
the modes being eliminated. In the present case, these are given by the four thick lines 
labeled a,b,c and d in Fig. 7 where each line stands for a region of width dK located at the 
cut-off , ie., a distance A from the Fermi points. The external momenta are chosen to be 
(4321) = {LRLR), at the Fermi surface. All the external fc's and tu's are set equal to zero 
since the marginal coupling u has no dependence on these. This has two consequences. First, 
the loop frequencies in the ZS and ZS' graphs are equal, while those in the BCS graph are 
equal and opposite. Second, the momentum transfers at the left vertex are Q = Ki — K3 = 
in the ZS graph, Q' = Ki — K4 = vr in the ZS' graph, while the total momentum in the BCS 
graph is P = Ki + K2 = 0. Therefore if one loop momentum 5 = K lies in any of the four 
shells in Fig. 7, so does the other loop momentum 6which equals K, K + tt or —K in the ZS, 
ZS' and BCS graphs respectively. Thus we may safely eliminate the momentum conserving 
5 function in Eqn.( |312] ) using / dQ. This fact, coupled with 

E{-K) = E{K) (313) 
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E{K'^K±n) = -E{K) (314) 

leads to 

_ r duodK u{KRKR)u{LKLK) p r duodK u{K' LKR)u{RKLK') 

du[ R R) - J^^ 47r2 {iu - E{K)){iuj - E{K)) ~ i-oo JdA Air^ {icu - E{K)){iuj + E{K)) 
1 r duodK u{-KKLR)u{LR - KK) 
~2J-ooJdA 47r2 {iuj - E{K)){-iuj - E{K)) ^ ' 

= ZS + ZS' + BCS (316) 

where J^j^ means the momentum must he in one of the four shces in Fig. 7. 

The reader is reminded once again that the names ZS, ZS or BCS refer only to the 
topologies of the graphs. To underscore this point, especially for readers who have seen a 
similar integral in zero sound calculations, we will now discuss the ZS graph. In the present 
problem the loop momentum K lies within a sliver dA of the cut-off. Both propagators have 
poles at the point lu = —iE{k = ±A). No matter which half-plane this lies in, we can close 
the contour the other way and the cu integral vanishes. This would be the case even if a small 
external momentum transfer {Q = K3 — Ki << A) takes place at the left vertex since both 
poles would still be on the same side. This is very different from what happens in zero sound 
calculations where the loop momenta roamed freely within the cut-off, and in particular, go 
to the Fermi surface. In that case, the integral becomes very sensitive to how the external 
momentum transfer Q = K3 — Ki and frequency transfer fl = uj^ — cui are taken to zero 
since any nonzero Q, however small, will split the poles and make them lie on different half 
planes for k < Q and the integral will be nonzero. It is readily seen that 

dudk 1 dk i 

J-00 J-A Att^ (iu;-k)iiu;-k-Q + in) ^ J-A 2^nTiq^^^ ~ ^ ^ ^ ' 

where the step function 9{k) is is simply related to the Fermi function: f{k) — 1 — 9{k). If 
we keep ^2 7^ and send Q to zero we get zero. On the other hand of we set ^2 = and let 
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Q approach zero we get (minus) the derivative of the (Fermi) 9 function, i.e., a 5- function at 



the Fermi surface . Thus reader used to zero-sound physics should not be disturbed by the 
fact that the ZS graph makes no contribution since the connotation here is entirely different. 

□ 

Now for the ZS' graph. Fig. 6b, Eqn.( |3T5D . We see that K must lie near L since 1 = R 



and there is no RR scattering. As far as the coupling at the left vertex is concerned, we may 
set K = L since the marginal coupling has no k dependence. Thus K + 7t = R and the vertex 
becomes u{RLLR) = —u. So does the coupling at the other vertex. Doing the uj integral 
(which is now nonzero since the poles are always on opposite half-planes) we obtain, upon 
using the fact that there are two shells (a and b in Fig. 7) near L and that \E{K)\ = \k\ = |A|, 



ZS' = [ 

Jd 



dK 



dAeL Att\E{K)\ 

(318) 

2n A ^ ' 



The reader may wish to check that the ZS' graph will make the same contribution to the 



/3-function in the field theory approach. 



The BCS graph (Eqn. (|315|) , Fig. 6c) gives a nonzero contribution since the propagators 
have opposite frequencies, opposite momenta, but equal energies due to time-reversal invari- 
ance E{K) = E{—K). We notice that the factor of | is offset by the fact that K can now 

lie in any of the four regions a,b,c, or d. We obtain a contribution of the same magnitude 

-'^-'^Before we move on the ZS' graph, let us notice another related fact. Suppose we choose to find the one 
loop /3-function using the field theory method. Then we will calculate the the four-point function Fin a cut- 
off theory and demand that it be cut-off independent. The same three graphs will appear in the expression 
for r, but the loop integrals will in fact go all the way up to the cut-off. Consequently the ZS graph will 
make a contribution that is very sensitive to how the small external momenta and frequencies are chosen. 
However since this contribution, when nonzero, comes from integrating a (5-function at the Fermi surface, it 
will not have any sensitivity to the cut-off and will make no contribution to the derivative with respect to 
the cut-off i.e., to the /3- function. The situation is parallel to what we saw in massless theory in Section 
IL There the expression for F had an infrared divergence when rp = 0, but that did not affect the derivative 
with respect to the upper cut-off A. On the other hand in the modern approach one never saw any singular 
behavior even in the intermediate steps. 
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but opposite sign as ZS' so that 

du = (319) 

dt 

du ^, ^ , ^ 

— = /?(u)= 0. (320) 

Thus we find that u is still marginal. The fiow to one loop for // and u is 
There is a fine of fixed points: 

u* arbitrary. (324) 



Notice that /? vanishes due to a cancellation between two diagrams, each of which by itself 
would have led to the CDW or BCS instability. When one does a mean-field calculation for 
CDW, one focuses on just the ZS' diagram and ignores the BCS diagram. This amounts to 
taking 

which, if correct, would imply that any positive u grows under renormalization. If this 
growth continues we expect a CDW. On the other hand, if just the BCS diagram is kept we 
will conclude a run-off for negative couplings leading to a state with < ipRipL >^ 0- 

What the /3 function does is to treat these competing instabilities simultaneously and 
predict a scale-invariant theory. 

Is this the correct prediction for the spinless model, which as we saw, had the marginal 
interaction in its interaction? Yes, the exact solution of Yang and Yang (1976) tells us there 
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is no no gap till u is of order unity. If the RG analysis were extended to higher loops we 
would keep getting /? = to all orders. This follows from the the Ward identity in the cut-off 
continuum model (De Castro and Metzner 1991) which reflects the fact that in this model, 
the number of fermions of type L and R are separately conserved. How do we ever reproduce 
the eventual CDW instabihty known to exist in the exact solution? The answer is as follows. 
As we move along the line of fixed points, labeled by u, the dimension of various operators 
will change from the free-field values. Ultimately the umklapp coupling , {RR LL) , which 
was suppressed by a factor {ki — k2){k^ — k^) , will become marginal and then relevant. If we 
were not at half-filling such a term would be ruled out by momentum conservation and the 
scale invariant state, called a Luttinger liquid, (Luttinger 1961, Haldane 1981), will persist 
for all u. While this liquid provides us with an example of where the RG does better than 
mean-filed theory, it is rather special and seems to occur in c? = 1 systems where the two 
Fermi points satisfy the conditions for both CDW and BCS instabihties. In higher dimensions 
we will find that any instability due to a divergent susceptibility is not precisely cancelled 
by another. 

As an aside, note that in the ZS' and BCS diagrams, the integrand is a function of just 
uj'^ -\- k"^ so that we have rotational (Euclidean) invariance. In this case we can, if we wish, 
work with a disk of radius A in a; — A; space rather than the rectangle of width A and infinite 
height. You may check that if we integrate out a shell of thickness dh. in the u — k space, 
we get the same contribution to the /3-function. 

The main results from this section are summarized in Table V. 
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V THE RG IN d > 1: ROTATION ALLY INVARIANT 
CASE AT TREE LEVEL 

We now proceed to apply exactly the same approach to spinless fermions in o? > 1. The 
nontrivial geometry of the Fermi surface will play a profound role and the application of 
the RG leads to some phenomena not seen in the usual applications to critical phenomena, 
which is of course what makes it interesting. We start with the simplest case of a circular 
Fermi surface in o? = 2. The extension of the analysis to spherical surfaces in d = 3 is very 
direct and will be explained. Only in Section X will we take up the nested Fermi surface 
that leads to CDW formation. 

A Tree Level in d = 2 

Let us begin with a square lattice containing spinless fermions at very low filling. In this 
case we can approximate the free-particle dispersion relation as follows: 

E = - cos K^- cos Ky (326) 
~ -2 + K^/2 (327) 

Since the problem now has rotational invariance, it is isomorphic to the problem of electrons 
in free space with a dispersion relation 

E=— (328) 
2m ^ ^ 

We will therefore study the latter since this allows us to make contact with Landau 's work on 
it. Let us introduce a chemical potential so that the ground state is a Fermi circle of radius 
Kp — ^/2mii. Next, we linearize the dispersion relation near the Fermi surface : 

s(K) = E{K)-n (329) 
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2m 

= '^ + 0{P) {k=\K\-KF) (331) 
m 

= vk (332) 

where v is the Fermi velocity. The free-field action now becomes 

IT TT ;^k(^^A;)(z.;-t;A;)^(u;^A;) (333) 
-oo /TT JO 2tc J~a Ztt l -I 



To obtain this, we must replace the measure Kdk by Kpdk, as the diff'erence will prove 
irrelevant under the RG and absorb a factor into each of the two Fermi-fields. Mode 

elimination proceeds just as in c? = 1: we eliminate all modes obeying A/s < \k\ < A for all 
u and 6. The same scaling of k and uj and the fields as in Eqn. (|278|) leaves 5*0 invariant. The 
only difference is that the internal index i which took just two values (L = left and R = right) 
is now replaced by a continuous parameter 6. The d = 2 theory thus looks like an integral 
over one dimensional theories, one for each direction, each with infinitesimal weight. 

Next we dispense with rotationally invariant quadratic interactions as in d = 1: either 
they modify the chemical potential, rescale existing terms, or are irrelevant. (Also irrelevant 
is the difference between K and Kp in the measure.) Let us then move on to the really 
interesting case of the quartic interaction. This has the general form 



= ^ / V^(4)V^(3)V^(2)V^(l)w(4, 3, 2, 1) (334) 

2121 JKu)9 



where 



V'(z) = %l){Ki,Ui,9i) etc., (335) 
2-K dO- /-A dk- r°° duj ' 



KojB 



n Jq 27r J -A 271 J-oo 271 



ki = \Ki\-KF. (337) 



Much of the new physics stems from this measure for quartic interactions. Let us understand 
it in some detail focusing on the factor ^^(A — \k/^\) which plays a crucial role. 

We start with a quartic interaction invariant under space-time translations and Fourier 
transform it, getting an integral over four cj's subject to a 5-function constraint and four 
momenta Ki subject to a momentum conserving ^-function. Let us now eliminate one of the 
four sets of variables, say the one numbered 4, by integrating them against the ^-functions. 
The cu- integral is easy: since all w's are allowed, the condition U4 = ui + 002 — uj^ is always 
satisfied for any choice of the first three frequencies. The same would be true for the momenta 
if all momenta were allowed. But they are not, they are required to lie within the annulus 
of thickness 2 A around the Fermi circle. Consequently, if one freely chooses the first three 
momenta from the annulus, the fourth could have a length as large as 3Kp. The role of 
9{A — \k4^\) in Eqn.( p36D is to prevent exactly this. 

Now such a ^-function will arise in the 0^ theory also if we eliminate A;4 by integrating it 
against the momentum conserving 5-function. Its effect is however quite different there. For 
one thing, even if we ignore it, nothing very serious happens since the first three fc's lie in 
the tiny ball of size A and can never stray too far, being bounded by 3A. In particular, 
it will be controlled by A and decrease with it. In the present case, even if the first three 
momenta lie on the Fermi surface, the fourth can be off by an amount of order Kp rather 
than A. Secondly, even if keep the 6 function in the 0^ case, its response to renormalization 
is very simple. Under the action of A ^ A/s, 



^(A-|A;4|) 



9{A/s-\k^\) 



(338) 




(339) 




(340) 




(341) 



8. 



.8 



Thus the 6 function of the old variables goes into exactly the same function of the new 
variables. Since the rest of the integration measure goes into itself upon rescahng from k to 
k' (and absorbing factors of s into the new fields), we get the usual result 

u'{k') = u{k'/s). (342) 

Upon doing a Taylor series in its arguments, we get the familiar result that the constant 
part Uq is marginal and the other Taylor coefficients are irrelevant in d = 4. 

Let us try to do the same here starting with Eqns.(^3^ - p37| ). We first reduce the range 



of each ki by a factor s. Then we rescale all momenta to bring the range back to the old 
value. We must finally see if the 6 function responds as it did above in the case of the 0*^ 
theory. If it does, we could conclude that 

u'{k',u',e) = u{k'/s,u'/s,e). 

But it does not! The problem is that k4 is not a function of just the other three little k's but 
also of Kp: 

^4 =\{Kf + ki)ni + {Kf + k2)n2 - {Kf + A;3)4| - Kf (343) 
where fij is a unit vector in the direction of Ki. 

fli = i cos 9i + j sin 9i (344) 

where 6i is the orientation of the unit vector along momentum K^. (In d = 2, we will use Qi 
and 6i interchangeably. ) 

It is now easy to check that if we carry out the manipulation that led to Eqn. (|341| ) we 
will find: 

eiA - Mki, k2, k3, Kf)\) ^9iA-\k'^{k[, k'^, k's, sKf)\). (345) 
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Thus the 6 function after the RG transformation is not the same function of the new variables 
as the 6 function before the RG transformation was of the old variables due to the fact that 
Kp — > sKp- As mentioned earlier, we cannot ignore the d function, since unlike in (fA 
theory, it is possible for very large /C4 (order Kp) to arise, even if the first three are of order 
A. We have a real problem implementing the RG program: how are we to say what the new 
coupling is if the integration measure does not come back to its old form? 

Before describing the solution this impasse, let us restate the problem in more geometric 
terms. Imagine that we have renormalized with a large s and are down to a shell of very 
small thickness, i.e., h^jKp is very tiny. Thus all the momenta are essentially on the Fermi 
surface and the only freedom is in their direction, or Q^. Now the point is that we cannot 
choose three of these angles freely, but only two, if all vectors are to lie on the Fermi circle. 
For example, if we choose angles 9i and ^2, the sum of the corresponding vectors lies along 
the bisector of the angle between them. The only way this initial state momentum can equal 
the final state momentum + i?4 is for the final angles to equal the initial angles: 

Case I: 63 = ^1 (346) 

^4 = 02 (347) 

or 

Case II: 63 = 62 (348) 

^4 = Oi. (349) 

In the case of identical spinless fermions Cases I and II are physically equivalent. 

There is only one exception: if the initial angles are exactly opposed to each other, leading 
to a total momentum P — 0, the final momenta are free to point in any direction as long as 
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they oppose each other: 

Case III ■.62 = -ei (350) 
^4 = -^3- (351) 

In summary either 3 and 4 are slaved to (he equal to ) 1 and 2 or 2 and 4 are slaved to (be 
opposite to ) 1 and 3. 

Let us now back off from the hmit A/Kp — and discuss the problem when A/Kp is 
small but not zero. Figure 8 depicts the situation. First let us ask for all pairs of momenta 
that 

• Lie within the annulus, and 

• Add up to some P. 

The construction in Fig. 8 gives all of them. First we draw two annuli with centers separated 
by P. They intersect in two regions (called I and II in the figure) of size of order A in each 
direction. If we start at the center of the left annulus and draw a vector to any point in I or 
II, and then a vector from that point to the center of the right annulus; we get two vectors 
that meet the twin requirements listed above. For example, the initial vector Ki and K2 
correspond to choosing this point from region II. Since + K4 — P, the latter pair must 
also stem from this construction. The figure shows them linked to region I. It is clear that 
the direction of the final vectors is within A/Kp of the initial vectors, with 63 ~ 62 and 
64 ~ 61 for this choice. (Had we chosen the final vectors to come from region II also, the 
other paring would have occurred. )The figure makes it clear that if choose the orientations 
of the first three vectors without paying attention to this restriction, i.e., choosing a point 

— * — * 

outside regions I and II (but inside the annulus) for terminating Ks, K4 can end up being 
much shorter or longer than Kp. 
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The problem of the changing ^-function may be stated as follows in terms of this figure. 
Let us begin with some cut-off and a choice of four angles that is allowed at that cutoff as per 
this construction. There is some coupling u for this choice of angles. If we reduce the cut-off, 
the allowed choices of angles shrinks and a coupling that was previously allowed may no 
longer allowed. Since the angles 6 play the role of internal (isospin- like) degrees of freedom, 
we have a situation in which the range of internal labels is changing under mode elimination. 
This makes it impossible to make meaningful comparison of the couplings before and after 
since certain processes are no longer allowed. (Imagine an SU (4) theory renormalizing to an 
SU{3) theory.) 

Since we do not want the momentum or internal labels to change their allowed range of 
values (if we are going to follow the standard RG procedure of comparing the couplings before 
and after to see how they are flowing), we will take the view that all 6's are allowed, but that 
some u 's are abruptly renormalizing to zero as the cut-off is reduced. 

The situation is a lot clearer if we use a smooth cut-off for ^4: 

^(A- |A;4|) ^ e-l'^*!/^ (352) 

for now no process is disallowed, but only exponentially suppressed as we renormalize. We 
may interpret that as the exponential decay of the corresponding u under renormalization. 
Q Now recall 

ki = \KF {Qi + Q2- Qs) +^1^1 + k2n2 - fcsOsI - Kf (353) 

v ' 

A 

where is a unit vector in the direction of In what follows we shall keep just the A piece 
and ignore the 0{k) terms. This is because the only time the latter are comparable to the 

former is when both are of order A, in which case k^ ~ Kf, and this regions is exponentially 

-"^^Nothing is gained by using a soft cut-off for the rest. 
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suppressed by the smooth cut-off Eqn. ( p52| ) anyway. 
Under the RG transformation at tree level, 

TL'^J-A ZTT JO iVT J~oo ZTT 

p <^^-[sK,/m^\-i\^^k'i, J I, e)H^ij^. (354) 

^J-A 2tt Jo 27r J-oo 27r 
Let us write 

^-{sKf /A)\\A\-1\ ^ ^-(Kp/A)\\A\-l\^-as-l)KF/A}\\A\-l\ (^^^^^ 

SO that the measures before and after have the same factor e"*^^-^/^-"'^'"^'. Now that the 
measures are identical before and after the RG transformation we can compare apples to 
apples and identify the new quartic coupling: 

u'ik'uj'e) = e-«^-i)^^/^)ll^l-^lM(A;Vs u'/s 9). (356) 

We may conclude that the only couplings that survive the RG transformation without any 
decay correspond to the cases where 

|A| = \VLi + - = I (357) 

In (i = 2 this equation has only three solutions: 

Case I : = Qi {Hence Q2 = ^a) (358) 
Case II : = (Hence Cli = ^4) (359) 
Case III : Qi = {Hence fig = -^4) (360) 

This result which was anticipated earlier, can be deduced algebraically or seen geometrically 
by considering Fig. 8 in the limit of zero shell thickness since this gives the allowed region in 
the hard cut-off scheme after infinite amount of renormalization i.e., at the fixed point. 
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Notice that we have an extra conservation law at the fixed point KjKp = 0: not only 
is total momentum conserved, the set of individual momenta is also conserved. The only 



exception is when the incoming momenta are equal and opposite. Then they add up to zero 
and the final momenta can be any two opposite points on the Fermi circle. In all cases, two 
of the momenta determine the other two rather than three of them determining the fourth. 
For couplings that do obey this condition, Eqn.( |356] ) becomes 



Performing a Taylor expansion in k and u and comparing coefficients of separate pow- 
ers, we conclude that the leading term, with no dependence on either variables is marginal, 
while all the rest are irrelevant. We shall refer to this term as u hereafter. 

We see that the tree level fixed point is characterized by two independent functions and 
not a handful of coupling. They are 



uMe,e29s)=92;03 = 9i;02;0i] = ^2) = m - ^2 = M (rot tnv), (362) 



uM9,929s) = -9,;93;-9,-9^;] = V {9,- 9^) = V {9, - 9^ = 9,^) {rotznv). (364) 



Note that the manifestation of the Pauli principle on F and V is somewhat subtle: F will 
not be antisymmetric under 1^2 since the way it is defined above, we cannot exchange 

1 and 2 without exchanging 3 and 4 at the same time. On the other hand, since 3 and 4 

^■^A common concern frequently expressed by those familiar with many-body physics of these fermions is 
this: is there not some ambiguity on how the limit a; —> is to be taken as ki —^ 7 If so, how does u even 
have a Taylor expansion at the origin of ivk space ? The answer is that u is just the bare coupling that goes 
into the cut-off theory. It is not the full four-point function T which is calculated in the cut-off theory by 
summing all Feynman diagrams with u as the coupling, A as the cut-off, and propagators whose momenta 
can go right up the Fermi surface . Indeed the T so computed has all the above mentioned singularities. In 
contrast u is obtained by taking some analytic (short range) interaction in the full momentum space and 
eliminating the modes outside the cut-off. This procedure cannot produce any nonanalyticity. The situation 
is just as in (j)'^ theory where the bare couplings in the action are analytic functions of momenta while the 
full Green's functions are plagued with infrared singularities coming from the soft modes. 




(361) 




Pauli , 



(363) 
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can be exchanged without touching 1 and 2 in the definition of V, V must go to —V when 

A concrete example is useful here. Let us begin with the nearest neighbor interaction of 
spinless fermions in d = 2, transcribed into momentum space, 



2!2! 



sm( ) sm( ) cos( ) + x ^ y 



(365) 

and expand it to second order in (as we did in the kinetic energy) obtaining a rotationally 
invariant interaction: 

M ~ (i?i - i?2) ■ (i?3 - i?4) (366) 
If we evaluate this with all \K\ at Kp and use Eqns. (p62| , p63| , p64| ) we obtain 



^(^12) ~ f/o|i?i-i?2|'-f/o(l-cos^^i2) (367) 
V(^13) ~ [/oi?i-i?3-t/ocos^i3 (368) 

Observe the effect of particle exchange on these functions. 

Let us understand the significance of the coupling function F . If we calculate the four- 
point vertex F, it will be given to lowest order by u. Since the only m's that survive require 
that the final directions ^3 and ^4 equal the initial ones 9i and 62 up to a permutation, 
only forward scattering exists at the fixed point. The only exception occurs when the 
incoming particles have opposite directions: then they can scatter to another pair with equal 
and opposite momenta, the corresponding amplitude being given of course by V . 

When we study noncircular Fermi surfaces, we will find yet another coupling function W 
that survives if there is nesting. This will correspond to processes in which the momentum 

transfer equals the nesting vector Q^. 

^"^This is true at lowest order when T = u. But higher orders in the forward scattering couphng can only 
give forward scattering. 

95 



B Tree Level analysis in d = 3 

Let us repeat the preceding analysis in ci = 3 where the Fermi surface is parametrized by two 
angles 6 and (p. All integrals over 6 get replaced by integrals over solid angle. The process 
of mode elimination and rescaling of fields and momenta k proceeds exactly as in c? = 1 or 
d = 2 since the coordinates on the Fermi surface play the role of an internal variables (like 
isospin) and are unaffected by the RG transformation. We will then end up with Eqn.(^S^). 
The equation for the couplings that survive is still the same as Eqn.( |5F7D , but now the unit 
vectors f2j can point anywhere in three dimensions. This in turn means that the condition 
I A| = 1 not only has the solutions given in Eqns.( p58| , |359| , |360| ) but a continuum of others. 



First consider Eqns.(P^- 1359|) which tell us that ^3 and ^4 must coincide with and ^2 up 



to a permutation. In d = 3, the former can rotate about their sum. In other words, Figure 8 
with zero shell thickness must now be viewed as depicting not the intersection of two circles, 
but two spheres. Thus the vectors 3 and 4 do not have to coincide with 1 and 2 (up to a 
permutation) but can rotate around a cone with opening angle equal to that between 1 and 
2. The planes containing 1 and 2 can have an angle 0i2,34 with the plane containing 3 and 
4. (The cases I and II considered in d = 2 correspond to 0i2;34 = or vr.) 

Let us review this. The incoming particle momenta 1 and 2 lie on the Fermi sphere. Their 
sum lies in the plane they define and bisects the angle between them. The final particles also 
on the sphere, can give the same sum by lying anywhere on the cone generated by rotating 
the incoming pair around their sum. 

Although the individual momenta are no longer conserved, we have the additional re- 
quirement that the angle between the final pair is the same as the angle between the initial 
pair: 

Qi-n2 = n3- Qi. (369) 
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Once again the only exception is when the incoming momenta add up to zero. In this 
case the final momenta are free to point in any direction as long as they are mutually 
opposite. Thus couplings corresponding to non-forward scattering (in which the initial and 
final directions are not the same) do not vanish under the RG transformation, but survive 
as marginal couplings. The fixed point is then characterized by a function 

F = F(Z12,012;34) (370) 

= (371) 

of two variables: 

and 0i2;34 which is the angle between the planes containing 12 and 34 respectively. 

In addition to this we still have the V function coming from Eqn. (|360|) with V = V{Qi ■ 

If we go back to nearest neighbor coupling evaluated in this case we find: 

F = [/o(l-zi2)cos0 (372) 
V = Uoz,, (373) 

Notice that since 3 and 4 are not slaved to 1 and 2 anymore, we can exchange just the latter. 
This causes — > + vr, which in turn changes the sign of F. As for V , it changes sign under 

Zl3 —^13- 

Table VI and VII contain a list of fixed points, couplings, and their fiows including ones 
to be discussed later in the paper. 
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VI RG IN d > 1 : ROTATIONALLY INVARIANT 
CASE AT ONE LOOP 



We have found that the couphngs F and V are marginal at tree leveL The next thing to 
do is to buckle down, as we did in d = 1, and go to the one loop graphs and see if they tilt 
the marginal ones towards relevance or irrelevance or preserve marginality. We first do the 
analysis m. d — 2 and then discuss briefly the changes encountered upon going to d = 3. 

A RG for the two-point function 

Just as in d = 1, mode elimination produces at one loop (tadpole graph as in Fig. 5) a nonzero 
change in the quadratic term in the form of a chemical potential. To retain the old Fermi 
surface, we must find an input 5^* which will reproduce itself under the RG transformation. 

Carrying out the same analysis as in o? = 1 we find 



where the last equation signifies that is a constant independent of a; , or ^. 

We must begin with this quadratic term if it is to reproduce itself under mode elimination. 
We may also see this as the counter term one must add as the interaction is turned on to 
maintain the same Fermi surface. The same holds in c? = 3: 




(374) 



(375) 



5n* 



(376) 




(377) 



(378) 



Note that only the forward scattering F ( with = 0) enters this equation. 
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B The one loop /3-function for F 

Now we turn to the real issue: the renormahzation of the quartic couphngs. Since some of 
these couphngs (caUed F and V) proved to be marginal at tree level, we must go, as we did in 
c? = 1, to one loop to see if they flow. The analysis will be done using the modern approach, 
though a passing remark may be occasionally made about the field theory approach. The 
same diagrams that appeared in = 1 will appear, the integrals will be of the same form as 
Eqn.(RTH), but the range of loop momentum integration will reflect the higher dimensionality. 



Eqn.( |312D gives in this case (suppressing vector symbols) 

^ r f dtodK f^- deuiK + Q,3,K,l)ui4:,K,2,K + Q) 



oo 

oo 



dA 47r2 Jo 27T {iuj - E{K)){tuj - E{K + Q)) 

dujdK r^^ de u{K + Q' , 4, K, l)u(3, K,2,K + Q') 
-ooJdA 47r2 io 2^ {iuj - E{K)){iuj - E{K + Q')) 
1 r dtodK r^^ dO u{P - K, K, 2, 1)m(4, 3,P-K, K) 
~2J-ooJdA Att^ Jo 2^ {lu - E{K)){-iu - E{P - K)) ^ ^ 

= ZS + ZS' + BCS (380) 

where Q = — Ki, Q' = K4 — Ki and P = Ki + K2. The subscript dA on /^^ reminds us 
that both loop momenta must be in the shells being eliminated. Although K is restricted 
to this region explicitly, the restriction on the other momentum, K + Q in the ZS graph; 
K + Q' in the ZS' graph and P — K in the BCS graph is implicit. The Feynman diagrams 
are shown in Fig. 9. 

On the left hand side of the above equation, we choose the variables 4321 such that the 
corresponding coupling survives renormahzation at tree level, which means it must be an F 
or a V. 

Let us begin with the renormahzation of F. We set all external legs at zero frequency 
and on the Fermi surface {k = u = 0) since the dependence on these variables is irrelevant. 
As for the angles, we chose fli = lis. Consider the ZS diagram in Fig. 9 given by the first 
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integral in Eqn.( |380| ). Since Q = here, we know that if K lies in the shell being eliminated, 
so does K + Q for any direction of K. In other words, 6 runs over the full range. However, 
this diagram vanishes for the same reason it did in the d = 1 case: both poles in the u plane 
are on the same side. Even if we put in a small momentum transfer Q << A at the left 
vertex, it will not change anything. This is because the loop momenta are near k = ±A and 
it takes a minimum momentum transfer of order A to knock something from below Kp to 
something above it and vice versa. This is in accord with the general statement that the 
bare coupling has no singular limit at small u or Q. 

Now consider the ZS' diagram. Due to the momentum transfer Q' of order Kp at the 
left vertex, not only is the magnitude of the loop momentum restricted to lie within the shell 
being eliminated, its angle is also restricted to a range of order dKjKp. This is clarified 
by Fig. 10. The dark circles now represent the thin slices being integrated. The intersection 
regions, of order rfA^, give us the allowed loop momenta for the (5 function calculation: they 
lie in the shell and have the right momentum transfer Ki^ — Ki. Of the eight intersection, 
only the four marked ones satisfy the condition of being on opposite sides of Kp so that the 
uj integral survives. Since the uj integral gives a denominator of order A, the contribution to 
du is order so that the /3-function vanishes as we take the limit dt = 0. 

It is also clear that if Q' = and Q is not, (which obtains if = Qi) we can repeat the 
above argument with ZS ^ ZS'. 

You may check that in the field theory approach we would get a /?- function that went 

-"^^ What if we choose to follow the field theory approach? The situation is exactly as in d = 1 . We must 
evaluate the one loop graphs contributing to F which is sensitive to whether we take the the Q << uj ot 
vice versa. The limit which gives a nonzero F , {uj << Q), makes a contribution of the form J duj —> 

f{K)~f{K+Q) ^Yieie f is the Fermi function. As Q ^ 0, this becomes a S function at Kp and makes a 

E{K)-E(K+Q) 

contribution that is A- independent and hence irrelevant to the /3-function. 
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as Since we are ultimately going to send 7^ ^ 0, this will not matter in that limit. But 
note that the two schemes do not agree in detail except for the marginal flows. 

Finally for the same kinematical reason, the BCS diagram does not renormalize F at one 
loop. Consider Fig. 8 with and replaced by the two momenta in the BCS loop, K and 
P — K. In each annulus keep just two shells of thickness dh. at the cut-off corresponding to 
the modes to be eliminated. The requirement that K and P — X lie in these shells and also 
add up to P forces them into intersection regions of order dis?. This means the diagram is 
just as ineffective as the ZS' diagram in causing a flow. Thus any F is a flxed point to this 
order. Notice that unlike the d—1 /3-function which vanished due to a cancellation of two 
terms, the present one vanishes trivially. 

C The one loop /3-function for V 

Let us now look at the evolution of V . We choose the external momenta equal and opposite 
and on the FErmi surface. The ZS and ZS' diagrams do not contribute to any marginal 
flow for the same reason BCS and ZS' did not contribute to the flow of F: since Q and Q' 
are of order Kp, they are kinematically suppressed by an extra factor in the fleld theory 
approach and by in the modern approach. But the BCS diagram produces a flow in 
either approach. We follow the modern approach. The flow is due to the following factors: 

• The loop angle can run freely over its full range because no matter what value K takes 
in the shell being eliminated, the other momentum P — K — —K automatically lies in 
the shell. The two energies E{K) and E{—K) are equal while the two a;'s are equal 
and opposite. 

• The u and k integrals behave as in d = 1 and produce a factor dt — 
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We find 

dViOi - 63) 1 /-s- de 



/ —^iei-e)vie-es) [ssi] 

Jo ivr 



dt 87r2 Jo 2tt 

This is an interesting example of tlie /3-function for tlie coupling function. Fortunately we 
can simplify the picture by going to angular momentum eigenf unctions: 

Vi = / ^e^'V(^) (382) 

JO iTT 

and obtain an infinite number of flow equations, one for each angular momentum /: 

The flow tells us that the couplings Vi , marginal at tree level, are marginally relevant if 
negative, and marginally irrelevant if positive. The former is just the BCS instability. As 
for the repulsive case, if we integrate the flow we get 

Viit) = , , (384) 

l + t[V,(0)/47r] ^ ^ 

which is just the observation of Morel- Anderson (1962) that inter-electron repulsion is loga- 
rithmically reduced if we develop an effective theory for the modes close to Kp- 

If this analysis were to be repeated in c? = 3, the only difference would be that the BCS 
/3-function for V{z) would be decoupled using the Legendre polynomials Pi{z) with I odd. 
The decoupled equation would have the same form as Eqn.( p8^ ) with the same implications. 

D Fixed point structure at one loop 

Let us take stock: the tree level fixed point is characterized by two marginal functions F 
and V . The function F is marginal at one loop also, while V is marginally relevant in an 
infinite number of ways, one for each angular momentum /, if attractive, and marginally 
irrelevant if repulsive. It appears that even if a single Vi is negative, we run off to some 
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other massive fixed point with a BCS gap. What if all the VJ's are positive? It appears that 
these couplings will renormalize to zero logarithmically and we will end up with a fixed point 
characterized by F . This turns out be incorrect, atleast in principle, due to a reason first 
pointed out by Kohn and Luttinger (1965). This is tied to some irrelevant operators which 
cannot be ignored. Here is the point. An irrelevant operator by definition is something 
that renormalizes to zero, not something you can set equal to zero at the outset without any 
consequences. Before it renormalizes to zero, it can modify the flow of the relevant couplings. 
Recall the case of mq, the scalar coupling in = 4. Although it was irrelevant, adding it to 
the gaussian fixed point generated a mass term tq which then took off. A very similar thing 
happens here: an irrelevant term produces a small negative BCS coupling , which then takes 
off. This subtle issue is discussed in the next subsection. 

We discuss here another irrelevant term which does not destabilize the fixed point, but 
modifies our description of it. 

Consider the sunrise diagram. Fig. 11. In mode elimination this diagram comes from 
taking two quartic terms and seeing how they feed back on the quadratic term. Though it 
is also of order v?, it has two loops and may be ignored in the one loop discussion we are 
having. But if one evaluates it, one finds it is irrelevant due to phase space restrictions in 
the limit of small K/ Kp- However, before we reach this limit, it can produce interesting 
effects that the one loop analysis did not include. 

Let us write its contribution as —Tj{kuj)ipip where E is called the self energy. ( It has no 
dependence here due to rotational invariance.) If we Taylor expand S as follows: 



The details of the evaluation will not be provided here. The interested reader is asked to do the phase 




(385) 



as 

diuj 



space analysis. 
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we see that 



E(00) is to be eliminated using a counter-term S/i of order li^.This is to ensure that 
Kp is unchanged. 

• (1 — Z~^) changes the coefficient of icu to Z~^. 

• 1^ changes the coefficient of k from Kp/m to (1 + ■^^)KF/m. 

At the free -field fixed point we rescaled the field to keep the coefficient of both the quadratic 
terms fixed. Now we have seen that there is a manifold of fixed points parametrized by 
nonzero u. In this case iu) and k may receive different contributions from eliminated modes 
and there is no rescaling which will keep both terms fixed. So we will keep the coefficient of 
iu fixed at unity , i.e., define 

ij;' = s-3/2z-V2^ (386) 

This has two effects. The coefficient of k changes as follows: 

Kp m Kp 

m Kp ok m 



(387) 



m* 

which defines the effective mass m*. 

Next, the new quartic coupling is given by 

u'={u + 5u)Z'^ (388) 

where 5u is the contribution we have already discussed. Question: does this modify the {3 
function we calculated? Answer: not to order v? since Z deviates from unity at order v? 
and this produce changes of order v? in the equation above. 

How about the fact that m* is now moving as we renormalize ? Upon looking at the 
kinematics of the sunrise diagram one can tell that as A goes to zero, its contributions will 
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vanish. Thus although m will evolve to m* in the early stages of renormalization, a fixed 
point characterized by some m* will emerge asymptotically. 

How can we have a fixed point in a theory where there is a nonzero dimensionful parameter 
Kpl The answer is that the fixed point theory described above has no knowledge of Kp: 
the ZS graph gets all its contribution from the delta function at the Fermi surface, while 
the other two, which know about Kp, are suppressed by the factor A/Kp and vanish at the 
fixed point. This situation is not changed by going to higher loops, as will be shown in the 
next section on the expansion. 

E The Kohn-Luttinger Effect 

The fiow of Vi was such that attraction led to instability, while repulsion meant downward 
renormalization to 0. There is no doubt the former is true, but the latter is incorrect in 
principle. Now the fault is not with the solution to Eqn. (|383|) , but with the equation itself. 
In deriving it, we ignored the contribution from the ZS and ZS' channel graphs on the 
grounds that they were finite and down by powers of A/Kp, which made them unimportant 
in the limit A/Kp 0. But there is a surprise waiting for us if we go ahead and compute 
their contribution to the fiow. As shown in the Appendix , the modified fiow in c? = 3 is 
(upon setting all positive numbers independent of / and t to unity): 

''"i'^)^'^^ (389) 

dt An /15/2 ^_X7/4 + /-7/2]2 ' ^ ^ 

where V{7i) is the backward scattering amplitude in the BCS channel and A = A/Kp. 
Notice that as A/Kp — > 0, the second term vanishes as {A/KpY^^ which is why we ignored 
it. (The strange power comes in because the intersection region scales in a special way when 
the momentum transfer is ~ 2Kp, which is the region that dominates this piece; note the 
backward scattering amplitude in the answer.) Why do we care about this piece? 
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Let us imagine that we are just beginning our RG . The input potential is the projection 
on the angular momentum I channel of some short range potential. It follows that V;(t = 0) ~ 
e~' as Z — > oo in order that the sum of such coefficient times the Pi{z) and all the derivatives 
of the sum with respect to z converge to given analytic function V{z). By contrast, the 
second term, at fixed A goes as as Z — > oo. It is clear that as soon as the flow begins, 

the exponentially small initial coupling V;(0) will very quickly be driven to negative values 
by the second term. Thereafter both terms will drive the instabihty. 

This is the RG version of the famous Kohn-Luttinger argument, which is discussed at 
some length in the Appendix . The argument implies that at T = 0, the flxed point we 
studied always faces the BCS instability. However one is still interested in this fixed point 
characterized by F. The reason is that for the Kohn-Luttinger superconductor could 
be very low, or the I value for pairs, absurdly high. Thus if we imagine a tiny nonzero 
temperature above this Tc, the instability will disappear. (For a recent survey of the Kohn- 
Luttinger effect see Baranov et al 1992).) Recall from our analysis in Section III that a 
temperature T = 1/P leads to an imaginary time coordinate of range < r < /3. In other 
words a quantum system infinite in space and at inverse temperature (3 is mapped by the 
path integral method to a system in d+1 dimensions which is infinite in the spatial direction 
and of width (3 in the (imaginary) time direction. As we renormalize, the thickness in the 
time direction (in the sliding units) will get reduced by s just like the correlation length, and 
just unlike the momenta. Thus we will fiow to smaller /? or larger T, the fixed point being 
T = oo. But if the cross-over is very slow, then in the interim the fixed point described by F 
will control the physics. Interestingly enough there are many real world systems described 
in exactly that fashion. More on this in the section on Landau's Fermi Liquid Theory Table 
VI contains a summary of results from the one loop analysis. 
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VII THE 1/N PICTURE, LEAP TO ALL LOOPS 

So far we have followed the RG program to one loop in some detail. We must now see 
what happens at all loops. In general this would be a formidable problem. Luckily in the 
present problem a great simplification arises in the limit — > 0. The presence of this small 
parameter allows us to relate the sum over all loops to the one loop result. 

Now a similar thing happens in theories where N, the number of species of fields, becomes 
very large. In this case, in the limit — > 0, it is possible to sum over all loops and the 
answer is expressed in terms of the one loop graphs. It is also possible to correct this answers 
in powers oi 1/N. This is called the 1/N expansion. 

We will begin with a review of the approximation. It will then be shown that A/Kp 
plays the role oi 1/N. We begin as usual in d — 2. 

A 1/N in d=2 

Consider a 4>^ theory with action 

N ,2 1 JV AT 

S^--ll^t^^r-l^Y.Y. 4>l4>^V^,ct>]ct>,. (390) 
i=l ^ i=lj=l 

Lots of integrals are suppressed and only the internal index is highlighted since what we are 
about to say is independent of dimension. All we need to note is that there are N species 
of fields (or particles) and they have a quartic interaction V . The interaction has a factor 
1/N \n front of it to ensure that we have a nontrivial limit as A?" — > oo. Note that in the 
interaction vertex, if an i and a j come in, the same indices also exit as shown in Fig. 12. 

Let us look at a four-point function to one loop, as shown in Fig. 12. Among the one loop 
graphs, only the first is of the same order as the tree level graph: this is because the extra 
factor oi 1/N coming from the extra vertex is compensated by a sum over the loop index 
which is free to take all values. This in turn was because the index i that came in went 
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out at once, leaving the loop index I free to roam over all values. By contrast, the external 
indices have insinuated themselves into the loops in the other two diagrams (ZS' and BCS) 
and there is no sum over the indices there. These graphs are then of order and hence 

suppressed by a factor 1/N relative to the tree graph. It is clear that the sum over iterated 
ZS loops ("the bubble sum") gives the leading behavior (in of the four-point function. 

Note that the /9-function of this theory is completely given by the one loop answer. This 
is because iterates of the ZS loop ("the bubble sum") merely produce the n-th power of In A 
when iterated n times while higher order terms in (3 come from subdominant logarithms. 
This conclusion is also evident from Weinberg's (1993) discussion o the graphical content of 
the effective action. 

Let us look at our theory now and consider a four-point function r(^2, ^i, ^2, ^i) in obvious 
notation with all external cu — k = 0. Return to Fig. 12, this time using integrals over ^'s 
in place of a sum over discrete indices. In the ZS graph, where the incoming "index" 9i 
immediately exits, there is no momentum transfer at the left (and hence right) vertex. Thus 
the loop angle runs over the full range — 27r. In the ZS' and BCS diagrams on the other 
hand, there is a large momentum transfer Q' or large total momentum P (Figures 8 or 10) 
tell us that loop angles must he within A/Kp of the external angles. In other words, the 
external angles have insinuated themselves into the loops and frozen the loop sum. Since 
the ratio of the ZS' and BCS graphs to the ZS graphs is A/Kp, we expect that this ratio 
will play the role oi 1/N. 

It should not be too surprising that we have a 1/N description available for us here: 
the noninteracting d = 2 theory was written as an integral over the internal index 6 which 
labeled pseudo two-dimensional theories, (with a phase space dkdcu), one for each direction. 
But what should be the N assigned to this integral (rather than sum) over 97 Is it infinity. 
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since there are infinite directions or is it of order unity since each has infinitesimal measure? 
Let us sharpen the analogy with 1/N analysis to answer this question. As a first step let us 
write the free-field action with all factors of Kp intact: 

5'o = / — / — Kp / — i){i0J - v*k)il) (391) 



-A 2-K Jo 2tt J-oo 2tt 
Now chop the angular integration into regions of width = 2A/ Kp so that the annulus 
breaks up into N — 2t:Kp/2A. cells labeled by an index i. The momentum of a point within 
a cell i is 

K = Kp Qi + kiQi + ku U (392) 
= Kpni + ki (393) 

— * — * 

where is the unit radial vector at the center of cell i, is a unit tangent, ki and k±i are 

radial and angular displacements from the cell-center. We refer to the first piece of order 

Kp as the large momentum and the other as the small one. The measure per cell is 

A dk r^/KF ^ de _ /-A dk dk^ _ r d^k 
-a2^J-a/Kp ^2^ ~ J-a2^ J-A^ " J (27r)2' ^ 

In this notation 

N r ^ ^ d'^kduj- 

So^Ylj i^iikiLOi) [iuj - v*ki\ ipi{ki ^i)-^^^- (395) 

We now write down an interaction term 

Sp = -— / i^jih^^d'^jih uj2)Fiji/ji{k3UJ3)ipi{ki uji) (396) 
Fij = FiQi-n^) (397) 
/ = / n ^^^|^27r%i + - - (398) 

Notice that the interaction conserves momentum: the fact that the index i appears once in 
a ijj and once in ijj means the large momentum is conserved. The small one is also conserved 
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because of the explicit S^'^^ function. If we express all frequencies and momenta in units of 
A, it will be found that the only place where A appears will be front of the interaction term 
in the form A/Kp exactly playing the role of The proof is left as an exercise. 

Now, this interaction is not exactly the one we have dealt with so far; it forbids certain 
processes that were previously allowed. Consider a process where the initial particles have 
momenta Ki K2 from cells ii 12 respectively adding to a total P. If we draw Figure 8, we 
will indeed find that K4 lie with the same intersection regions as Ki K2 , hut this need 
not mean they are in the same cells. In other words, the intersection region can straddle 
more than one cell, if we imagine these cells permanently etched on the annulus. The old 
interaction would allow the full intersection region, while the new one would only allow the 
part of it in which the cell indices match in pairs. But notice that the difference between the 
two interactions shows up only in situations which are kinematically suppressed by a power 
of kjKp. On the other hand, if we concentrate on just the sum over iterated ZS loops, with 
cu, Q << A, it can be verified that the two give identical answers. 

The careful reader will find one problem with the analogy. In the usual examples, 
u in Eqn.( p90D is held fixed as varies. Here A^ is related to the cut-off A. As we lower 
the latter to increase the former, we must follow the flow of -F(A) = F(l/A^). However, we 
have seen that there is no flow within the F couplings in the asymptotic region. Thus Fij is 
essentially constant independent of A^ for large A^. 

In view of the above, we make the employ the following two-stage attack on the fermion 
problem: 

• Reduce the given problem in the full K-space to a small A (large N) theory using the 
RG . 

• Solve the resultant theory using the smallness of 1/A^ . 
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Consider now another type of coupling corresponding to . In schematic form and the 
same notation used above, 



This interaction leads to a bubble sum (iterates of the one loop graph) in the BCS channel. 
Here the coupling V grows as increases since there is a flow now. If we want to increase 

keeping V^(A^) fixed, we must start with weaker and weaker Vs. If we do this, we can 
handle the BCS problem also by summing over bubble diagrams. 

Let us hereafter assume the BCS amplitudes V are absent in the spirit discussed earlier. 
Then we are left with iterated ZS loops and F's. The physics of this is Fermi liquid theory, 
to be discussed in the next section. 



We have seen that RG allows nonforward amplitudes to survive in d = 3 as marginal in- 
teractions. If we divide the spherical Fermi surface into patches of size A^, and label them 
by an index i, this will run over roughly {Kp/KY values. The interaction Fij will not be 
of the form Eqn. ( |398| ) since non-forward scattering amphtudes are allowed. Let us however 
divide the possible interactions into a set involving just forward scattering and the rest. If 
we consider a forward scattering four-point function F*^ (with cell index conserved) it will 
be given as the bare vertex plus a sum of iterated ZS diagrams all involving only forward 
scattering amplitudes. The insertion of a large angle scattering anywhere will produce a 
factor A/Kp. (The suppression factor is not [A/Kp)'^ because kinematics only restricts the 
the angle between 1 and 2 to be that between 3 and 4, the plane containing the latter is still 
free rotate by the angle we called 0i2;34- )If we consider instead Ffj, which is nonforward, 
we will find that it is given by just the tree level coupling we put in the action. All loop 




(399) 
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1/N ind = 3 
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corrections will be down by powers of KjKp. If we consider any response function to a soft 
(low wave number) probe, these amplitudes will not enter the physics. If we compute the 
lifetimes of particles, these amphtudes will play a role, but phase space will again introduce 
powers of KjKp. Thus the nonforward amplitudes live in shadow world, perhaps large in 
magnitude but small in their effects. 

As for the BCS amplitudes, we can find the stable state with a gap, using the 
expansion to limit ourselves to summing bubble graphs. Since nonforward amplitudes never 
enter the computation, this explains why it is permissible to use the reduced BCS hamiltonian 
in which only scattering within pair states is kept. (These amplitudes are just our V 's of 
course.) 

Note that as in any theory, the /5-function at one loop (which we have calculated) 
is all there is in the large or small A limit. 

C Two-point functions at large N 

Let us close by asking what happens to two-point functions at large A^. Let us look at 
some of the graphs contributing to shown in Fig. 13 if we use the one loop fixed point 
action Eqn.( p7^ ). We see that all iterates of the tadpole are exactly canceled by the one loop 
counterterm or fixed point chemical potential whenever we can draw one more loop, 
we can also use the counter term which exactly kills it. The sunrise diagram on the other 
hand, brings in new corrections. But this diagram is suppressed by 1/A^ and ignored in the 
limit we are interested in. This means that in this theory, we know the exact self- energy of a 
particle for the given interaction F. This will play a big role in Fermi liquid theory to which 
we will turn our attention shortly. (One way to understand the above result is the following. 
When the shell thickness goes to zero, the limiting theory acquires particle-hole symmetry. 
We are quite accustomed to this symmetry determining the chemical potential, even on a 
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lattice. This will continue to be true even for nonspherical Fermi surfaces. ) 

The picture we developed in the large limit agrees with rigorous calculations of Feldman 
et al (1990,1991,1992), as explained to me by E.Trubowitz. According to these authors: 

• The system always goes to a BCS state at T = 0. 

• If the BCS diagrams are eliminated the rest define a convergent series with a finite 
radius of convergence. 

Further details may be found in the references given. 

VIII LANDAU'S FERMI LIQUID THEORY 

Nearly four decades ago Landau (1956, 1957, 1959) attacked the problem of interacting 
fermions at very low temperatures T « Kp. Assuming the system evolved continuously 
from the noninteracting limit, he developed a phenomenological theory which proved very 
successful, for example in the study of Helium -3 (Pines and Nozieres (1963), Leggett (1975), 
VoUhardt and Wolfle (1990), Baym and Pethick (1991), Lifshitz and Pitayevskii 1980.) The 
picture he arrived at, called Fermi liquid theory , may he described in the terminology of 
this paper as the the fixed point described by F . For many readers of Landau's work there 
was an element of mystery surrounding some of the manipulations. This had to be so, 
since he substituted forty years of subsequent developments (the RG in particular) with his 
remarkable intuition. 

Following his work, a diagrammatic derivation of Fermi liquid theory was provided by 
Abrikosov et al and is described in their book (Abrikosov al (1963)). While the details were 
rather tedious, they did a lot to clarify where everything came from. I believe the approach 
developed here, using the RG , provides an even simpler route to Fermi liquid theory at least 
for those with a certain background. 
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So let us pretend we do not know what Fermi liquid theory is and ask how we would 
arrive at it. (This paper specializes in T = 0, and one expects the results to work also 
for the crossover region at low T with minor changes.) Using the RG developed here, we 
would find that after repeated renormalization, we would have mapped the initial problem 
to one with K/Kp — > 0. Setting V — Q eliminates the BCS instability, i.e., implements 
Landau's requirement that there be be no phase transitions, leaving us with a fixed point 
theory parametrized by the marginal couplings m*, F{z(f)) . The physics of this fixed point is 
Landau's theory. The excitations of this effective theory are the quasi-particles of Landau, in 
contrast to the "bare" particles created by the fields we began with prior to mode elimination. 
The fact that the quasi particles have infinite lifetimes was established by Landau using phase 
arguments. In the present analysis, the lifetime terms, which appear as 0{u!^) terms in the 
self-energy, are irrelevant under the RG transformation. (At the free field fixed point, they 
fall down by under a factor s reduction of the cut-off.) 

Is there any interesting physics in this limit A/Kp ^ ? If the cut-off is going to zero 
in laboratory units, are there any Feynman diagrams at all or, do we just read off all the 
physical scattering amplitudes from the vertices in the action, with no loop corrections, there 
being nothing left to run in the loops? In other words, is the full vertex function F the same 
as the (duly antisymmetrized) couplings F in the action? 

The answer depends on whether we are looking at = or 0, ie.e at processes in 
which the final pair lies or does not lie in the same plane as the initial pair. As for the latter 
it is indeed true that tree level amplitudes in the action would not be dressed by any loops 
due to the kinematics in c? > 1. For example at one loop, the loop momentum would be 
restricted in magnitude and angle to size A and A/Kp respectively. This, coupled with the 
fact that these diagrams have no singularities, would make them negligible and we would 
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have r = F. 

For forward scattering however, the iterated loops in the ZS channel, with only forward 
scattering couplings (F) appearing in all the vertices, would have no restriction on the loop 
angle, which is why the graphs survived in the large picture. However the magnitude of 
loop momentum k would still be bounded by A. Why would these graphs survive in the 
limit of vanishing cut-off? The answer is that the integrand has a S- function singularity 
at the Fermi surface (derivative of the Fermi function), if the external frequency transfer is 
zero, rendering the integral insensitive to A as long it is non-zero. 

Any attempt to introduce a nonforward amplitudes F{(f) ^ 0) into the iterated ZS loops 
would bring in a suppression factor, call it 1/N or A./Kf. Thus the forward and nonforward 
amplitudes do not mix. If we focus our attention on computing responses to soft probes 
((a; Q) « A), we can ignore the nonforward F's. The resultant theory is just Landau's 
Fermi liquid theory . The function F{z, (p — 0) is called the Landau parameter F{z). 

Landau's F[z) can be introduced in another equivalent way. Let us begin with our fixed 
point theory (hereafter in d = 3) 



then there are no self-energy corrections in the limit we are in. Imagine now a state with a 
macroscopic number of (quasi) particles added to the ground state so that all states up to 




(400) 



where u can contain nonforward amplitudes as well. We have seen that if 




(401) 



momentum k — r(n) are occupied in the direction VL. We can make such a state the ground 



state by modifying the hamiltonian to : 




(402) 
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The action now becomes 

Let us now calculate the energy e'{k Q) (associated with H') of a particle in a state labeled 
by (Qk). It is found from G~^, which we can calculate exactly in the large limit since we 
just need to evaluate the tadpole. This gives: 

e'ik n) = v*{k - rm - / — \ - 5/.*. (404) 

Doing the uj integral gives us 6{—k' + r{Q')) whereas the integral in 6fi* gives 6{—k). As a 
result, the energy e of a particle, now measured with respect to H (which differs by v*r from 
€') is : 

e{kn) =v*k + J j^F{n ■ n')r{n') (405) 

It is evident that the integral involving F represents the interaction between the (quasi) 
particle in question and the rest. Consider now a state with 6n{k Q) quasi-particles at the 
point {k Q), with 6n = —1 if it is a hole. In terms of r{Q), 

Sn{kQ) = e{k)e{r{Q) - k) - e{-k)e{k - r{Q)). (406) 

The energy of such a state, with reference to the ground state follows from Eqn.(|405|): 



E [6n] = v*k6n{k ^) + ^Y.Yl ^^(^ ^)F{Q ■ n')6n{k Q') + 0{6n^) (407) 

k a kd k'Vt' 

where V is the volume, kept finite so we can do a sum rather than integral over momenta, 
and F is proportional to F, which is how Landau introduced his F. 

Some texts (Mahan 1981) devote some time to why Landau went on to keep the quadratic 
term in 5n, and if he did, why he did not go to higher orders. The RG approach provides its 
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own version of the answer. Both terms (coming from ip{iuj — v*k)il) and ipipipip ) are marginal, 
whereas higher powers of ifj would bring in irrelevant operators. In particular the term with 
four powers of competes with a term with just two, since the latter has an extra k or 
multiplying it and these renormalize downwards under the RG transformation. 

Of course we did not have to wait for this analysis to understand why Landau did what 
he did. If we go from Eqn.(?n3 ), which gives the energy of an excitation, to the energy of 



all of the excitations we find: 



E 



V (27r) 



(40^ 



(27r) 

We may view the above expression is representing the elastic energy of a membrane or rubber 
band representing the Fermi surface with r(fi) as the deformation parameter. (This equation 
may be found in page 54 of Nozieres and Pines.) It is now clear that both terms are of the 
same order in the deformation. Haldane (in his lectures, yet to be published) emphasizes 
this aspect of Landau theory. He writes down the effective H as a quadratic function of 
some current densities obeying an algebra familiar in Conformal Field Theory. It is clear 
that we are discussing a solvable theory. The fact that only forwards scattering interactions 
{F{z., = 0)) enter this theory is what makes it solvable and also endows it with additional 
symmetries. A very concrete application of Haldane's approach maybe found in the work of 
Houghton and Marston (1992). 

More recently Castro Neto and Fradkin (1993) used a coherent state path integral to sum 
over the configurations of the Fermi surface. 

This concludes the link between the present formalism and Fermi liquid theory . Once 

117 



we have the concept of the Fermi hquid theory there is no need for the RG . However, for the 
sake of those RG minded readers who have followed all these arguments, but are not familiar 
with Fermi liquid theory , three sample problems will be discussed, not only to provide some 
instant gratification but also because each of them tells us something very instructive. 

A Landau Theory for the masses 

Our fixed point theory parametrized by v* and u* F has evolved from some bare theory 
with mass m, coupling U and so on. In the RG approach, one does not usually attempt 
to reconstruct the bare parameters in terms of the final fixed point quantities due to the 
unavoidable loss of information that accompanies mode elimination. There is however an 
ingenious argument due to Landau which does exactly that by relating m to m* and F. The 
reason not all information about bare quantities is irreversibly lost is due, as always, to a 
symmetry, Galilean invariance being the operative one here. Let t/ be a unitary operator 
that acts on a state > and gives it an infinitesimal boost with velocity 

5v — 5p/m. 

Under this active transformation , the energy of the eigenstate 1-0 > changes as follows: 

5E ^<U'il:\H\U'il:> - <'il:\H\il:> . (409) 

Now Galilean invariance is the statement that the boost affects only the kinetic energy of the 
particles and not the interaction energy. Thus the response of H to this transformation is 
the same as in free-field theory and given by : 

U^Hu = H + P- — + --- (410) 
m 

where P is the total momentum operator. Thus to first order 

< Uil}\H\Uil) >-< il}\H\il} >=< il}\P ■ > . (411) 
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Let \ip > he a state containing one extra particle at the Fermi surface Fermi surface . Since 
the ground state has zero momentum, this is a state of momentum Kp in the chosen direction. 
Let the boost be in the same direction. The energy change according to the right hand side 
of of the above equation is 

5E = Kp—. (412) 

m 

As for the left hand side, the active transformation has three effects which change the energy 
to first order in the boost: 

• The quasi-particle momentum goes up by 5p and by the very definition of effective 
mass or Fermi velocity, its energy goes up by Kp^- 

• The sea gets bodily shifted by 5f). This does not affect the sea kinetic energy since the 
total momentum of the sea was zero. 

• The interaction of the quasi-particle with the shifted sea changes its energy as per 
Eqn. ( ^05|) with r{9) = 6p cos 6, where the angle 6 is measured relative to the boost. 

Adding all the pieces and equating the result to what was given in the previous equation, 
we get the famous relation: 

K/jl = K,^^SpfF(z)z^ (413) 

m m* J Att^ 

where z = cos 6. In terms of the dimensionless function 

^ m* F 

and the expansion in terms of Legendre polynomials 

^{z) =Y.<^iPi{z) (415) 
I 

we obtain 

rrf 1 , ^ 

— = 1 + -$!. 416 
m 3 
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Zero Sound 



Zero sound refers to natural oscillations in a Fermi liquid resulting from interactions between 
the particles. We find it just as we would find the natural frequencies of an oscillator: by 
looking for poles in certain response functions. This means the density-density response 
function or compressibility in our problem. 

Let us imagine an external probe (f){Quj) which couples to the density p{Quj) of the 
fermions producing density fluctuations. The compressibility x is given by 



{2nY5^'\Q)x{Qoo) = -< p{Quj),p{-Q - uo) > 



(417) 



We are interested in the limit Q,uj « A. We will compute the correlation function by 
using the diagrammatic rules we have employed so far. The first few diagrams are given by 
Fig. 13. Because of the fact that the radial part of the 5-function {5{k — k')) that we pull out 
of our graphs differs from the traditional one {6{k — k')/Kp) used in Eqn.( p7| ), we find: 



[-KM 



dkidTtiduJi 



iuo — v*ki iuoi + iu — v*ki — v*Qziq 



+ 



dkidVtiduJi 

QziqdZiQ 1 



iuj - v*QziQ Atc 



I{uJiki 
+ 



/(t^ifci) 
dk2dQ2d^^2 

QziqdziQ 1 
iui — v*QziQ 47?^ J iuj — v*Qz2Q 47r 



i{to2k2)F{n^ ■ n2] 



Qz2QdZ2Q 



+ ■■■ (418) 
' -F(fii ■ n2) 4(419) 



where ZiQ is the cosine of the angle between Q and the direction of the loop momentum Ki. 

Now, we would like to study the simplest problem of this kind and therefore like to choose 
a constant F. Unfortunately, F, although not antisymmetric under the exchange of 1 and 2, 
still vanishes when the initial angles coincide. (See the nearest neighbor example wherein 
F ^ (1 — cos ^12).) We will compromise and for once introduce spin. We will assume only 
up and down particles scatter, with a constant Fq. If we now look at Figure 14, we see that 
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the first loop gets an extra factor of 2 due to the spin sum. The second also gets only a 2 
since spins at the vertex must be opposite. This restriction continues down the chain. Now 
we find that the series is geometric. The sum gives 

X(Q-) = (420) 

where we have chosen to look at real (rather than Matsubara) frequency uj since we are 
looking for real propagating excitations. 

Notice how we get an answer that is very sensitive to whether Q/uj or vice versa. 
In the former case Jq vanishes whereas in the latter case it equals — ^^2^, . This is the kind 
of sensitivity that plagues the four point function F also, as alluded to earlier. 

Clearly a pole occurs occurs in x when 

1 If zdz 



where 



Fq Att'^v* J s — z 

' ^ ^ (422) 



2 s-1 



Qv* 

In terms of the dimensionless $; introduced in Eqn.(^T^): 



(423) 



1 

$0 



2 s-1 



(424) 



For the solution to exist, we require s > 1, i.e., the velocity of propagation uj/Q must exceed 
the Fermi velocity v*. 

We will not discuss the extensive physics of this phenomenon. (For example, Mermin 
(1967) has showed that we will always have a zero-sound mode, for any reasonable F. )The 
main point was to show that the narrow cut-off theory has a lot of life in it. 
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C Static Compresssibility 

To find the equilibrium compressibility, we simply set a; = , Q — > in the preceding 
calculation. This means that 

Qz 1 

^ (425) 



uj — v*Qz V* 

The iterated integrals then simplify to the point that we no longer have to introduce spin or 

assume F is a constant. Going back to the spinless case, we find: 

X = ^^^^^ (426) 
^ l + $o 

Had we computed x free-field theory we would have found 

xo = (427) 

Thus 

X_ ^ m*/m _ l + $i/3 
Xo l + $o l + *o ■ 

Now, no one will dispute that this is indeed the ratio of < ipipipip > correlation functions in 
the fixed point theory to the free-field theory. But in Landau theory one equates this to the 
ratio of compressibihties. This is not so obviously correct and I thank N.Read for forcing 
me to clarify this point. The problem is this. Let us begin with the full path integral over 
the bare fields prior to the RG transformation with a coupling of the bare charge density to 
some external field A: 

Z{A) = J diPocbPoe'^^^^ (429) 

1 
2!2! 



S{A) = I ^o(^a; - E{K))^o + 7^ I ^o^o^o^oU + A I ^o^o- (430) 



Now, it is certainly true that 

^2 \nZ 



dA^ 
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< '0o'0o'0o''/'o >= X (431) 



Suppose we now perform the RG transformation. In the process we rescale the field: 
TpoZ'^^"^ — ip. (I apologize for using Z to denote two different things. Hereafter we will 
only see the above definition, as the wave- function renormalization factor.) This means 

Now the partition function is preserved by mode elimination, and we can take its second 
logarithmic derivative with respect to A after the RG transformation to find But this will 
give < ipipipip > whereas we computed the operator without the Z's in what we called x 
of the fixed point theory. 

So this is the mystery. The resolution lies in the fact that besides the rescaling, a term of 
the form Aipip is generated (in addition to what was already there) as we carry out the RG 
transformation, and this precisely cancels the effect just discussed. This is an example of a 
Ward identity based on charge conservation. (See Abrikosov et al (1963) for a discussion.) 
Here is a glimpse of how it works with no numerical factors. Let us look at two graphs that 
cancel. (All loop momenta he in the eliminated region and correspond to fast modes.) Take 
the sunrise diagram (Figure 15a) whose iu; derivative at a; = contributes to Z, the field 
rescaling factor: 

icu iuj(l - — ) = iujZ-\ (432) 
iou 

Imagine routing the external momentum through the upper line. Taking the icu derivative 
clearly squares that propagator. This is shown in Figure 15b, with the cross denoting the 
place where the second propagator joins the first. 

Consider now the other phenomenon: generation of new terms. In the mode elimination 
scheme the coupling between A and ipip can take place via the fast modes as shown in Figure 
15c. Notice that this diagram coincides with that in Figure 15b in the limit when the probe 
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brings in zero momentum and frequency. Consequently the field rescaling effects ( due to the 
the self-energy diagram) precisely cancel induced terms effects (due to the vertex correction 
diagram). Although we took just a pair of diagrams, the result is exact. It reflects the fact 
that even though the quasi-particle can break up into many particles, (so that its chance of 
being a single particle is reduced), the field can couple to the fragments now, and the total 
charge of the fragments ( which is all the field couples to in the limit of zero frequency and 
wavelength) is that of the quasi- particle. (The careful reader will ask: what about the c- 
number term of the form in the action that comes from integrating fast modes? These 
contributions from the fast modes drop out as the external momentum and frequency vanish, 
which is the limit we are interested in.) 

Thus a lot of Landau theory acquires its power due to the fact that not only are many 
quantities (like < ipipipip >) computable in the fixed point theory, they directly correspond, 
with no intervening, unknown prefactors, to physical observables (like compressibility) due 
to Ward identities. 

D Notes for the experts 

Here are some notes for readers who are familiar with the details of one or other ideas invoked 
earlier. 

• In the diagrammatic treatment of Fermi liquid theory , one organizes the graphs as 
follows (in the notation of Abrikosov et al.(1963) First, one looks at the theory in the 
limit where the external transfers Q,u; — > 0. The full four-point function is called F'^ 
and corresponds to the limit lu/Q ^ 0. It is given as a sum of diagrams where the 
bare vertex is called F'^, which corresponds to the limit Q/co — > 0, and is irreducible 
with respect to a pair of particle hole lines which are singular at the Fermi surface . It 
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is assumed that the bare vertex is analytic in its arguments, and the trouble makers, 
the particle-hole lines which produce all the singularities in the small {Q,uj) limit, are 
explicitly displayed. 

In the RG approach, the bare vertex u contains all the safe modes which include 
particle-hole lines with at least one of them outside the cut-off. The only lines shown 
explicitly in the cut-off theory are particle -hole lines both within the cut-off, these being 
the modes yet to be integrated. 

Couplings corresponding to non-forward scattering, called F{z, 7^ 0) in this paper, 
are very important for the study of lifetime effects and transport properties. A nice 
discussion of this may be found, for example in Mahan's book (1981). Note in particular 
the discussion of the work of Dy and Pethick in page 947. These authors asked how the 
forward scattering Landau parameters F(z) may be extended to non- forward scattering 
and argued for a certain dependence based on symmetry under exchange. The 



additional factor of cos0 they came up with is exactly what we find in Eqn. (|372| ) of 
this paper. 

Even though F is a marginal coupling, there are no anomalous dimensions for the 
operators, in contrast to the fixed line in d = 1 along which the fermion field and 
other operators have continuously varying dimensions. It is worth finding out if the 
interaction is a redundant operator in the sense of RG . 

Landau theory appears very much like a classical self-consistent theory. We understand 
this as coming from the large N saddle point, which, like any saddle point, represents 
a form of classical limit. 

Some readers familiar with 1/A^ expansions may ask how the F's manage to change 
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physical quantities like x by factors of order unity, when one always thinks of interac- 
tions as producing changes of order (In other words, why are the dimensionless 
numbers $, appearing in, say, Eqn.( ^2^ ), producing corrections of order unity? ) The 
central feature of expansions is that loops of a certain kind are order unity since 
the loop sum pays for the extra factors of the interaction. If we compute a scattering 
amplitude, which exists only due to interactions, we get a term of order since that 
is the strength of any one coupling. This is a sum of a tree level or bare term and 
iterated loops, all of the same order. In other words the quantum corrections in the 
iterated loops are of the same order as the tree level term. But if we compute some- 
thing like the density- density correlation function (not often done in field theory), this 
is given by the polarization bubble which is nonzero in free-field theory and hence of 
order unity in the series . Loop corrections to it (ZS graphs) are of the same order 
as well. 



IX NON-CIRCULAR FERMI SURFACES: GENERIC 



We now discuss a Fermi surface which has no special symmetries other than time-reversal 



extension to = 3. A surface that meets these conditions is an ellipse and is depicted in 
Figure 16. 

The first step is to set up the RG transformation for the noninteracting problem. Since 
\K\ is no longer a measure of energy, we must draw contours of constant energy e (measured 
from the Fermi surface) and retain a band of width A in either side of the Fermi surface. 
Thus our starting point is the action: 




Once again we focus on = 2, discussing in passing the 




(433) 
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where 9 parametrizes the Fermi surface, and J is the Jacobian for (K^Ky) — > {9e). We will 
expand 

J{e9) = J (9) + eJi(9) + ■■■ (434) 

around the Fermi surface and keep just the first term; the rest will prove irrelevant. 

The RG transformation is exactly as before, with e in place of k. It is clear that higher 
order terms in J renormalize to zero with respect to this RG transformation. The interaction 
term is 



5^4 = ^ / ^(4)^(3)^(2)^(1)^(4, 3, 2, 1) (435) 

where 

e-l^^l/^. (436) 



A Tree Level analysis. 



J{9i)d9id£iduji 



The analysis proceeds exactly as in the rotationally invariant case. As A, (the cut-off in 
energy now) is reduced to zero (in fixed laboratory units) by mode elimination, we find once 
again that the set of initial momenta 1 and 2, must coincide with the final set 3 and 4. To see 
this, we must simply replace figures with intersecting circles with say, intersecting ellipses 
(Fig. 16) or whatever may be the shape of the Fermi surface. In smooth cut-off only such 
couplings will be spared exponentially downward renormalization. 

So once again the marginal couphngs at tree level obey ^3 equal to 9i or ^2 unless 9i — —^2 
in which case 6*3 = —6*4. (The figure shows the second possibility). The tree level amplitudes 
are labeled as before: 

u{929i929i) = -u{9i92929i) = F{92 9i) (437) 
u{-9s9s - 9,9i) = V{93 9,) (438) 
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but are no longer functions of the differences of their arguments. 
B Tadpole graph. 

We now see a new feature with the tadpole. Let us understand this without any reference 
to the RG . Suppose we begin with the action Eqn. ([435|) , and go to one loop. The tadpole 



graph, Fig.5, makes the following contribution to the self-energy 

neue) = AjF{ed')^-^^^so{e). (439) 

Thus 

= {iu}-e- eo{e)) (440) 

and the Fermi surface has moved to e = —eo{6). This generally involves a change in shape. 

Now even in the rotationally invariant problem the Fermi surface moves; turning on 
interactions changes the Fermi surface radius from Kp = ^^2m^ to Kp{fi) such that 

^ + S(«) = (441) 

In that case we can add a counter term Sfi* = eo{6) = Eq to the action to restore the old 
Fermi surface radius. Recall that the action (schematic) 

S = J i^iiuj - + ^0 / V^V' + ^ / uiJ^iji; (442) 

was invariant under the RG at one loop and order u. 

Now, it was pointed out that adding a counter term to maintain the radius of the Fermi 
surface was not fine-tuning since it corresponded to maintaining a fixed density, which is 
experimentally viable. It was also pointed out that if we did not add the counter term, the 
system would not acquire a gap, as in 0^, theory, but would simply move to the new radius 



defined in Eqn.(441). 
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In the non-rotationally invariant problem clinging to the old Fermi surface is certainly a 
case of fine tuning. There is no experimental way (such as sealing off the system) to preserve 
the detailed shape of the Fermi surface . We could add a constant 6fi to keep the density 
(i.e., volume enclosed by the Fermi surface ) constant, but we will not; we will let the system 
find the true Fermi surface for the given fi and u. We can use the RG to determine the final 
Fermi surface as follows. We rewrite the initial action by adding and subtracting a presently 
unknown term £o(^)-' 



^{tuj -e)^+ I eo{e)^^ + / ui/ji/jiPij (443) 



2!2! 

- / nil 

2!2! 

and demand that this be a fixed point. To one loop and order u this means: 

rA r2n roo dl^dO' duj' J{e')F{e 6') 



em 

which leads to 



em - 



A/sJo J-oo {2TxY{iuj — e') 



(444) 



^m = A / . (445) 



To obtain this we have to do the uj integral in Eqn.( ^44D . Note that in that equation 



everything in the integral is evaluated at order m*^ due to an explicit F in it. This means 
that the propagator in the integral has an angle-independent e' . 

Since we now have a fixed point it must be true that we have found the correct Fermi 
surface . From the knowledge of e* we can reconstruct the new Fermi surface . In principle 
one could go order by order in this " renormalized perturbation theory" . However in the large 
limit which appears here also, the one loop answer gives the full self-energy correction or 
change in Fermi surface . 
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C One loop at order 

If we do the mode elimination as in the rotationally invariant case, we will find once again 
that The ZS and ZS' graphs once again do not contribute to the /^-function for F the same 
reason as in the rotationally invariant case: either there is not enough momentum transfer 
to knock the internal line at —A to A or, there is a kinematical suppression factor dK/Kp. 

If we do mode elimination as in the rotationally invaraint case, we find once again that 
the ZS and ZS' graphs do not contribute to the fiow of F for the same reason: either 
there is not enough momentum transfer to knock an internal line at —A to A, or there is 
a kinematical suppression factor dK/Kp- The fiow in the BCS channel is unaffected by 
non-rotational invariance. As long as E{K) = E{—K) the BCS diagram, given by the third 
term in Eqn.( pT5D , will make a contribution and we will obtain: 

dv{ei;es) 1 f-^^de. 



/ vie^;e)vie;e,)Jie) (446) 

Jo 2tc 



dt 87r2 Jo 27r 

though we can no longer use rotational invariance to decouple this equation using angular 
momentum eigenfunctions. It is however possible to do a double Fourier expansion. This 
deserves further study analytically and numerically. 

Weinberg (1993) has recently derived such a flow equation for superconductors whose 
Fermi surfaces obey just time-reversal invariance using the notion of effective actions from 
quantum fleld theory. Besides the flow, his paper has a careful derivation and analysis of 
the effective action for superconductors. All the important properties of the superconductor 
may be derived from its effective action and the notion of broken gauge invariance (Weinberg 
1986s). 

If we ignore the BCS interaction, we expect the Fermi liquid, 1/A^ etc., to work as 
before except for lack of rotational invariance. The F function will now be a function of two 
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variables. We expect to find zero sound. We do not expect any simple relation between m 
and m* due to lack of Galilean invariance. 

Finally if we consider a Fermi surface without time-reversal symmetry, we can get rid of 
the BCS amplitudes even at T = 0. (However, in drawing the analog of Fig. 8, we must draw 
the time-reversal inverted version of the second Fermi surface (displaced by P ) since the 

— * — * 

previous construction assumed that if K is an allowed vector, so is —K. 

X NON-CIRCULAR FERMI SURFACES: NESTED 

We are finally going to discuss spinless fermions on a square lattice at half- filling. A specific 
model for this problem is the one with nearest-neighbor interaction; 

H = Ho + Hi (447) 

= E ^l^Kmf) + h.c. + Uo E (^^0>W-^)(^njW)-^), (448) 

<jj'> <jj'> 

where j labels sites on a square lattice and the subscript < jj' > on the sums means j' is 

restricted to be the nearest neighbor of j in the direction of increasing coordinates. (Thus 

if j is the origin (0, 0) , j' is restricted to be (1, 0) or (0, 1).) The chemical potential (found 

by opening up the brackets in the interaction term) is 

= 2Uo (449) 

where the the factor 2 comes from the number of nearest neighbors. 

At C/q = 0, the half- filled system will once again be a perfect conductor as can be seen 
by going to momentum states. Likewise at ^/q = oo there will be a charge density wave with 
more charge on one sublattice, say the one whose x and y coordinates (which are integers in 
lattice units) add up to an even number. 
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Once again we will focus not so much on the fate of this one model but rather on a class 
of models described by the same free-field fixed point and its perturbations. Let us therefore 
find the fixed point describing the noninteracting problem. 

Let us take as the free-fermion dispersion relation: 

E ^ - cos - r cos Ky. (450) 

which corresponds to the problem with unequal hopping in the two directions. The reason 
for the choice r 7^ 1 will follow shortly. Notice that we expect to see a CDW state at large 
repulsion in the nearest neighbor model even if r 7^ 1 (since in this limit the hopping term is 
not the deciding factor, the interaction term is). Notice also that for any r, E still has the 
symmetry: 

E{K + Q'n) = -E{K) Qm = (tt , tt) (451) 

which is all we will need. The chemical potential that gives rise to half-filling aX U = is 
= which means that E — e and that the latter also changes sign when we add the nesting 
vector: 

e{K + Qn) = -e{K) (452) 

The Fermi surface for r > 1 is sketched in Figure 17. The dark line shows the Fermi surface 
which now has two branches a = ±1. Each point on the Fermi surface goes to another 
point on the Fermi surface upon adding Qat- This means that if we shift the figure by 
Qjv, the shifted one (in the repeated zone scheme) will fit perfectly with the original like 
something out of Escher's drawings. This means that if the momentum transfer is Qat, the 
analog of Fig. 10 will show the complete overlap of the two displaced surfaces rather than 
their intersection. Also shown contours of energy e — ±A. This where where the modes 
to be eliminated under infinitesimal renormalization lie. Points at ±A are knocked to ^A 
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upon transfer of Qn. This in turn ensures that if one loop momentum hes in a shell being 
eliminated, the other (in the graph with momentum transfer also does. This will lead 
to a flow of the coupling in which final momenta differ from the initial ones by Q^. 

For nesting to take place, we need half- filling and the symmetry in Eqn.(HS^. The latter 



comes if we assume that hopping is always from one sub-lattice to the other in a bipar- 
tite lattice. Since hopping is usually just nearest neighbor to an excellent approximation, 
studying the effects of nesting once again do not constitute fine-tuning. 

Let us write down the action for the noninteracting problem. We will use as the final 
coordinates e and 6 = together with a discrete index a = ±1 which tells us which of the 
two branches we are on. Thus 

- ? £ I £ I L - («^) 

where 

J{ee) = ^ . (454) 

Jr"^ — {e + cos 6y 

It is clear why we introduced r 7^ 1: if it equaled unity, J would be plagued with (van Hove) 
singularities on the Fermi surface and the expansion in e would be impossible. Since nesting, 
and not van Hove singularities is what we are interested in here, we will study r > 1. In this 
case the Fermi surface value of the Jacobian is 

yr^ — cos^ o 

Henceforth the subscript on J will be dropped. 

Mode elimination of the action in Eqn. ([45 3|) and the rescaling of fields and e go as before 
to render the action the fixed point. 

As for the quartic term, there are now three sets of couplings that are marginal at tree 
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level. Besides F and V, we also have 

u [62 + TT, -a2; 61 + TT, -CKi ; 62a2; 9iai\ = -u [61 + n, -ai, O2 + tt, -a2, O^oci, OiOii] 

= W [62 a2 ; ^1 ai] (456) 

which corresponds to processes wherein the momentum transfer between 1 and 3 or 2 and 3 
equals Qjv- In this case, because of the nesting property of the Fermi surface , we are assured 
that the fourth momentum will lie on the Fermi surface if the first three do, and this, you 
recall, is the condition for the coupling to survive the tree level RG transformation. Note 
also that W describes umklapp scattering if both particles start out on the same branch and 
hop to the opposite one. 

As for the one-loop /3-function, the F's do not flow and y s have the usual flow from the 
BCS diagram 

1 foo du; f (PK u[-KsKs - K K]u[-K K - K^] 
du [-K, K, - K, K,] = -- - J^^^^^ ^ (iu; - e{K)){-iu; - e{-K)) 

(457) 

where "shell" refers to the shell being eliminated around both branches. The two contribute 
equally to give, in terms of V: 

Since we do not have rotational invariance, we cannot separate this using the angular 
momentum variables. We can however use the double Fourier transform and reduce it to 
discrete coefficients which will be coupled in their evolution. 

Let us now look at the flow of u [K' 2K' iK2K]\ where K'i = + Qjv^due to the ZS 
diagram: 

^ ^ ^ ^ J -00 2ti J shell {2tiY {iuj - s{K)){iuj - e{K')) ^ ' 
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where "shell" means both e{K) and £{K') lie in the thin shells of width dA near ±A. Q Due 
to the nesting property two interesting things happen leading to a flow: 



• If e{K) lies in the shell, so does e{K') = —e{K). 

• the u integral never vanishes since the poles always lie on opposite half-planes. 
Doing the u integral, we get, in terms of W 

= - / E W la,A a' »'| W lea O.a.l ^ (460) 

at J a l^^j 

where (a' 9') refer to K' . 

To get a feel for this problem, let us evaluate the nearest neighbor interaction on the 
present Fermi surface to obtain: 



. 2/^1 -^2x , 1/, cos6'icos6'2 aia2 rz — n; — 

sm ( — - — ) + -(1 — ~ cos^t^iyH - cos^6^2 



(461) 

and 



V [^3^3 Oiai] = Uo 



sin 9i sin 63 H — y'r^ — cos^ Oi\Jr'^ — cos^ 6^ 



(462) 



It is readily verified that W always has the same sign for all values of its arguments; 
opposite to that of Uq. Specializing to the repulsive case, it is clear from the above equation 
that > 0. Thus we have proven the instability for this initial condition. (This is a weak 
coupling argument. To describe the nearest neighbor problem in the small A theory it is 
not enough to simply restrict the full interaction to within the cut-off, we must take into 
account induced terms and renormalization due to elimination of modes. These are however 



-•^^We have defined W with i^s = and for this choice only the ZS diagram contributes to flow. Had we 
reversed the role of and ^4 in the definition of W , an extra minus sign would have entered it deflnition, 
and the ZS' (rather than ZS ) diagram would have contributed to its flow. The latter would also have had 
a minus sign relative to the ZS diagram. 
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higher order effects.) As for some other interaction, if it has any overlap with this direction 
we have found, it will be unstable. Note that unlike in the rotationally invariant problem, 
we have not explicitly displayed an infinite number of unstable directions (one for each /). 

Unlike in the rotationally invariant case, the present problem has many open questions. 
Here is list some of the more important ones. 



The study of flow equations Eqns.( ^5q - 460 ) is a very important follow-through to 



the work presented here. Calculations are being performed with G. Murthy in which 
the couplings from the nearest neighbor model (Eqn.( ^61 - 462 ) ) are taken as initial 



conditions. In studying the flows, it was important to remember that there are some 
couplings which flow for more than one reason. An example is when the initial momenta 
are K and —K and the final ones are K' and —K' where K' = K + Q this coupling 
is equal to a and a W. The general idea is to run the flow till the cut-off is small and 
then solve the theory by summing over diagrams that survive in the limit of vanishing 
A (our 1/A^.) 

Although the formalism was discussed in connection with Fermi liquids with a 
rotationally invariant Fermi surface, it applies to all problems discussed here: one 
can always reduce the cut-off (keeping track of the evolving couplings) and then use 
the smallness of A to do a sum over diagrams that lead in the expansion. For 
rotationally noninvariant problems, this can however be complicated by the constant 
motion of the Fermi surface as we renormalize, even if we keep its volume constant by 
modifying the chemical potential. In other words, as the modes are eliminated, the 
new Fermi surface and new one-particle energies must be used in defining contours of 
constant energy and choosing modes for the next round of elimination. At the one 
loop level considered here, this was a nonissue. 
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• It is important to consider the problem just below half-filling: here we expect that 
there will be an initial growth of interactions as we lower the cut- off, which will then 
freeze once the cut-off is comparable to the deviation of the Fermi surface from nesting. 
One must see if in this case an attractive interaction is generated in the meantime in 
the BCS channel. If so, this couphng will continue to grow, since it does not rely on 
nesting to do so. This will be yet another test of the notion that attraction can lie 
hidden in models that started out repulsive. 

• Another topic worthy of further study is the coupling u{—A, A, — A, A) where A — 
^Qn — (7''/2,7r/2) and lies on the Fermi surface for any r see Fig. 17. (The present 
remarks apply equally well to the point A' = {71/2 — n/2).) This coupling is a 
V & W and an F . It receives flow contribution from the ZS' and BCS diagrams 
and is forward scattering amplitude which will control particle self-energies. Is the 
distinguished nature of this point (from the point of view of the RG) related to why 
so many investigators ( Sachdev 1989, Trugman 1988, Elser et al 1990, Boninsegni and 
Manousakis 1990) find holes occurring at the point A when the half-filled systems are 
doped? Once could pick the initial value for the flow the coupling generated by the 
nearest neighbor interaction, follow the flow and look at the effective theory at very 
low energies to see if an answer comes out. 

• Consider the problem on a square lattice without nesting, say because of same sub- 
lattice hopping. At very large nearest neighbor repulsion, we can see that a CDW will 
result, with more particles in one or the other sub-lattice. However, there will be no 
instability at inflnitesimal repulsion. One expects from continuity that the transition 
will take place at small coupling for small nesting violations. One should then see this 



137 



phase transition at weak coupling from the RG. 

• Notice that our analysis depended on the nested Fermi surface. While the Fermi 
surface was nested in the absence of any interaction, don't the interactions cause it to 
move? Won't the shape change from perfect nesting even if we change the chemical 
potential to sit at half-filling? What happens to the CDW instability then? First of 
all this question does not affect our one loop calculation which uses the zeroth order 
propagators with their zeroth order formulas for e) . Whether or not the nested surface 
will stay nested at higher orders in the interaction is an open question which could 
control the higher terms in the /9-function and decide the ultimate destination of the 
flow. It will not however change the fact that the free-field fixed point is unstable since 
that has been established close to the fixed point by our one loop calculation. However 
if the CDW instabihty is to really take place, the fiow must keep going till we hit a 
fixed point with a gap. Now, it is clear in coordinate space that at strong coupling in 
a bipartite lattice, there will be a CDW. For this to come out of the RG, nesting must 
be preserved as we renormalize. This is however a conjecture and has not been proven. 

Another interesting question is the following: if we set the hopping coefficient r — 0, 
we seem to decouple the chains. Will we then get a Luttinger hquid? No. This is because 
the interaction terms couple the various decoupled chains. The exact cancellation that took 
place between the BCS and CDW instabihties in a one chain model with just one couphng 
will not repeat itself anymore. 

XI NON-FERMI LIQUIDS IN d =2 

So far we have seen two means by which the Fermi hquid could be destroyed: BCS and 
CDW instabilities. In both cases, the flow came about because individual Feynman diagrams 
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had logarithmic singularities. Thus the perturbation series had zero radius of convergence. 
(There were essential singularities of the form e~^/" in, say, the CDW order parameter. ) 
If perturbation theory can tell us about the instabihties, why follow the RG route? The 
answer, from the d = 1 example is that even if individual diagrams diverged, it was possible 
for the /3-function to vanish, producing novel scale invariant behavior. We are looking for 
such a state in o? 2. I see no evidence for it if 

• the Fermi surface is spherical, 

• the coupling is weak, 

• the input interaction is short ranged. 

• we work in infinite volume from the start. 

In particular, I have examined, together with A.Ruckenstein and H.Schulz the channel ana- 
lyzed by Anderson (1990): the incoming particles of opposite spin at the same momentum 
on or near the Fermi surface . We found that there was no flow in this coupling as the cut- off 
went to zero. Indeed there was no singularity in the diagram when the external (euclidean) 
frequencies vanished. Setting them to nonzero values did not help. Of course setting them 
equal to some real frequency did cause singularities, but these correspond to propagating 
modes and not instabilities of the ground state. (Recall that the usual instabilities were 
seen at zero external frequency.) But it must be pointed out that this is not in variance 
with Anderson's arguments which rely very much on doing things in flnite volume and then 
carefully taking the infinite volume limit. Unfortunately the RG in a finite volume is not an 
easy prospect and we were hoping that if (despite Anderson's cautionary note about going 
to infinite volume too quickly) the effect showed up in our calculation, it would be additional 
corroborations to Anderson's argument, but within the standard infinite- volume machinery. 

139 



Engelbrecht and Randeria (1991) who studied this problem in the low density expansion did 
not find any instability. The same remarks apply to their calculation also. 

We must therefore relax some of the above conditions. As mentioned above, Anderson 
dropped the infinite volume condition. Another possibility (Ruckenstein et al 1989) is that 
at strong coupling a new possibility, the marginal Fermi liquid liquid , which has impressive 
phenomenological success, arises. The present weak-coupling analysis has nothing to say 
about it and surely cannot exclude it. For example, it is possible for the bare coupling that 
enters the action of the narrow cut-off theory (or in the Abrikosov et al (1963) treatment) 
to be singular. Earlier it was stated that this would not happen because the bare couplings 
were obtained from the input parameters by integrating out safe modes. But this only 
assures us that the individual diagrams that add up to give the bare coupling are finite. It is 
certainly possible for the infinite sum to diverge beyond some maximum coupling. This does 
not contradict Landau's analysis since the basic assumption that the physics in question is 
a continuation of the noninteracting problem is invalid. The stability of the marginal fermi 
liquid has studied by Zimanyi and Bedell (1991). 

Another possibility is that even before mode integration, the input coupling is singular. 

Stamp(1992, 1993) has taken a pragmatic approach and considered the effect of singular 
interactions, setting aside the question of their origin. 

I decided to look at the Coulomb interaction in this fight. We cannot simply say that 
it gets screened; this is a picture that makes sense when a subset of diagrams in standard 
perturbation theory are resummed in a certain way to produce the screened propagator for 
the Coulomb potential. The book- keeping is different in RG : the coupling that goes into 
the action has not been screened by particle-hole pairs at the Fermi surface . It is easy see 
(within the sharp cut-off scheme) that at any stage Q << A is unscreened. The correct 
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procedure is to follow the evolution of the bare coupling as the modes are eliminated and 
see where it ends up when no more integration is left over. It is shown in the Appendix that 
the final potential is screened. (This is a smooth cut-off version of a sharp cut-off calculation 
devised with G.Murthy.) This analysis however assumes that the fermion propagator has 
the standard Fermi liquid theory form. Ideally one should let the fermion propagator also 
evolve as modes are eliminated and see if we still end up with a screened interaction. This 
has not been done. 

Now there are concrete examples of non-fermi liquid behavior if we are willing to consider 
impurity problems. Consider for example the example provided by Affleck and Ludwig 
(1991a, 1991b, 1992) from the Kondo problem. More recently Perakis et al (1993) have 
given another example from the Kondo problem which shows non-fermi liquid behavior for 
a range of parameters. 

Although the search for non-fermi liquid did not yield anything at weak coupling, it is a 
worthwhile goal since the copper-oxides seem to call for something different. As pointed out 
by Anderson, these may not be connected to the fermi liquid fixed point. Rather than reach 
the novel fixed point from the fermi liquid, one could attempt writing down different fixed 
points. They may require additional fields besides fermions, e.g., gauge boson (Nagaosa and 
Lee (1990), loffe and Kotliar (1990), Polchinski (1993)) : in a strong couphng field theory 
the low energy physics may bear no simple relation to the microscopic theory. Another 
route is to start with one or more one- dimensional systems, which can have Luttinger liquid 
behavior, and couple them perturbatively (Wen 1990, Schulz 1991). 
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XII SUMMARY AND OUTLOOK 

The main aim of this paper was to find a way to apply the RG methods to interacting 
nonrelativistic fermions, in particular to understand the various instabilities for weak per- 
turbations. In particular we wanted to see if the system remained gapless. Since the RG 
was so successful in dealing with critical phenomena, it was decided to follow a path where 
one relied heavily on analogy to this prior application. Since gapless systems corresponded 
to critical systems the idea was to use the language of fixed points and their perturbations. 

We started with a brief historical review starting from the original formulation of the RG 
for use in field theory, and ending with the modern approach pioneered mainly by Kadanoff 
and Wilson. It was pointed out that while the RG has always expressed the invariance of 
the theory under a change in cut-off followed by a suitable change in the parameters, the 
emphasis has shifted from viewing the cut-off as an artifact to be sent to infinity (where it 
belongs in a continuum theory like quantum electrodynamics) to viewing it as a dividing line 
between interesting and uninteresting degrees of freedom (the slow and fast modes) even in 
a problem where the cut-off is finite to begin with. 

In Section II we discussed the charged (complex) scalar field in four dimensions. It 
was shown how if one wanted to study physics at long distances, one could work with just 
the "slow" modes in the functional integral and how the unwanted fast modes were to be 
eliminated and the new couplings deduced. The gaussian fixed point was studied in detail 
and the fiows to one loop were deduced. It was shown how, upon reducing the cut-off to 
very small values, i.e., reducing the phase space to a tiny ball, the coupling functions reduce 
to a handful of couphng constants, which were just the first few terms in the Taylor series 
of the coupling functions about the origin. A comparison between the modern way to find 
the fiow equations {(3- functions) using mode elimination and field theory methods (trying 
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to get rid of cut-off dependence by proper choice of bare parameters) was provided. It was 
shown how although the two had very different book-keeping schemes, they gave the same 
answer for relevant and marginal couplings. 

In Section III the functional integral method for fermions was introduced. First a few 
simple problems in the thermodynamics and dynamics of fermionic oscillators were solved 
using operator methods. Then the Grassmann integral formulation was introduced and the 
same results were regained. The rules for calculating correlation functions (Wick's theorem 
for fermions) were derived and used in the calculations. 

The stage was now set for deahng with the nonrelativistic fermions. The strategy would 
be the following. We would first start with noninteracting fermions and write a hamiltonian 
that faithfully described the physics near the Fermi surface, i.e., within a cut-off A on either 
side of it. The logic was that the questions we were interested in, were decided by these 
modes at least at weak coupling. The corresponding functional integral would then be 
written. A mode elimination process that reduced the cut-off and rescaled momenta and 
fields and left the action invariant would be found. Given the RG and its fixed point, it would 
then be possible to classify the perturbations as relevant irrelevant or marginal. While this 
closely paralleled the scalar field problem one major difference was anticipated: since we are 
renormalizing towards the fermi surface and not the origin, the remaining phase space will 
be infinitesimal perpendicular the Fermi surface, but of fixed size in the tangential direction. 
Thus the fixed point would be characterized by coupling functions that depended nontrivially 
on the angles used to parametrize the Fermi surface. The case d — 1 was clearly exceptional 
since the Fermi surface consisted of just two points. This made the problem very similar to 
continuum field theories wherein the phase space for bosons and fermions is a ball centered 
at the origin. 
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The study was limited to spinless fermions. At this point the reader can surely see that the 
inclusion of spin really is straightforward. Since all the interesting flows were due to special 
properties of the shape of the Fermi surface, and since spinless fermions could display these 
shapes, they were the clear pedagogical choice. To keep the discussion concrete, a nearest 
neighbor spinless fermion problem at half-filling was frequently invoked. The question was 
whether it developed a CDW gap at arbitrarily small repulsion. 

The RG scheme worked remarkably well in the d — 1 warmup, described in Section IV. 
We took two slices (of width 2A) near each Fermi point (L/R) and found a mode elimination 
scheme that left the action of the free theory invariant. We then turned on a quartic inter- 
action. At tree level it was found that only the frequency and momentum [k — K — Kp) 
independent part of this coupling was marginal, the rest were irrelevant, as were all couplings 
with six or more fields. The marginal coupling did depend on the discrete internal index 
L/R but due to the Pauli principle, reduced to just one independent number. We had to 
go to one loop to resolve the fate of this marginal coupling. The /9-function route avoided 
the CDW and BCS mean-field instabihties by playing them against each other and giving 
the correct answer: a scale-invariant system, the Luttinger liquid. (This correctly describes 
the nearest neighbor model, which can be solved exactly.) We also discussed the quadratic 
terms that were induced by the quartic term upon mode elimination. These meant the Fermi 
surface was moving to take into account the interactions. There were two options. We could 
add counter term, determined order by order, to keep Kp fixed, or we could let the surface 
move to the new Kp. The latter would be found by doing renormalized perturbation theory: 
we split the original // into two parts, one which gives the correct Kp m. the propagators 
and a counter-term, which keeps it there. The division would be revealed to us order by 
order. It is a remarkable fact that Fermi systems, unlike Bose systems can remain critical 
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for a range of ji; they do this by using the fact that their fixed point is characterized by a 
whole surface, which can wiggle around as interactions are turned on. This idea applies with 
minor modifications in higher d, and we will not speak of it further. 

The scheme was then extended in Sections V and VI to rotationally invariant Fermi 
surfaces m. d — 2 and d — In d = 2 an annulus of width 2 A and radius Kp was used. 
The RG transformation was essentially the same as in o? = 1 except for the fact that we 
had an integral over the "internal" variable 9 rather than a sum over the Fermi points L 
and R. At tree level only two couplings, F{6) and V{6), which corresponded to forward 
scattering and Cooper-pair interactions survived. They had no dependence on external k 
or uj. Going to one loop we saw that there was no flow of F. This was because of the 
kinematics ol d > 1. There was however flow in V. We decoupled the flow into an inflnite 
number of equations, one for each angular momentum /. The flow was marginally irrelevant 
for repulsion (corresponding to the flndings of Morel- Anderson) and marginally relevant for 
attraction, corresponding to the BCS instability. We then formally set V = and identifled 
F with Landau's F function. Thus RG had led us automatically to Fermi hquid theory. The 
situation in c? = 3 was essentially the same, except for one difference that was signiflcant: 
while va. d — 2 only forward scattering amplitudes survived, in d = 3 nonforward ones did 
also. However the latter did not have any effect on low energy, low-momentum transfer 
physics. 

Having identifled the Fermi liquid as the flxed point of our RG, we next asked why it 
is solvable, i.e., why, despite the quartic interactions, one is able to calculate a lot response 
functions. It was pointed out that as we eliminate modes, a 1/N expansion emerges, with 
N — Kp/K. Landau theory is the N — oo limit. (By using certain collective coordinates, 
this limit can be made to correspond to a saddle point which gives the exact answer.) This 
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is the first example where really is large in the original problem. However a two stage 
assault was needed to ensure this: first use RG to bring A./Kf to a small value and then 
use 1/N. Of course the 1/N is not a new way to do Landau theory, but just a new way 
to understand it: it has always been known that bubble graphs (ZS loops) dominate the 
low momentum response. However, after the application of the RG there was nothing left 
but ZS bubble graphs due to the kinematics. The bubble graphs were ones in which the 
loop angle runs over all values. However these need not have amounted to much since the 
width of k integration was going to zero. But the graphs did survive since the integrand had 
a (5-function singularity on the Fermi surface and the integral was oblivious to the cut-off. 
(This is why the graph survived but did not contribute to the /9-function, which probes the 
sensitivity to cut-off.) 

An RG version of the Kohn-Luttinger instability was given. The details were relegated 
to the Appendix. 

We briefly looked at three effects in Fermi liquid theory since each taught us something 
and also gave the readers familiar with RG but not Fermi liquid theory, some instant grati- 
fication for their efforts. 

By embedding Fermi liquid theory in the framework of RG and 1/N we not only automate 
the process that Landau had to finesse with his genius, we also prepare ourselves better to 
study variants of the problem he attacked, such as the problem with impurities which is a 
big field (Lee and Ramakrishnan 1985). 

We then studied generic nonrotationally invariant problems. We found that we were 
led to F and V which no longer depended on just the difference between their arguments. 
The flow equations for V were derived. By going to a Fermi surface with no time-reversal 
invariance, we could eliminate the BCS instability and have a real Fermi liquid theory at 
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T — 0. The new feature here was the changing shape of the Fermi surface as we renormahzed. 
Moving on to the case of a nested surface in o? = 2 , we found that a new couphng 

— * 

W corresponding to momentum transfer Qn — {n n) survived at tree level. At one loop 
this coupling began to flow. For the case of the perturbation corresponding to the nearest 
neighbor interaction (truncated to modes within the cut-off) , we saw the CDW instability. 
It was pointed out that the points (±7r/2 ±7r/2) where holes seemed to appear upon doping 
the half-filled system are exceptional from the point of view of the RG: the are forward 
scattering amplitudes that flow due to BCS and CDW diagrams. Several problems were 
proposed for further study. 

Finally we discussed the possibility of singular Landau parameters and non- fermi liquids. 
This was surely a possibility at strong coupling, but did not seem to happen at weak coupling 
within the scheme we were employing. 

Although we used the modern Kadanoff-Wilson approach to renormalization we fre- 
quently made contact with the old field theoretic scheme for computing the /9-functions. 
The main difference between the two schemes was that in the former, all loop momenta were 
at the cut-off, while in the latter, one was at the cut-off and the rest at or below it. In section 
II it was pointed out that in the case of the scalar field, the difference did not show up in 
the flow of the marginal couplings at one loop, since we could set all external momenta to 
zero and momentum conservation meant that if one propagator was at the cut-off, so was 
the other. It was pointed out in the fermion problem even if we set all /c's to zero, there 
could still be large momentum transfers of order Kp- Yet we saw repeatedly that the two 
approaches gave the same one loop flow for marginal couplings. Let us recall why. First, 
F never flowed in either scheme since it got its contribution at the Fermi surface and did 
not know about A in either scheme. As for V, when we set the incoming momenta on the 
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Fermi surface, at opposite angles, the total momentum vanished, just as in the scalar field 
theory, and the two propagators in the loop were equal and opposite and at the cut-ofi^ in 
the BCS diagram. The other two diagrams did have large momentum transfers. They were 
however suppressed in both schemes by kinematics, but not equally: the suppression was 
by dK/Kp\n the modern scheme and by K/Kp in the field theory scheme. This difference 
did not matter at the fixed point since both factors vanished. Finally W fiowed due to the 
ZS diagram. Even though the momentum transfer was large {Qn ) it was such that if one 
propagator was at the cut-off so was the other due to the condition E{K) = —E{K + Qn)- 
The other two diagrams were suppressed by phase space. 

The field theory scheme was the best choice for studying screening and the Kohn- 
Luttinger effect discussed in the Appendix. 

There are many possible extensions. Inclusion of spin will produce new effects such as 
spin-density waves, but no new formalism is required to incorporate it. Inclusion of bosons 
is very interesting, but not discussed here due to lack of space. Ye and Sachdev (1991) have 
used the RG ideas espoused here to study a boson-fermion system describing the metal - 
superconductor transition. Polchinski discussed phonons in his TASI article (1992) and in a 
more recent preprint (1993). 

We can study disordered systems using the replica trick. Finally we can go to finite T. 
The fixed point discussed here will describe the crossover to T = oo in the early stages. 

To conclude, the analysis in this paper has shown that the RG, which has proven its worth 
in critical phenomena, chaos etc., is just as effective in helping us understand interacting 
fermions. Conversely, the Fermi systems, with their novel phase space and fixed point 
structure, offer us a far from ordinary manifestation of the RG at work. 

With various phenomena such as Fermi liquid theory , BCS , CDW and SDW instabilities. 
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screening, and the Kohn-Luttinger phenomena all under the auspices of the RG, we have a 
better chance of solving extensions and generalizations of these problems. 
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A APPENDIX 



Here we will use the field theory approach to study two phenomena: screening of the Coulomb 
potential and the Kohn-Luttinger effect. A smooth cut-off will be used on loop momenta. 
Coulomb Screening 



We know that the instantaneous Coulomb interaction, before any mode elimination, 
is given by ^"^q^^ . The extra Kj, comes from the way our fields are normalized. The 
corresponding bare vertex is 



y(4321) = ATre'K'p 



(463) 



Since we will be focusing on Ki, we will drop the second term. 

It is generally agreed that the Coulomb potential gets screened. Should we be using a 
screened version here? No! Screening comes from organizing diagrammatic perturbation 
theory in a certain way by first summing a class of (RPA) diagrams. In RG the organization 
is different. Since the Coulomb potential or bare coupling in the action is unambiguously 
known before any mode elimination, we must see what it evolves into as we carry out the 
RGT. In the field theory approach, using a smooth cut-off e""'*^! where a = 1/A, we have to 
second order : 

W) = ViQa) + V^iQa) j (2.)^ - EiKM. - E^K')) ^'''^ 

where K' = K + Q, and where k and k' = \K'\ —Kp run from — oo to oo. We have not shown 
the ZS' and BCS diagrams since they do not dominate as the ZS diagram does at small Q. 
Our plan is to find the /3- function by setting the a-derivative of V to zero. But first let us 
evaluate the former more explicitly. Doing the uj integral we obtain two equal contributions 
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from processes where hole gets promoted to a particle and vice versa and end up with 

We can now do the k and z integrals. In doing so we replace K by Kp wherever appropriate 

and use the fact that k, k' and Q are all much smaller than Kp. For example we set 

k' = \K + Q\-Kp^k + ^^^^^ (466) 

and so on. The result is 

where in the last equation we have used a simpler function with the same limits at small 
and large values to facilitate the analysis. The /3-function is now computed to be 

jy ^ _ V^m aQ 
da KpT:^{l + aQf ^ ' 

Notice the flow is strongest at aQ ~ 1 i.e., Q ~ A, which makes sense in the sharp cut-off: 
we need a minimum momentum Q ~ 2A to knock a particle from the shell at —A to the 
shell at A. When Q is too different from this range, the flow is essentially nil. (Had we used 
sharp cut-off the (3 function would have had a string of 9 functions, which is why we do not.) 
Integrating the flow from a — Qioa — a 

V{QO) Kftt'^ l+aQ 

where V{QO) is the input potential before mode elimination. The final answer, in terms of 



A is that 



47rp2 7^2 
^ Q+A w 



^""^ (472) 



151 



where G is the inverse of the Thomas-Fermi screening length. 

This is the bare potential that goes into the action when the cut-off is A. To regain the 
potential before mode elimination, we must set A = oo. You may ask how we can have a 
A > Kp. The point is that A is defined as the inverse of the a in Eqn.(|46^) and not as the 



real cut-off for the integral. Indeed with a = 0, the integral is still finite because it is limited 
by Q: for very small Q, we cannot knock states far below Kp to states above it and only 
such processes contribute to the u;-integral. 

Now screening refers not to the bare coupling in the action but to the full physical four 
point function. It is however clear that at fixed Q if we send a to oo we kill all loops and 
the bare coupling itself gives the full answer, which is 



r(Q) = T^(Q,oo) 

(474) 



Q2 + e2 

Let us examine what we have above. If we fix Q at some value and lower A, this is 
what happens to the bare charge in the action. First, at A = oo, we have the unscreened 
potential A-Ke^Kp/ Q'^ . It has this form as long as A >> Q^/Q. When Q « Q^/Q, it goes 
as Ane'^KpA/lQ'^Q). Finally when A << Q, we get the screened form in Eqn.( |47^ . 



We must be clear about what is done here: we tried to understand how screening takes 
place in the RG scheme as we eliminate modes. Since we assumed the fermion propagator 
had the Fermi liquid theory form throughout, we haven't really verified that non- Fermi 
liquid theory is ruled out. In other words, we must study the evolution of the fermion 
propagator as modified by the Coulomb potential as we go along. This has not been done; 
we have assumed Fermi liquid behavior. 
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The Kohn-Luttinger Effect 



Years ago Kohn and Luttinger pointed out that in principle any system will face the BCS 
instability at low temperature, even if the initial coupling is repulsive. Let us recall their 
argument with no reference to the RG. Consider the BCS amphtude to one loop as shown 
in Figure 18 : 

r{-K3 Ks - Ki Ki) = r{Q = - Ks) = V{Q) + BCS + ZS + ZS' (475) 
Let us compute the coefficients of the Legendre expansion 

r, = l\{z)r{z)dz (476) 

— * — * 

where z — Q,i ■ Qs and I will be odd due to the Pauli principle. 

The bare potential will make a contribution Vi, assumed to be positive. Since V{z) is 
assumed to be analytic in the interval — 1 < z < 1, we must have Vi ~ e~' in order that 
the infinite sum over polynomials and all derivatives converge. Now look at the one loop 
corrections. They are nominally smaller, being of second order. However their dependence 
on I is very interesting. This is because the ZS and ZS' graphs have singularities when 
Q = Ki — Ks — 2Kf or Q' — Ki — K4 — 2Kf respectively, which correspond to z = ^1. 
Let us focus on just ZS, since the Pauli principle will determine ZS' for us later on. Due to 
the singularity, the Legendre expansion coefficients fall as 

S^i ^ -f (477) 

where V(7r), which enters both vertices is the backward scattering amplitude. (Let us see 
why. At the left vertex 1 gets knocked into 3 = —1, so that we have a momentum transfer 
2Kf in the direction of 1. The loop momentum K must be nearly —1 getting knocked into 
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+1 if it is to lie on or near the Fermi surface and obey momentum conservation. A similar 
argument applies at the other vertex.) It follows that if we hold V fixed and look at large I, 
the second term, which is attractive (for odd I, which is all we have) will dominate over the 
first, which is falling exponentially. There is no question of hoping that V{Tr) — 0, since it 
is (up to an overall minus sign) the sum of all the VJ's, all assumed non- negative. 

We are trying to reproduce this in the RG language. The procedure will be just as in 
the screening calculation, except now we do not assume Q is much smaller than Kp, indeed 
it is nearly 2Kf- Now we get, upon doing the cu integral and setting all unimportant factors 
to unity but paying attention to the sign, 

POO p1 Q(y\p-oik -ak' 

m)^V(Q)-V'Ml ,,/_^,.±J-^ (478) 

where, once again « = the primed quantities refer to K' = K + Q. Let us define 

X = 2Kf - Q. (479) 

By drawing a sketch of the Fermi surface you may verify that since all the action is near 
Q = 2Kf, 

and that — l<2;<lforA;>a; while < z < 1 for k < x with z^ being the point where 
A;' = 0. Putting all this in 



r = 



V - V^n) 

f 

Jo 






X 2KQ 'Q-2K 
+ 2KQ 



2KQ 



2KFk 



-2ak 



(481) 



By subtracting the second integral from to oo, which does not affect singularity in 
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question, and shifting the origin and rescahng k, we get 



y + x 



(482) 



We will drop the e since it does not modify the singularity at a; = 0. Next we approximate 
as follows: 



X = 2Kf-Q 

= 2X^.(1 -sin ^13/2) 

~ -sin ^13/2) (1 + sin ^13/2) 

= Kf{1 + z) 



(483) 



Hereafter we will use Kp — 1. We now do an integration by parts, throw out the surface 
term and obtain 



r = y 



a Jo 



-ay 



1 



y + l + z 2y + l + z 



dy 



Now do the angular momentum transform using 



Qi{zo) 



1 /•! Pi{z)dz 



2 J-i Zq — z 



to obtain 



(484) 



(485) 



r = y + 



a Ji 



Qi{y) 



dy. 



Next we use the result that as Z — > 00, 



Qi{y) ^ V2/7rXe 



-X 



(486) 



(487) 
(488) 
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to obtain for large I and a 



r ^ V + 



aP Jo 



g-Q;X2/2/2 _ g-Q;X2/4/2 



dX 



= V + L. 



(489) 
(490) 



Note that 



L{a/l — >• oo) ~ 
L{l/a — > oo) ~ 



We fit L with a simpler function with the same limits: 



r 

R(x) 



-R{a/l') 



X 



1 + 



(491) 
(492) 



(493) 
(494) 



If we calculate the /3-function (including a factor of 2 due to the ZS' diagram) and also 
include the usual contribution from the BCS diagram we get the result quoted in the text: 



dt 



47r 115/2 [^7/4 + /-7/2] 



2 ■ 



(495) 
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C TABLES 



Table I Definitions of terms frequently used in connection with the RG 


Symbol 


Meaning 


RG 

Z 
A 

s and t 

Slow modes 0< 
Fast modes 0> 
/9-f unction 
S 

s 

Sa or Sb 
Critical system 
Critical exponents 
Critical surface 

S* 

Relevant variable 
Irrelevant variable 
Marginal variable 


Renormalization Group 

Classical partition function or Feynman's path integral for quantum problem. 

The cut-off, the maximum allowed value of momentum k. 

The parameter in the RG: A = Aq/s = Aqc"*, where Aq is fixed. 

Modes to be retained. 

Modes to be integrated out. 

The rate of change of couplings with t, the logarithm of the cut-off. 

The action or hamiltonian.The Boltzmann weight is e^. 

The space of all hamiltonians. Each axis is used to measure one parameter. 

Any two actions or hamiltonians describing two different systems. Points in S. 

A system tuned to be at a phase transition. Has power law correlations. 

Exponents for the above power laws. Are universal. 

The locus of all actions or hamiltonians that describe critical systems. 

The fixed point of the RG transformation in the space of hamiltonians. 

Any deviation from S* which gets amplified under the RG action. 

Any deviation from S* which gets renormalized to zero. 

Any deviation from S* which remains fixed under the action of the RG. 
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Tal)le II Fixed point couplings and flows: a snnniiary 


Item 


Comments and relationships 


< 0*(A;2)0(A;i) >=< 21 > 
G{k) = 

Wick's Theorem 

< 4321 > =< 42 >< 31 > + < 41 >< 32 > 

SS = - J,^ ^<P\k)r{k)m 
r{k) = ro + r2k^ + ■ ■ ■ 

-2^/a0*(4)0*(3)0(2)0(1M4321) 
ii(4321) = -uo + 0{k) 

RG action 

c 

Cumulant expansion 
Tree graphs 
Loop graphs 

Connected graphs 
Tree level RG 

RG flow to one loop 


Complex scalar field and its conjugate. 
Two point function. 

Propagator. < (t)*{k)^{k') >= {27r)^S\k - k')G{k). 
Gives N- point functions in terms of 2-point function 
in gaussian model. 

Four point function evaluated using Wick's theorem. 
Gaussian model action with cut-off A. 
Quadratic perturbation. 

Quadratic coupling function. Only ro is relevant; 

r2 is marginal and the rest are irrelevant. 

Quartic perturbation in schematic form. 

Quartic coupling function. m(4321) = m(3421) = m(4312). 

Only uq is marginal, higher Taylor coefiicents are irrrelevant. 

Reduce A by s by integrating out fast modes. 

Rewrite result in terms of 4>'{k') = s~^4>{k), where k' — sk. 

Field rescaling factors defined by (j)'{k') — C0(^)- 

Chosen to make gaussian action the fixed point. 

< >= e<^>+5(<^^>"<^>^)+--- 

Graphs with no closed loops. 

Graphs in which there are closed loops. 

One works to a given number of loops. 

Graphs in which there are no disjoint parts. 

Calculation with zero loops in Feynman diagrams. Reduced 

to ignoring fast modes and reexpressing the perturbation 

in terms of new momenta and fields. 

5 = 2ro + auo and ^ = -bu^, a,b>0. 
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Tal)le III Suiiiiiiary of tlie Foriiiioii (Jscillators 


Symbol 


Comments or definitions 


N = 
|0 > 

1 > 

Ho = Qo*^* 
/i 

P = l/T 
H^Ho-iiN 

/ O ^ Trne-l^<^"o-^^N) 

A 

"Hubbard model" 

<N> 

T(^(r)*t(o)) 

G{t) 
-G(O-) 

G(r) 


Fermion destruction and creation operators, 
{^t^ ^} = ^t^ + ^^t = 1 ^ ^2 ^ Q_ 

Number operator. N = orl. 
State with = 0. 
State with = 1. 
Oscillator hamiltonian. 
Chemical potential. 
Inverse temperature. 

Free energy operator. Also called hamiltonian. 

Grand partition function. 

Average of operator fl. 

Free energy. Defined by Z = e~^^^^'^K 

Ho = + ^N{N - 1) 

Average occupation. At /3 — >■ oo this becomes 

9{fi — f2o) for single oscillator, 

elu - fio) + 0{i^ -Qo-U) for "Hubbard model" . 

^(r)^(r)^t(o) - ^(-r)^t(o)^(r). 

< r(*(T)*t(o)) >. The Green's function. 

N 

foo f^l \ -iuirduj 
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Tal)le IV Foriiiioii Patli Integrals: Useful Relations 



ip and ip arc Grassmann numbers. They anticommute with each other and with 
fermion creation and destruction operators. Their differentials are also Grassmann numbers. 

|0 > -ijj]! > 
0|- < 

>= 4j\ip > 

< ^l^t 

< ^\tp >= 

jijdij = 1 

Jld^^ 

Jl(hp = 
J ipipdipdip = 1 = — / ipipdipdip 

J e-^.A//,,^,^- [H- di^idi/ji] = detM 

ijjii/jjijjki>i)^ = - =<ll xjk > - <Jk xjl > (Wick's theorem) 

I = J\iP ><^\e-^'^dHi^ 

< ^K)^(-2) >= 

< i^iuMu) >= 
27r5(0) = P 

^ iHz a„. ~ J-c 



PZ dfx J-oo 27r Vjj-Q.o+iJi 
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Tal)le V Spiiiless Feiiiiioiis in (1 = 1: Suiiiiiiary of syiiilxjls and fornmlae. 
CDW Charge density wave. 

H = -^j:n^Kn + + h.c. + Uo T.ni'^HnMn) - l/2)(V't(^ + + 1) _ 1/2) 

rij = "^p'^ jipj Particle number at site j. 

11 = Uq The chemical potential that ensures half- filling for repushon Uq. 
Tij = \ + |(— 1)-'A The mean field ansatz. H appears as follows in terms of A: 
H = -\Y.n i^Kn + l)V'(n) + h.c + Uo[\ Ej - A Eji-iYrij] + Uo [Ej ■ rij :: rij+i : 
: Uj : = charge density at site j with mean value subtracted. 
UoJ2j '■ ^j+i '■ Part neglected in mean- field calculation. 
A CDW order parameter. 

E{k) = — cosX Dispersion relation for free fermions. 
Kp = 7r/2 The Fermi momentum at half filling. 
K Momentum measured from origin. 

k = \K\ — Kp Momentum measured relative to the Fermi surface (points). 
L and R: Names for left and right Fermi points K = 
Ho = - ^i'KK)i){K) cos K Free field hamiltonian. 
Hq = Ei I-A ^'4'i{k)'4'i{k)k Free field hamiltonian hnearized near Kp. 

Zq — I Y[i=L,RY[\k\<Adtpi{^f^)dtpi{^f^)^^° F*r6e fermion partition function. 

-^0 = Ei=L,K/-A f /~oo - k)i>^{u;k) Free fermion action. 

k' = sk u' = suj '^[{k'uj') = s''^l'^'4>i{kijj) RG transformation. 

= Ei /-A H SZo ^f^{kuj)'^i{ook)ilJi{u;k) Quadratic perturbation. 

Ii{k,uj) = yUoo + f^mk + fioii^ + ■ ■ ■ Taylor expansion of qudratic coupling. 
Only is relevant. The rest are marginal or irrelevant. 

SS4 = 2!iT /ft-o; '0(4)'0(3)'0(2)'0(1)'^(4, 3, 2, 1) Quartic perturbation in schematic form. 
lK.= lnU I-n /~oo [2vr?(iri + K,-Ks- K,)27r6{u, + - - u,) 

5: Delta function for momentum conservation modulo 27r. 

Umklapp process A process where momentum is changed by a multiple of 27r. 

u{4, 3, 2, 1) = Uo sm{^^^^) sm{^^^^) cos( ^^+-^^~^'^~^" ) Nearest neighbor model. 

Ui^i^i^i^{kl, uj'j) = Ui^i.^i2i-^{k'J s,u'Js) Renormalization of quartic coupling at tree level, i = L ,R. 



Uo = ulrlr = urlrl = -urllr = -ulrrl Symmetries of marginal couphng constant. 

dUQ 

dt 



^ = Flow at one loop. 
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Tal)le VI Phase space for various fixed points: Suiniriary of syiiil)ols and formulae. 
One particle phase space is as follows: 

d = 1 : The Fermi surface is a pair of points, called L(left) and R (right), 
phase space = EIl I-a f I-oo ^^ere k=\K\- K^. 

d — 2: Circular case. The Fermi surface is parametrized by angle 9. 
phase space = J^'' f /_\ f J^^ where k ^ \K\ - K^. 

d = 2 : Noncircular case. The Fermi surface is parametrized by an angle 6 if connected 
and by an additional label a = ±1 if it has two branches. 

phase space = J2a lo'' If '^(^) /-A ^ I^oo if' where J is the Jacobian on the Fermi surface 
and e measures the energy relative to the Fermi surface. 

d = 3: Spherical case. The Fermi surface is parametrized by 2: = cos^ and (j). 
phase space = J ^ j\§JZo't ^^^^^ = ^"^^^^ = ^■ 

In all cases the RG transformation reduces A (limit on k or e) by a factor s — e~*. 
The quartic interaction is written in schematic form as 

5S^ = /^V^(4)V^(3)V^(2)V^(1)«(4321) 

where the labels 1 to 4 stand for all the attributes of the fields: 
momentum, frequency and spin if included. 

The coupling obeys m(4321) = -m(3421) = -m(4312) = m(3412) due to Fermi statistics. 

stands for the integral over the above phase space for each of the four 
fields, restricted by delta functions that impose u) and momentum 
conservation (modulo 27r if appropriate). 
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Tal)le VII Fixed point couplings and flows: a snnnnary 



The marginal part of coupling u depends only on the angles that parametrize the Fermi surface 
Subscript "NN" denotes values in nearest neighbor model 



Fermi Surface 



Couplings marginal at tree Level and their flow at one loop 



A: d=2: 

Circular Fermi surface 



u{e, = 02, 9s = 9^, 92, 9^) = -u{9, = 9u 9^ = 9^, ^2, ^1) 
= F{9r,92) = F{9i- 92) = F{9i2) 
Fnn = Uo(l - cos6'i2) 



dt 



u{9^ = -93, ^3, 92 = -9i, 9i) = V{9s - 9i) = V{9n) 
'^^^ = -8^ r iVi9, - 9)V{9 - 9s) 
Vnn = Uo cos 9i3 

^--i,Vi- Jo" ie'^'V{9). (BCS instability if 1^ < 0) 



B: d=2: Noncircular 
Time reversal invariant. 
Asumed connected 



F same as in A, but F{9i, ^2) ^ F{9i - ^2) 
V same as in A but ¥{9^, 9i) ^ V{9s - 9i) 
= -8^ t f ^(^1; 0)V{9- 9,)J{9) 



dF 
dt 



0. 



C: d=2: Non circular 
No time-reversal 
invariance 



F is as in 5. ^ = 0. 

at 

No V. No BCS instability 



D: d=2: nested 
a labels branch 

r: hopping anisotropy 



F defined as in B and ^ = 

dt 

Fnn = t/o(sin2(^) + i(l - - ^Vr' - cos^ 9iVr'' - cos' ^2)) 

V defined as in S. Same fiow as in 5 

VNN{93a3 ^lai] = Uo[sm9i sin ^3 + ^Vr' - cos' 9Wr'' - cos' ^3] 

W[92a2; 9iai] = u[92 + tt, -a2; 9i + n, -ai, 920L2\ ^iQii] = u{9'^, d^, 9[a[; 92OL2, 6'iq;i) 

'^^'"1,^ "^'^^ = - / Ea [«2^2 «' W {9a 9^a^\ ^ 



W^iViv = -C/o(sin2(^) + i 



i?: d=3: Spherical 
f2j is direction on sphere. 
0i2;34 is angle between 
1 — 2 and 3 — 4 planes 



w(^4,f^3,^^2,^^l)ln,.0,=n4.Q3 =^(^l-^2;</>12;34) = f 

0) is Landau's F function. 
Fnn = Uo{l — z) cos (p. 

— * — * 

Vnn ~ UqZi3 — UqQ^i ■ Q3. 

^ = -cVi^ where c> and = / y(z)P,(z)dz 



= 0. 



167 



D FIGURE CAPTIONS 

Figure 1 Four types of diagrams that appear in the cumulant expansion to lowest order, 
corresponding to the 2^ = 16 choices for the legs on the vertex to be fast (F) or slow (S). 
Figl.(a) corresponds to all slow modes, denoted by S. This a tree graph with no integration 
over fast (F) modes. Fig. lb typifies graphs that vanish since they involve an odd number 
of fast lines. Figure Ic is the tadpole graph. Here two fast lines have been joined to form a 
loop. This term with two external slow lines will renormalize the quadratic term. Figure Id 
is a two loop graph. It comes from a vertex with all fast lines upon joining them in pairs. 
Its contribution is a constant as far as the slow modes are concerned. 

Figure 2 The three one loop graphs that renormalize the quartic couphng of the scalar field. 
The names ZS, ZS' and BCS are used to label the topology of the graphs and do not imply 
the corresponding phenomena (like BCS superconductivity). In part (a) lines 1 and 3 meet 
at a vertex, in (b) 1 and 4 meet while in the last lines 1 and 2 meet. All loop momenta 
lie in the shell being eliminated. The external momenta can be chosen to vanish if we want 
the renormalization of the marginal part of the quartic coupling. When this choice is made, 
both propagators have the same momentum. This is true for all three graphs. 
Figure 3 One loop flow of the quadratic (ro and quartic Uq couphngs in the scalar field in 
d — 4. Notice that although uo is irrelevant, a point on the uq axis does not flow to the 
gaussian fixed point at the origin. For this this to happen, we must tune the parameters to 
lie on the critical surface, shown by the arrow flowing into the origin. 

Figure 4 Graphical representation of the perturbative calculation of the type 1 fermionic 
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propagator in the toy Hubbard model. Figure 4a shows the two terms we kept in the analysis 
to one loop. Figure 4b shows the corrections that are graphically disconnected and ignorable 
since they get cancelled by the partition function that comes in the denominator of all aver- 
ages. Fig 4c denotes a correction which embellishes the type 2 propagator in the loop we did 
consider. The effect of such graphs is to turn the free propagator in the loop to the exact 
propagator for 2 which in turn means the integral over the loop equals the exact density 
A'2. The last part shows a two loop contribution that vanishes upon to u integration. We 
have not shown the iterates of the connected diagrams which are part of the geometric series 
which we assumed and summed in the text. 

Figure 5 The tadpole graph which renormalizes the fermion at one loop. It has no de- 
pendence on the deviation of the external momentum from Kp or uj We have used this 
freedom to set both these to zero on the external legs. The effect of this graph may be 
neutralized by a counter-term corresponding to a change in chemical potential. One may do 
this if one wants to preserve Kp. 

Figure 6 The one loop graphs for /3(it) for the spinless fermions md — 1. The loop momenta 
lie in the shell of width dh. being eliminated. The external frequencies being all zero, the loop 
frequencies are equal for ZS and ZS' graphs and equal and opposite for the BCS graph. The 
ZS graph does not contribute since both loop momenta are equal (the momentum transfer 
Q at the left vertex is 0) and lie a distance A above or below the Fermi surface and the uj in- 
tegration vanishes when the poles lie on the same half-plane. The ZS' graph has momentum 
transfer tt at the left and right ends. This changes the sign of the energy of the hne entering 
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left the vertex. The u integral is nonzero, the poles being on opposite half-planes. The BCS 
graph (c) also survives since the momenta loop momenta are equal and opposite (since the 
incoming momentum is zero) and this again makes the poles go to the opposite half-planes 
because the hues have opposite frequencies. The labels 1 ... 6 refer to the master Eqn.(312). 
Figure 7 The figure shows the regions of momentum space being integrated out in the d — 1 
spinless fermion problem. The thick lines stand for the slices of width dA.. They lie a distance 
A from the Fermi points L and R. In the ZS graph which has zero momentum transfer, both 
lines lie on the same slice and the a;-integral gives zero. In the ZS' graph, the momentum 
transfer tt, connects a and c (which have opposite energies) and h and d similarly related. In 
the BCS diagram the loop momenta are opposite and correspond to a, d or 6, c. 
Figure 8 The geometric construction for determining the allowed values of momenta. If 1 
and 2 add up to P, 3 and 4 are constrained as shown, if they are to add up to P and lie 
within the cut-off. Note that both the magnitude of 3 and its direction 6*3 are constrained 
to lie within ~ A/i^T^,' of 9i. Had we chosen a that terminated in intersection II, O^, would 
have been within KjKp of Qi- At the fixed point, the angels become equal pairwise. If the 
incoming momenta 1 and 2 are equal and opposite,the two shells coalesce and 3 and 4 are 
free to point in all directions as long as they are equal and opposite. 

Figure 9 The diagrams that renormahze the marginal quartic couplings for a d = 2 circular 
Fermi surface. All external frequencies are chosen to vanish and all external momenta are 
on the Fermi surface at angles given by ^2^. The marginal couplings obey (^3 = and 

— * — * — * — * 

Jl4 = Q2 {F couphng), or the above with 4 <-> 3 ( — F, by Fermi statistics) or Q2 = — f^i 
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— * — * 

and Q4 — — Jla, (coupling V ). All loop momenta must lie within a shell of width dA. In 
the first case (F), there is no momentum transfer at the left vertex {Q — 0) of the ZS graph. 
No matter what direction we choose for K, the other loop momentum K + Q = K will also 
lie in the shell. However the two propagators have the same energy and the a;-integral gives 
zero. In the other two graphs the loop angle is restricted to to lie in a region of width dA/A. 
Figure 10 shows this for the ZS' graph and Fig. 8 can be used to show this for the BCS graph 
if we replace the annuli of thickness A with shells of thickness dA. Thus none of the diagrams 
cause a flow of F at the fixed point. The flow of V, receives no contribution from the ZS 
and ZS' diagrams because there is no correlation between incoming and outgoing momenta 
and both Q and Q' will be of order Kp and the diagrams will be down by dA/A. The BCS 
diagram however does cause a flow of V exactly as in d = 1: the loop frequencies are equal 
and opposite but the loop energies are equal due to time-reversal invariance. 
Figure 10 Construction for determining the allowed values of loop momenta in ZS'.The 
requirement that the loop momenta come from the shell and differ by Q' forces them to lie 
in one of the eight intersection regions of width (dA)^. The diagram does cause a flow of F 
or V for this kinematical reason. (In the field theory approach, the suppression factor will 
be A/Kp instead of dA/ Kp and the diagram is once again unimportant at the fixed point.) 
Figure 11 The sunrise diagram contribution to self-energy. This has two loops and is sup- 
pressed by a factor A/Kp, or in the language of the 1/N expansion, by a factor 1/N. 
Figure 12 Flow of indices in the 1/A^ analysis for the problem with an F type interaction. 
The tree level graph is order 1/N. The ZS graph is of the same order since the sum over I 
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pays for the extra vertex. The other two graphs are order 1/N'^. 

If we consider interactions of the V type, the BCS diagram and its iterates will be favored. 
Figure 13 Role of the counter-term (shown by a cross) in The tadpole (order u) is 

cancelled by the next diagram with just a cross. The value of the cross is found by first 
evaluating the tadpole and then arranging the cancellation. The next one with two loops is 
cancelled exactly by the one following it, where the cross still stands for its order u value. 
The sunrise diagram (last) will require a change in the counter-term, but this will be of 
relative order 1/A^ and hence ignorable in the large N or small K/Kp limit. 
Figure 14 The compressibility at small expressed as a correlation function of two 

densities. It is assumed we have renormalized down to a very small cut-off A which is still 
larger than probe momentum. The p's stand for the external probes or fields coupling to the 
particle density. The first graph on the right hand side is what one gets in free-field theory. 
The rest of the graphs are corrections due to interactions. The bubble sum, i.e., the iterates 
of the ZS loop are all of the same order in 1/N . If we try to introduce any large scattering 
amplitudes (these are marginal in d = 3) into the graph, (say across one of the ZS bubbles), 
there will be a suppression by powers of A./Kf. Only the ZS bubbles survive the kinematical 
restrictions that all propagators lie within the narrow cut-off. 

Figure 15 Ward identity at work. The idea is show that the coupling of the charge density 
to a field A is unaffected by interactions in the limit of vanishing probe frequency and mo- 
menta. Figure 15a represents the renormalization of the quadratic coupling by the addition 
of the self-energy E(A;, u). (The loop momenta are all fast.)If we expand E in a power series 
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in ioi), the first derivative term will modify the coefficient of the iu term in the quadratic 
piece by E' = The field must then be rescaled to neutrahze this. The graph in Figure 
15b represents E', the (external) frequency derivative of that in Figure 15a. The derivative 
introduces an extra propagator. The corresponding diagram represents the change in the 
field rescaling. However the graph in Fig. 15c, which refiects the fact that the slow fields can 
couple indirectly to the external field A via fast modes denoted hy ipp, equals the graph in 
Figure 15b when external momentum and frequency of probe vanish. Consequently the field 
rescaling and modification of the couphng to the external field exactly cancel for these two 
graphs. 

Figure 16 A Fermi surface with time-reversal symmetry but no rotational invariance. The 
dark line may be taken to represent a very thin shell left after a lot of renormalization. The 
figure shows that once again the set of initial vectors 1 and 2 coincides with the set of final 
vectors as the cut- off goes to zero. All the results of the circular case hold except for one 
thing: F{6i, 62) and ^(^3, ^1) are no longer functions of the differences of their arguments. 
The only instability is the BCS instability which only requires time-reversal invariance. In- 
deed the construction in the figure assumes this symmetry: the solutions it gives for K3 and 
K4 actually point from the surface to the center of the second ellipse instead of the other 
way around. However time-reversal invariance assures us that if we continue these vectors 
past the origin of the second ellipse, by an equal amount, we will hit the surface, so that 
and K4 as shown are acceptable solutions. In general, to find the final set using the initial 
set, we must draw the time-reversed version of the second surface and displace its center by 
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P relative to the first. This did not matter so far due to time-reversal invariance of the 

— * 

surfaces. If this condition were forgotten, one would erroneously conclude that when P — 0, 
the two surfaces will coincide and the Cooper pairs can roam over all angles, leading to the 
BCS instability. The correct construction would show in this case that the two surfaces, 
even with P — would intersect only in a few places of area ~ A^. 

Figure 17 A nested Fermi surface with Qn — (tt, vr) and hopping anisotropy r > 1. The 

filled states go from the origin to the dark lines, the two branched Fermi surface. Any point 
on Fermi surface (defined here by zero energy £ = 0) goes to another on the surface upon 
addition of Q which reverses the energy. Consequently any point just (above) below the 
filled sea goes to a point just above (below) the sea, leading to the failure of perturbation 

— * 

theory at second order when a perturbation of momentum Qjv is introduced and to the flow 

— * 

of the coupling W through the graph where the momentum transfer isQn- The thin lines are 
equal-energy contours at £ = ±A and stand for the infinitesimal shells being integrated out 
in the RG program. The point A at (7r/2 7r/2) is privileged. It lies on the Fermi surface for 

— * 

all r, scatters into minus itself under the addition of Qn- The coupling for A, —A — > A, — A 
is an F, a V and a W and fiows for more than one reason. It is conjectured that this is 
probably why holes are found at these points (and two more, obtained by refiecting on the 
y-axis) upon doping. 

Figure 18 The Kohn-Luttinger diagrams. F is the full BCS amphtude and V is the bare 
vertex. The ZS and ZS' diagrams have singularities when the momentum transfer equals 
2Kf whereupon X ~ — 1 in ZS scatters into X -|- Q ~ 1. Thus V — V{7r). 
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